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Shyamal Kumar Ray 
Department of Civil Engineering 
Indian Institute of Technology, Kanpur 

COMPUTER IMPLEMENT ATI OK OP OPSET - A SEIF CALIBRATING 

WATERSHED MODEL 

The present study is "basically related to the implemen- 
tation of OPSET, a self calibrating watershed model program, 
to the computer system IBM 7044 and IBM 370/155 and the appli- 
cation of this to an Indian watershed to find out the optimum 
set of some parameter values which influence the runoff m a 
watershed greatly. Test data are available to implement the 
program. For the Indian watershed the measurable watershed 
parameters and obher climatological data are collected from 
different places. 

Physical description of the watershed under study is 
given and the definitions and magnitudes of the measurable 
watershed parameter values are also presented here. While 
optimizing the parameter values the main objective was to 
synthesize the streamflow and hydrograph peaks as closely as 
possible with those of recorded values respectively. Synthesized 
and recorded values of flows are tabulated for comparison. 

Changes are made inside the program where necessary and 
are given in this study. Ultimately the optimized parameter 
values obtained from OPSET program can be used m the Kentucky 
Watershed Model to use in ungauged watersheds. 



CHAPTER I 


INTRODUCTION 

1.1 Hydrologic Cycle and Processes ; 

Hydrology is a branch of physical geography. It deals 
with the origin and distribution of the waters of the earth. 

The study of hydrology starts with the concept of hydro- 
logic cycle. ViTater evaporated from the oceans and the free 
water surface of the land mass is transported by moving air 
masses. Under proper conditions the vapour is condensed to 
form clouds, which, m turn, may result m precipitation. 

The precipitation which falls upon land is dispersed m seve- 
ral ways. Part of the water is returned to the atmosphere by 
the processes of evaporation and transpiration. Of the remain 
mg part, a portion passes over and through the surface soil* 
to the stream channels and other portion penetrates deeper 
inside the earth to become part of the earth’s ground water 
supply. Under the influence of gravity those surface waters 
and ground waters always move towards the lower elevation b”t 
a substantial quantity of these return to tho atmosphere 
through the processes of evaporation and transpiration. 

A review of the hydrologic cycle, particularly of 
those processes which are important in governing runoff rates, 
provides the best background for understanding the model m 
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this study namely OPSET, the basic logics of which are 
dependent on the Stanford Watershed Model. The processes are 
described briefly in the following paragraphs. 

Moisture held on vegetative surfaces is known as 
Interception Storage. This is dependent on typo and density 
of vegetative cover. This reduces the runoff but its effect 
is less m large storms compared to small storms. 

Procipitation not held by vegetative cover but hole m 
hollows and behind ridges on soil surface is known as Deprc- 
ssion Storage. This is also more effective m reducing run- 
off during small storms than during larger storms. Depression 
storage capacity is generally larger than interception storage 
capacity and is governed primarily by the roughness of the 
soil surface , human activity and watershed topography. Steep 
slopes reduce the depression storage capacity. In arid climate 
the less frequent runoff does not wash out the ndgos and 
depression storage is greater. 

A portion of the water in contact with the soil surface 
infiltrates into the soil which may return to the stream channel 
as Interflow or Baseflow or may be used by vegetation as ^va- 
potranspiration. Infiltration into tho ground depends generally 
on the soil permeability, surface characteristics, presence or 
absence of vegetation, etc. But freezing of the upper zono, 
rising of the underground water table and low permeability at 
moderate depths are some factors for reduced infiltration rate. 
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High ovapotranspiration rates will deplete the soil 
moisture between storms end consequently increase the infil- 
tration rate. Evaporation from soil surfaces and exposed 
water surfaces depends on the degree of exposure, air tempe- 
rature, water temperature, wind velocity, atmospheric pressure 
and percent solids m the water. Moisture loss through trans- 
piration depends upon the degree of vegetstive cover and other 
climatological ffetors. Urban development greatly mcieas^r 
the relative imperviousness of a watershed. 

Water which doesnot percolate to the water tabUe (to 
later appear as baseflow) map proceed toward tho stream in 
two ways ; (1) it mav travel over the land surface ( overland 
flow)* or 2 ) it may travel partially through and partially 
above ground (interflow). The time lag of overland flow is 
governed by tho slope, distance to the nearest collector 
stream and roughness of the soil surface. Overland flow is 
very rapid m impervious surfaces. 

If tho watcrtable is higher or if th.ro is an impervious 
layer at a shallow depth, the interflow increases. Bapid flow 
through the soil layers is soon m most pcrmeaolc soils. 
Quantity of interflow is dependent on tho local infiltration 
capacity. 

Overland flow, interflow, baseflow enter the channel 
system at different, scattered points m a 


watershed. The 
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channel inflows combine as they flow downstream and untimately 
roach the grugmg statioii m the stream. This channel flow 
rate depends on the slope, size, shape, and hydraulic rough- 
ness . 

A schematic diagram of the hydrologic cycle is shown 
m Pig. 1. 

1.2 Mathematical Modeling and Analysis ; 

full synthesis of the hydrologic cycle is practically 
impossible duo to vast complexity of the systems involved in 
the study, inadequacy of the knowledge now available and the 
knowledge likely to exist m the foreseeable future. The 
other way is to solve the practical technological problems by 
establishing workable relationships between measurable parame- 
ters m the hydrologic cycle. 

'Systems investigation’ is needed m hydrology to o«V 
lish quantitative relationships between precipitation an 1- ' m - 
off, which cm be used for ' reconstruction or prediction (R0±) 
of flood seouences and watershed yields. System mvestigaii 'C 
methods are (a) parametric and (b) stochastic. In the 
present work parametric study is done. 

Parametric hydrology is the development of relationships 
among physical parameters involved In hydrologic events and 
the use of these relationships to generate, or synthesize, 
non-recorded hydrologic sequences [l]. Among all the studies 
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made in parametric hydrology only the general non-linear 
analysis is independent of detailed knowledge of physical 
hydrology where emphasis is centred on topics of study. In 
the method of correlation analysis -various combinations of 
variables are tested to explore the significance of their 
effects m the hydrologic system. The combination, among 
those investigated, that yields a relationship most closely 
approximating the recorded output function m terms of the 
recorded input function and other arbitrary parameters is 
adopted as the bent prediction equation. 

In system analysis, the synthetic models m hydrology 
are dependent on the continuity equation of the form, 

I = Q +A8, where I = total inflow, Q = total outflow and 
A3 = change m internal storages for any given time interval. 
Any system like this, defined by the continuity of matter, is 
known as 'closed* system. 

Through the general non-linear analysis, which started 
a few years ago, one can find out the relationships between 
partial inputs and partial outputs to be established me eper 
dently, subject to some conditions of mathematical continuity 
and boundedness . It is not necessary to satisfy physical 
continuity conditions among the total input, total output and 
inner storage. When input and/or output represents only par- 
tial form of the total input or total output of the system, the 
system is called 'open'. Considering a hydrologic unit as an 
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open system one has obvious advantage m that it bypasses the 
need to evaluate component inputs and outputs. 

Hyd fologic system cannot be treated as analytic system 
because current mathematical knowledge is insufficient to eva- 
luate the nonlinear kernels of the functions of series when 
multiple inversion is not possible due to complex input. 

Physical simulation, for a. hydrologic system, is not 
possible due to limited data. Ultimately digital simulation 
helps m this regard. 

1 . 3 Simulation . 

fo simulate moans to duplicate the essence of the system 
or activity without actually attaining reality itself. It is 
used to circumvent the difficulties of duplication of environ- 
ment, of mathematical formulation, of lock of analytical 
solution techniques, or of experimental impossibilities . [ 2] . 

In hydrology simulation is used to correlate the ooiurvod 
hydrological behaviours to the factors which influence these 
and to predict the behaviour of watersheds whore observed cat-’ 
arc unavailable or incomplete. For this, largo quantities o 
physical data, on rainfall, evapotranspiration, stream! low, < nd 
other related variables are collected. Qualitative knowled ~ 
about the hydrological cycle is reasonably complete and inter- 
action between major components such as rainfall intensity and 
infiltration rates can bo logically explained. Still the 
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extension of qualitative knowledge to quantitative procedures 
is challenging. The land phase of the hydrological cycle 
contains complex interactions that are difficult to document 
from observed data, and the volume of data that must bo analysed 
is formidable. Streamflow simulation attempts to develop 
algorithms of hydrological processes for quantitative calcu- 
lations. In the present study hydrological reactions are 
assumed deterministic. 

Since wo are dealing with a physical system, a digitally 
simulated model deals with partial differential equations for 
rainfall, evapotrnnspiration, streamflow etc. This is too 
complex. Further, detailed and accurate data for soil mois- 
ture, overland flow etc. are not available m our country. 
Therefore, hydrologic basin models with lumped elements have 
boon used with different levels of details. 

1.3.1 Streamflow simulation models ,; 

There are many models for streamflow simulation, few of 
which are described here. Throe types of models are generally 
used these days. These arc (a) Mathematical Models, (b) Lot •na- 
tion Models and (c) Conceptual Models. 

In the study of hydrology with mathematical modeling, 
run-off process can bo analysed. Hydrometeorological charac- 
teristics which d efine the behaviour of catchments and th<- 
construction of mathematical models which approxime-toly 
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reflect the process of the runoff ar~ complementary. In 
modeling runoff some difficulties arc encountered such as 
lock of basic data combined with considerable noise level 
and uncertainty of some links in the hydrologic process 
nrmely the characteristics of water movement m the drainage 
basin, its depth and the circumstances of its emergence on 
to the surface, followed by routing through the channel. To 
overcome all those, some hypotheses m models arc used to 
match piodul with the natural data. The simplest of this kind 
is a model of snow-melt runoff on the surface of a mountain. 

In Retention Model, a catchment is regarded as an inter- 
connected system of catchment units and channel units. Each 
catchment unit consists of a retention storage and a runoff 
storage. The retention s corage does not meet the stream How 
and is composed of capillary soil moisture, pr-cipitr tmr > 
intercepted by vegetation, and some parts of the tlcprossi • ■ , 
ground water and channel storages. Thu runoff stora.gc consists 
of the water that is likely to become streamiflow, including 
over land flow, ’free* water moving through soil or litter, 
ground water and storage m minor channels. Informations of 
leakage, channel transmission losses are taken into a.ccount, 
if available. This model was originally suggested for analys- 
ing the effects of vegetation and soils on runoff volumes but 
it seems just as suitable for short term forecasting of floods 
due to rainfall [5]. 
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The objective of a Conceptual Catchment Model is to 
simulate as accurately as possible the streamflow hydrograph 
under natural conditions and also to evaluate the modifications 
to natural flows resulting from man-made physical changes or 
management practices. Ideally, a conceptual catchment model 
should satisfy the following four conditions for maximum use- 
funless m river f orecas ting; 

(1) It should represent the significant hydrological 
processes m a rational manner; 

(2) It should contain the minimum number of parameters 
required for adequate representation of the processes ; 

(3) The values of tho parameters should be measurable or 

be significantly correlated to easily measurable characteristics 

and 

(4) Tho model should contain procedures for updating as 
new information becomes available, m a rational if not opti- 
mal manner. 

There is still much to be learnt about the various 
phases of tho hydrologic cycle and research is being errriee 
out to increase tho knowledge of them. 

1.3.2 Conceptual models ; 

Tho basic philosophy of some of the well known conceptual 
models are presented below. 
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1. Stanford Watershed Model (SWM) ; 

'lh.is is the "best known of the current conceptual models. 
A brief description only is given here and tho fuller descrip- 
tion is given in Chapter 2. The SWM simulates continuously 
the movement of water as it moves over, into and through the 
soil. It uses as input rainfall, snowmelt or both. Interception 
ana storage losses must he satisfied initially. A portion of 
tho runoff comes over tho impervious surface and joins tho 
streamflow and the rest passes over the pervious surfaco and 
is subject to tho soil moisture conditions. Infiltration 
capacities vary from watershed to watershed. Detention storrgu, 
overland flow and interflow are due to the accumulation of wru.r 
on the upper zone. Water moves from upper zone to lower Zone 
and ultimately joins the ground water storage. Water stored 
in the soil moisture zone is called lower zone storage, lake 
evaporation from class A pan records is assumed to be poten- 
tial evapotranspiration. Evapotranspiration is assumed to 
take place from interception storage, upper zone, lower zone 
and ground water storage but at different rates. Ground water 
runoff, mterfow and surface runoff arc routed soperatoly and 
combined to produce outflow hydrograph. Hydrological processes 
arc represented mathematically and the parameters involved m 
these mathematical quantities are obtained by the digital 
computer. 
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2. The Utah State Simulation Model ; 

This simulation model has been developed by the Utah 
State University. The computer needed for its solution is an 
analog computer. Its objectives are (1) the development of 
improved relations for describing the various hydrological 
processes and inter-connecting link between these processes, 
and (11) the development of an analog computer having a high 
degree of flexibility and capability for the solution of hydro- 
logical and related problems. 

The hydrologic balance of this model is given by, 

P = I + A M 0 + SRO + A G 0 + CTO + ET + AS 

b S o 

where, 

P = Precipitation, 

I = Intcrccution loss, 

A M g = Change m soil moisture storage, 

SRO = Surface run-off, 

= Change m ground water storage, 

GUO = Ground water run-off, 

BT = Ovapotranspirrtion, and 

AS = Change m surface storage, 
s 

Interception loss is a function of type and density of 
cover and precipitation. SRO is the excess of precipitation 
ov.r the infiltration rate. Soil moisture storage is a function 
of infiltration, evapotranspiration, interflow and percolation 
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to ground water storage. The relation of actual to potential 
evapotranspirstion rote is assumed to be a function of soil 
moisture deficiency. Percolation of ground water is dependent 
on soil moisture storage. Ground water runoff is a function 
of ground water storage. The increments of surface runoff, 
inter flow and ground water are routed seporately through 
different amount of storage. 

3 . The Agricultural Research idorvice Model : 

This model, responsible for agricultural research nn 
the U.S., is based on infiltration theory approach. In this 
model first the depression storage on the surfa.ee must be satis- 
fied. An assumption is made that some impeding stratum of the 
soil acts as a control on saturated flow. Therefore, the infil- 
tration capacity (f) of a surface soil approaches the low eo±. r - 
tanc seepage rate (f ) of the impeding stratum as the oveny. ng 

V 

storage (S) is exhausted by infiltration volumes (F ) . T ad- 
equation is as follows s 

f = a (S - F) 1 * 4 + f 

where, a = percentage of area occupied by plant crowns, 
because plant roots connect largo pores and provide continuous 
channel. 

For this model, soil data arc needed which permit deter- 
minations of soil porosities and occurrence of impeding strata. 
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Soils fire grouped in accordance with their constant rate 
of infiltration after prolonged wetting. This constant rate 
is assumed as the seepage rate through the impeding strata. 
Porosities are classified as (i) pores (h), free water drain- 
able by gravity, and (ii) pores, available water capacity (ARC) , 
drainable by evapotranspiration. Recover of available sto- 
rage, by this means, is ft the rate of seepage (f ) to the 
extent oi freely drainable porosity, further recovery of 
available storage capacity between the days of rain is at the 
rate of ovapotranspiration (PT) to the extent of plant avaiJ able 
water (AWC). Different values of a’s, plant root coefficients, 
arc prepared relative to various land use and treatments. Poten- 
tial ET is estimated from pan- evaporation with a seasonal con- 
sumptive use coefficient applied to obtain the actual ovapo- 
t rams pi rat ion. Water in excess of infiltration and surfa.ee 
depression storage is assumed available for run-off and is 
routed by successive routing techniques to reproduce tho out- 
flow hydrogra.ph. 

4. The Columbia River Basin Mode l; 

The input to this model is weighted rainfall plus 
computed snowmelt on a daily basis. This daily input is dis- 
tributed into specific time period (3,6 or 12 hours). Tho 
total runoff is calculated by multiplying rainfall, snowmelu 
or both by a percent coefficient which is a function of soil 
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moisture index- The soil moisture index at the ond of a 
period is the soil moisture index at the beginning plus pre- 
cipitation minus total run-off minus evapotranspiration for 
the seme period. -Evapotranspiration on rainy day is less 
and is obtained by multiplying the daily value by a coeffi- 
cient. The percent of the total runoff contribution to base 
flow is computed as s function of a base flow infiltration 
index. The remaining runoff is separated into surface and 
subsurface components as a function of the total input rate. 
Components of surface , subsurface and brseflow runoff are 
routed separately, using a successive reservoir routing tech- 
nique. It is also possible m this model to introduce rule 
curves for operation of reservoirs and/or to specify release 
schedules for all or selected reservoirs. 

Studying all the above models it is observed that the 
Stanford Watershed Model is the best suited m this study, 
considering the availability of data. The Kentucky version 
of the Stanford Watershed Model known as Kentucky Watershed 
Model (KWM) was available with us for analysis. 

KWM has a number of parameters. Estimation of those 
parameters is vory time consuming. 'A self calibrating model' 
known as OPSUT is available for finding the best fit of vain., 
for the parameters so that simultated flow matches the rocorood 
flow. This model uses hourly rainfall data, daily st re onflow 
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data, and yearly, daily or 10-daily average evaporation data 
besides several other types of data. 

In India, there arc nearly 400 self recording rain gauge 
stations. The data from these stations are yet to be used m 
detailed hydrologic modeling of basins. It seems possible to 
use KWM m order to model basins with hourly rainfall date. 
Hence, it is proposed to implement xhe OPSE'T program developed 
by the University of Kentucky and available with us and use 
it, if possible, to model a basin for which data are available. 

1 . 4 Statement of the Problem % 

The objectives of the study are 

(1) to implement the OPSET programme for estimation of 
parameters of the KWM using tho test data, given in the original 
report [4], and 

(2) to estimatotho parameters for a basin m India with ciu 
requisite data. 

1.5 Scope of the Study ; 

(1) The implementation of OPSET on IBM 7044 and IBM 370/155 
computer systems is done using the set of available tost data. 

(2) Estimation of parameter values for a basin m India is 
done on the basis of available data. 

Within the limitations of time, computer facilities, 
money and data availability, parameters are determined only 
for 2 years for the Indian basin. 



CHAPTER II 


DETAILS OP OFSET PROGRAM 

2.1 Kentucky Watershed Model ; 

2.1.1 General Description s 

The Kentucky Watershed Model (EH) is the translated , 
revised and expanded version of the Stanford Watershed Model. 
Originally, the Stanford Watershed Model -was written m a 
digital computer language Sb 3 ALGOL. James [5,6,7] translated 
it into FORTRAN IV and modified it. This version is known as 
KWM. As the “basic logic of SVM and KWM are same, the des- 
cription of Stanford Watershed Model is given here. 

SWM is a systematic mathematical representation of the 
hydrologic processes m tno hydrologic cycle. Mathematical 
formulations m the model are empirical for some of the 
processes . 

In this model the major input data are evapotranspiration 
and rainfall. If snow is a significant part of precipitation, 
then snowfall data should also he given as input data. In^ “ t ’’' 
addition, several other input parameters and control arrays 
must he selected based on Watershed behaviour, characteristics 
and data availability. The details of input parameters ar>„ 
given by James [5]. 



18 


Output consists of synthetic streamflow, overland 
flow, interflow, baseflow, stream evapotranspiration (not 
and potential) and ground water storage. Those hydrologic 
quantities can bo obtained daily, monthly or annually as 
appropriate for uso by the researcher. 

2.1.2 Plow Chart of KWMs 

A schematic flow chart of KWM is given m Pig. 2. In 
the above mentioned figure the boxes represent classifications 
of moisture storage. Arrows art the processes through which 
moisture moves from one type of storage to another type of 
storage. The vertical straight line from precipitation to 
subsurface flow with arrow downwards represents the flow of 
water through the porous media due to gravity. When moisture 
cannot move downward due to some restrictions, such as low 
permeability, moisture starts accumulating in the higher 
boxes. Water moves laterally and joins the stream when 
storage capacities of these boxes arc exceeded. Ultimately, 
it goes out of the watershed. Vertical straight line with 
arrow upward indicates the withdrawal of water from all 
levels of storages due to radiation energy, and moisture 
needs of vegetation. 

The path taken by incoming moisture is determined by 
the antecedent moisture storage, the magnitude of the assigned 
parameter values, the entry rsit o and the season of the yQci. 
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storage 


Inter flow 
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The quantities appearing along the flow lines are the 
parameters which control tho flow of moisture along those 
paths. 

2.1.3 Input Data to KtfM ; 

The input data required by SWM or Fortran EWM arc 
divided into 6 groups as follows; 

1. Data to specify the desired program options and 
required specific output. 

2. Data to initialize the watershed soil moisture storage 
condition. 

3. Data to establish climatological events. 

4. Time-area histogram for watershed understudy. 

5. Data to assign values to the watershed parameters, 

6. Recorded streamflov to compare with the synthesized 
flow. 

2 . 2 OPSE'J s 

2.2.1 Definition of QPSET and Purpose of KWM ; 

OPSET is a self calibrating model. It is called OP SET 
because of its objective to estimate tho OPtimum SET of para.- 
meter values. Input parameters arc those which can bo directly 
measured by the user. Other parameters are repeatedly adjusted 
by OPSB'i' to match tho synthetic hydrograph with observed hycro- 
graph until no further adjustment is possible. Optimum set o" 
parameters are obtained from OPSBT for different water yoai 
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for a watershed and then averaged. They can he used in KWM. 
iCWM, in order to synthesize streamflows, uses these parameter 
values with much larger number of control options. A schematic 
diagram of parametric optimization procedure followed by 
OPSET program is shown in Pig. 3- 

2.2.2 Parameters Estimated, by OPSET ; 

There are altogether 13 parameters to be determined by 
OPSET which are classified as follows. These parameters are 
underlined in Pig. 2. 

Recession Constants ; 

L. IPRC - Interflow Recession Constant 

2 . BPEC - Baseflow Recession Constant 



(a) Runoff Volume Parameters: 


1ZC - lower Zone Storage Capacity, 

BMIR - Basic Maximum Infiltration Rate Within Watershed , 
SUZC - Seasonal Upper Zone Storage Capacity Factor, 

ETLP - Evapotranspiration loss Factor, 

BUZC - Basic Upper Zone Storage Capacity Factor, 

SIAC - Seasonal Infiltration Adjustment Factor. 


(b) Surface Volume Parameters: 

BIVF ~ Basic Interflow Volume Factor 
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Channel Routing Parameters ; 

1. ICTRI - Number of Current Time Routing Increments, 

2. CSRX - Channel Storage Routing Index, 

3. PSRX - Flood Plain Storage Routing Index, 

4. CHCAP - Channel Capacity - Indexed to Basin Outlet, 

2.2,3 Physical Significance of Parameters ; 

Over a real watershed, the hydrologic processes are 
continually going on at rates varying with time and loca.tion. 
The movement of moisture at any point is m response to acting 
forces (gravity) counteracted by the resistance to moisture 
movement along a given flow path. The resultant moisture 
moves m low resistance zone. Physical factors governing 
the saturated and unsaturated flows include the size and 
shape of the particles, the porosity, orientation, and the 
moisture content of the area, viscosity and surface tension of 
water and information on how these factors vary over time and 
space. In this model lumped estimation process has "been used. 
Option is there m SWM to divide the whole area into a number 
of segments. The estimate should also be lumped over a fi'-u te 
interval of time (15 mms . m SWM). They also can be lumped to 
a group of related hydrologic processes such as interception 
storage is added to depression storage m the SWM. Baseflow 
recession constants represent the cumulative effect of 
noisturc movement along a large number of routes, lumped 
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estimates imply that watershed parameters somehow aggregate 
the effects over space and time. Best values are those which 
give the synthesized flow close to the recorded flow, hut 
this depends on the data measurement error, size of the area 
and length of time represented etc. One point is very 
important that the parameters which are best m 15 minutes 
time increments will not be the best m 60 minutes time 
increments [8]. 

2.2.4 Parameter Estimation Criteria : 

There arc two criteria to choose the parameters which 
are to bo estimated by the OPSET program. First, tho parameters 
should be difficult to measure directly. Parameters like drain- 
age are? should not be estimated by matching synthetic flow 
tfith the recorded flow as it is easy to measure this to a satis- 
factory degree of accuracy. Second, the simulated flow sequences 
3hould be sensitive to tho parameters. Considering the above 
3wo criteria the thirteen parameters mentioned before 
fere Ghosen. 

1.3 Parameter Estimation : 

The procedure followed by OPSET program to estimate the 
>aramctcrs is outlined below. 

'•3.1 Recession Constants (BFRC and IFRO) ^ 

Traditional graphica.1 techniques to find out BERC (base 
low recession constant) and IFRC (inter flow recession const- 1 t) 
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are time consuming, The recession sequence is selected by 
the subroutine RECESS. Maximum length, of each sequence is 
50 days and for each station year upto 20 flow sequences 
are selected. Minimum number of days needed to estimate a 
single recession constant is 2 and for 2 constants 4 days 
are needed. Longer sequences give better results. As reces- 
sions from very smaill flow do not give good results a, criterion 
is used that the second day flow should be either greater than 
10 cfs or greater than 0.4 x AREA (where AREA is the watershed 
area m sq. miles). Second day is the day after peak or the 
first day whose flow is actually used to find recession 
constants. The recession sequence is terminated by a flow 
rise exceeding 0.1 x SQRT (AREA) cfs. After finding out the 
recession sequences through the subroutine REOESS the sub- 
routine SEI2RC and SET1RC arc used to find the values of IER0 
and BPRC. 

Subroutine SET2RC is called first to estimate the two 
recession constants IPRC and BPRC, If the sequence starts 
with a relatively low flow and flat recession then it is 
assumed that only base flow is present and subroutine SET1RC 
is used. Here BPRC is found. If the value of BPRC does not 
come up between 0.6 and 1.2 the entire sequence is discarded. 

2.3.2 Runoff Volume Parameters (LZC. BMIR. SUZC, EPLF. BUZC.SIAC) 

In order of decreasing sensitivity these parameters are 
DZC, BMIR, SUZC , ETLF, BUZC and SIAC . The subroutine SETFYP 
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adjusts the values of five flow volume parameters LZC, SUZC, 
BT1P, BUZC arid SIAC during the process of estimating the 
best set of values for the six flow volume parameters. In 
order to see how and to what degree each of these parameters 
affects runoff, a sentivity study was made with the data for 
Elkhorn Creek, Frankfort, Kentucky, USA, for the year 1964. 
First of all a best set of parameter values were selected by 
trial and error. Then for each of six parameters 2 computer 
runs were made by varying that parameter keeping others 
constants. Each parameter had effect on simulated flow and 
differed with the changes of typo of flow and changes m time 
of the year. BMIR is adjusted through the subroutine SETBMI 
and other five through SETFYP. These two subroutines function 
as unit and adjusts the six volume po„rameters simultaneously . 

First of all, trial values are taken for these poraao- 
tors and a year of streamflow is simulated. Then each parame- 
ter is adjusted by its adjustment rule. The new sot of six 
values is used to simulate another year of streamflows, rnd 
this process continues until the simulated flows have smaller 
SSQM as computed by subroutine SETFDI. Recorded and simulated 
flows cannot match perfectly due to data and modeling diffi- 


culties . 
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Runoff volume parameters are defined as follows ? 


1. BUZO 5 

lb is an index for estimating the capacity of the soil 
surface (uppor zone) to store water m interception and depres- 
sion storage. Wetter the underlying layer more will be the 
amount of water stored m this layer. 

2. SUZCs 


It is an index for estimating soil surface moisture 
storage capacity. Seasonal storage capacity changes due to 
summer vegetation and cultivation is taken care of by SITZC. 
BUZC and SUZC arc used to compute the upper soil zone nominal 
stora >e capacity (UZC) by the following equation. 


= SUZC (AJBX90) + BUZC 


where , 

ASX90 = Antecedent evaporation index (K = 0.90), 


1ZS 

IiZC 


= An index of the moisture content of the 
underlying soil, 


o = Natural logarithm base. 


5. LZC : 

It is the soil-moisturo storage capacity index which 
approximately equals the volume capacity of the soil to hold 
/rater. Decreasing LZC will be reducing the ability of soil 
bo hold water and thus increasing the synthesized flow and the 
reverse when the value of LZC increases. Delayed infiltration 
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or percolation happens from the upper zone to the ground water 
and lower zone storages when UZS/UZG exceeds the value of 
LZS/LZC (Ref. Figs. 4 and 5). 

4. ETIFs 

BTLF is an index used to estimate the maximum rate of 
evapotranspiration which could occur .within the watershed under 
current conditions of soil moisture content. This maximum rate 
is then used to estimate current actual evapotranspiration in 
the manner depicted in Fig. 6. More trees in a watershed me an s 

higher value of ETLF and less simulated flows (Ref. Fig. 6). 

• . 

5. SIAC i 

It is an evaporation-infiltration factor relating infil- 
tration rates to evaporation rates to account for more rapid 
infiltration during warmer period. Higher the value of SIAC 
less is the runoff. 

6. BMIR: 


It is the basic infiltration index used to control the 
rate of infiltration. Infiltration rates vary from point to 
point in a watershed and most runoff will he from points with 
smaller ratep, A cumulative frequency distribution is used in SWM, 

i ' 

varying linearly from zero to a maximum value (Ref. Fig. 7). 

The equations relating maximum capacity to infiltration 
rates, applicable at any particular point in time (CMIR), arc- 
given by, 

CMIR 

SIAM * Function (SIAQ) 

CIVM « BIVF x Function (XZS/bZC) 


Constant x SIAM x BMIR 
Function (LZS/LZC) 
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Increase m the value of BMIR reduces runoff during 
storm periods and later on increases runoff due to baseflow. 

2.J.3 Surface Volume Parameter (BIYF) ; 

It is anx.iindex controlling the time distribution and 
quantities of moisture entering the interflow. Increase m 
BIVF means more interflow by increasing the volume of 01 VII. 
BIVF enters the streamflow simulation through the equation, 

CIVi'l = BIVP x 2 (1ZS/LZ0) (Ref. Fig. 7) 

In order to chock negative interflow Cl VII must be > 1. 
It is determined by the subroutine SBTBIV. 

2.3.4 Channel Routing Parameters (CSRX. FSRX. BCTRI. CHCAl-’b 

The subroutine SB CEE P estimates single best values for 
NCi'RI and SRX for each of the 5 hydrographs (maximum 5) speci- 
fied m the input data. Hydrograph ECTRI values are averaged 
m it to have single best value of BCTRI. The storage routing 
uses the equation, 

°2 - I - SBX (I - 0l ) 

whore , 

C >2 = routed outflow at the end of the time mtervjl, 

I = average inflow during the time interval, 

0^ = outflow at tho beginning of the time interval, 

SRX = storage routing index, 

3RX for low flows is denoted by CSEX and SRX for flood flowsx 
ts denoted by PSRX. 




SETSRP estimates two channel routing parameters CSRZ, 


FSRX. When the synthesized streamflow is less than one-half of 
CHOAP, CSKX is used for routing. If the synthesized flow 
exceeds twice GHOAP, PSRX is used. When the synthesized flow 
is m "between these two SRX is interpolated from the equation, 

SBX = CSRX + (FSRX - CSRX) x ( ? 

where, 

Q is the synthesized streamflow. 

2.4 Parameters Directly Estimated from the Watershed - 
Characteristics s 

1. AREA : Total area of the watershed in sq, miles, 

2. FI HP : Fraction of watershed which is impervious and 

contributes its runoff directly into a stream. 

3. FWTR : Fraction of total watershed covered by water 

surfaces , 

4. VINTMRs Watershed interception volume storage capacity. 

It is dependent upon the typo and density of vogotnciv. 
cover. A table to choose VINTMR is given below: 

Table Is Interception Values for Various Types 
of Cover. [4] 


Watershed Gover 

VINTMR (m) 

Grass land 

0.10 

Moderate Forest cover 

0.15 

heavy Forest cover 

0.20 
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5. G-WETF : It is a factor which, when multiplied by the 

current rate of potential evaporation times 
the current ground water moisture storage 
gives an estimate of the current rate at which 
phreatophytes or swamp vegetation are drawing 
water from below the water table. For mosb of the 
be.sinsthis value is recommended to be zero [4]. 

6. SUBWF s Fraction of tho moisture entering ground water 

s borage which leaves the basin through sub- 
surface flow not measured by the stream gauge. 
G-enerally its value is zero. 

7. OFSS : Average slope (in ft /ft) of the overland flow 

surfaces perpendicular to the receiving channel. 
Mean value is taken from randomly selected 
group of points on a topographic map. 

To find out the discharge from the over land flow by 
Chszy-ilanning formula this OFSS is used m tho following 
equation: 

, = iUM (y 5/3) (s l/2) 

inhere , 

q = discharge in cfs/ft. 
y = depth of flow m ft. 

S = slope of tho surface m ft/ft (OFSS) 
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8 , 


9. 


10 . 


11 . 


12 . 


0FS1 : 


OFMN : 

OFMNI S j 

CHCAP : 


DIV : 


Mean overland flow length m ft. 

It is the average distance that surf a.ce runoff 
m the watershed travels before reaching a 
channel. It is estimated as the reciprocal 
of twice the drainage density in sq. miles. 
Manning's roughness coefficient for overland 
flow on soil surfaces. 

Manning's roughness coefficient for overland 
flow over impervious surfaces. 

Channel capacity - indexed to the outlet of basin 
CHCAP may be estimated from hydraulic analysis 
of the profile and cross-section of the stream 
channel. 

This is the daily flowdiversion m cfs. 

Diversions into the streams are considered as 
positive . 


2.5 Earlier Studios Using OPSEf ; 


The original work on the development of Watershed Model 
started m 1959 at Stanford University in USA. Its name was 
given as Stanford Watershed Model (SWM). This model was deve- 
loped through the versions of I, II, III and IV and was 
programmed for the digital computer m SUBAIG01 language, 

B.L. James translated the model III into FORTRAN (G-Level) 

and simultaneously incorporated most of the features of model IV. 

Thas model was named as Kentucky Watershed Model (KWM) . The 



University of Kentucky m the year 1970 developed a self 
calib rating model named as OPSET. Then this model was tested 
m 20 watersheds most of which are rural. Liou [8] and 
Ross [4-] gave full reports of the application of this OPSET 
models m the U.S.A. 



CHAPTER III 


IMPLEMENTATION OP OPSET PROGRAM 


3 . 1 Important Characteristics of the P^og-paTri ; 

Important characteristics of the OPSET program are 
described below 2 

3-1.1 Program Size and xfemory ; 

OPSET is a very lengthy program. In all, about 2400 
cards are needed to be punched (excluding data). Available 
core memory m IBM 7044 at III Kanpur is less than 32 K 
words (1 word = 6 bytes m IBM 7044 system). This memory 
is not sufficient for the OPSET program. On IBM 7044 two 
compilers are available viz. ( 1 ) WATFOR, ( 11 ) FORTRAN. 
WATFOR compiler gives an extensive diagnostic report but 
execution is slow. Due to memory problem, and even for 
debugging purpose, the whole program could not be run at a 
time. Mam program and each subroutine had to be fed sepax 
for correcting the syntax errors. This did not solve the pi 
pose. FORTRAN compiler ignores the undefined variables and 
some erroneous part of the program. 

All these difficulties are overcome through the 
IBM 370/155 computer system at IIT Madras. Here WATFIV com 
piler which has a memory upto 512K bytes, has been used to 
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debug the program. The total memory requirement with one set 
of data is 250 K (bytes) (1 word on IBM 370/155 system equals 
4 bytes). 

3.1.2 Computer Time ; 

In IBM 370/155 system through WATPIV compiler with on° 
year's data and all the options, computer time needed is 
15 minutes (compilation and execution). This implies that 
IBM 7044 will take more than one hour to do the same job* 
because IBM 370/155 compiles and executes respectively eight 
and four times foster compared to IBM 7044/1401 system at 
II T /Kanpur. 

3.1.3 TRIPS ; 

The program is divided into three portions. In execution 
of the first portion 6 variables are taken into account. In the 
execution of the second portion other variables are considered. 

In the final step all thS variables are taken togeth jr 
to give final synthesis. Execution ef «&eh of these portions 
is called a 'TRIP. First TRIP gives the optimized values of 
Lb C, BMIR, SUZC, ETLF, BUZC and SIAC. In addition to these 
values, two other parameters 3IVF and LZS are also obtained. 

The second TRIP starts working when the above values are optimized. 
In this TRIP the other six parameters CSEX, F,3Rk, 1TCTRI, CMC., f , 
IPRC, BFRC are optimized. Here the paramecer values of the 
first TRIP remains unaltered. In the third TRIP all the para- 
meter values taken together give the final run with the 


simulated flow. 
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3.1.4 language; 

i 

The OPSET program is written m FGKIEAJtf language. Ia this 
program one subroutine BEAD is used to read data m free for- 
mat. But this subroutine is not readily available with us. 

HEAD subroutine is therefore replaced by the formated EE AT 
statements . 

3.2 Modifications ; 

1 . Wat or Year . 

In the United States the water year is counted from the 
1st October of the one year to the 30th September of the next 
year. In this study it is taken from the month of Juno and 
is modified accordingly. 

2 . E vaporation Estimates ; 

In the test data of the OPSET program the evaporation 
is taken as the total evaporation m a year m inches and 
through the subroutine EVPDAY the daily values were calculi! d . 

In the present study the 10 daily average values of evapora- 
tion are used. 

3 • flo od Peaks ; 

The flood peaks data for tho watershed, used m the 
present worm, are not available. These values are, hence, inter- 
polated from the stream flow data,. 



3 • 3 Addition and Alterations : 


To take care of the nonavailability of some data pres- 
cribed for the program and to work with the program in a diffe- 
rent computer system, some necessary changes have been made. 
These are described below; 

1. All the evaporation data available m the study are m 
mm, but the OPSET program uses them in inches. Through the 
following statements the evaporation data m mm are changed to 
inches , 

990 DPET (JKRD) = 0.03937 *DPET (KRD) 

991 DPET (KRD) = 0.03937 * DPET (KRD) 

2. In the original program the following statements were 
used to read the recorded stream flow for one year. 

118 DAY = 274 

119 CALL BEAL (DRSE (DAY)) 

CALL DAYEXT (DAY,DPY) 

IE (DAY. EE. 274) GO TO 119 

It is necessary to change the above statements as in the present 
case the first day of a water year is the 152 nd day of the 
calender year and not the 274th day as stated m the program. 

The following statements are, therefore, introduced. 

118 DAY = 274 

EE AD 1000, (DTJM(J),J = 1, DPY) 

PRINT 1000, (DUM(J), J = 1, DPT) 

ISEQ = 1 
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119 LRSF (DAY) = DUM (ISEQ) 

ISEQ = ISEQ + 1 
1000 FORMAT (13F6.1) 

CALL DAY EXT (DAY, LIT) 

IP (DAY, EE. (396-DPY) ) GO TO 119 
Here one dummy array DUM (J) is introduced. June 1 data is 
read as October 1 data. But these two dates differ m a year 
by 122 days, so the last line m the above statement 274 
was replaced by 396-DRY. 

3. The following, statements are given m the mam program 

of OPSET. 


DO 126 KRD = 1, NSGRD 

CALL READ (ISGRB) 

126 CALL READ (DRSGP (ISGRD) ) 

C READ RECORDING- RAIN GAGE HOURLY TOTALS 

127 CALL READ (I7J3G, YEAR, MOUTH, LATE, CN) 

C PUNCH NO NUMBERS AFTER CN ON YEAR. EQ. 98 CARD 

IP (YEiiR.GE.98) GO TO 130 
HRF = 12* (CN-1) + 1 
HRL = 12* (CN-1) + 12 
DAY = MEDCY (MONTH) + DATE 
DO 128 HOUR = HRF, HRL 

128 CALL READ (DRHP (DAY, HOUR)) 

IF ( DRY. NS. 366 .OR. MONTH. NE. 2 .OR, DATE .NE. 29 )G0 TO 127 
he above statements are changed to the following form for 


easons mentioned above. 



DO 126 KRD = 1, NSGRD 

REID 500, ISGRD 

READ 103, 3L.SGP (ISGRD) 
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126 DRSGP (ISGRD) = 0. 03937 * DRSGP (ISGRD) 

C READ RECOILING ?J T • GIGE HOURLY TOTALS 

127 READ 500, IW3G, YEAR, MONTH, DATE, ON 

C PINCH NO NUMBERS APT LI; GIT OlT YEAR.EQ.98 CARD 
IE (fEAR.GD.98) GO TO 130 
HRF = 12* (CN~1) + 1 
HRL =12* (CF-1) + 12 
DAY = FiEDCY (MONTH) + DATE + 122 
IE (DAY.GT.DPY) DAY = DAY-DPY 
DO 2000 IJN =1,12 
IE (LAY.GT.MEDCY (IJN) ) GO TO 2000 
MONTH = IJJi - 1 
GO TO 2050 
2000 CONTINUE 
2050 CONTINUE 

READ 1005, (DRHP ( DAY, HOUR) , HOUR = HRE,HR1) 

DO 1291 HOUR = HRE, HRL 

1291 DRHP (DAY, HOUR) = 0. 03937 * DRHP (DAY, HOUR) 

1003 FORMAT (12P 5.2) 

1? ( DPY . NE . 3 6 6 . OR . MONTH . NE . 2 . OR. DATE . NE . 29 ) GO TO 127 
10 129 HOUR = HRE, URL 

Ocher .''Iterations are mainly the changes m the FORMAT statemenxs 
one 1 corree tions ol some minting errors in the OPSET program. 
These are not mentioned here. 
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3*4 Time Spent on Collecting the Bata and Implementing OPSlif % 

The collection of relevant data* required for this 
program, constitute a large portion of the present work. Collec- 
tion of data was started first m September 1974 and ultimately 
all the data could be collected from different places by the 
end of June 1975. Punching and correcting of the program 
started long before and ultimately m July 1975 the program 
was m perfect working condition. Only after that the work on 
the basic objective i.c. to use OPSIT for quantifying an Indian 
watershed, could be started. The work has been finally completed 
m October, 1975. 

3 • 5 -Factors Affecting the Efficiency of the Program : 

1. Due to heterogeneous properties m the large watersheds 
it is very difficult to correlate the model parameters with 
watershed character. Watershed should be small enough to be 
homogeneous to minimise this difficulty. 

2. The basin area should not be very large. 

3- One run of OPSEI program bakes a large amount of computer 

time and the computer expense is high. So for economic reasons 
it is desirable to use OPSGT for as few years as possible. 

4. Two fa.c tors cause the variation of the parameter value 

for Oil feront water years. These are (i) differences m the 
pattern of' deviation from year to year between gauge and wot.r*-’:. d 
rami' r 11s and (n) the differences m the pattern of the stre 
flows among the years. 



CHAPTER IV 


MODELING OR TH^ INDIAN BASIN AND THE 

INSULTS 

4.1 Description of the Basin . 

TRe basin in bordered by two canals on two sides (Fig, 8)*. 
The whole arts of the basin is more or less flat. The climate 
is temperate during hot summer months from May to September. 
Thereafter air becomes drier and cooler every day. Hoarfrost 
i ° common . u January and February. For few nights m the year 
the ground temperature falls below 0°C. Average mean monthly 
maximum and minimum tempera cures vaxy between 60°F to 110°f 
and 35°E to 80°F respectively. Northeast part of the area 
receives more than 50 inches (1270 mm) of rainfall annually 
and the southwest p. rt receives between 15 inches to 20 incues 
(381 mm to 508 mm). This area hss been formed of alluvial 
'o posits brought by the rivers. Some part is formed oC c Lay. 
laid clay or clay mixed with kankar at the top. 

| . 2 Da ta Availability ? 

I . 2 . 1 m easurable Watershed Data ; 

Che test data for the year 1956 for the basin Woot-e^ok 

'Che OPSb'T pro L i‘ ■ is 


it London* Xy. USA are available with us, 
Lirst run witii tins sot of data. 
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•For the basin m this study the following data were 
measured as accurately as possible. The description of all 
these are presented m Chapter 2. 

1. AREA? 

The maximum value of this is 650 sq. miles. The streamflow 
data were obtained from different places along the stream in 
different years. The areas under considerations are also 
different . 

2. PIMP ; 

The impurvious part of the area is signif icantly sma~ 

-L'his value has been taken as 0.0005. 

3. PWTR: 

There is no lakes or swamps inside the area. Hence thi ' 
value is taken as zero. 

4. VINTMR; 

Ho forest is observed inside the basin. The value of 
VINTER is taken from the Table 1 (Chapter II) as 0.10. 

5. GTkETP; 

As a significant number of phreatophy tes or swamps arc 
not known to exist inside the area, this value is taken as zero. 

6. SUBWF; 

It is assumed that no water leaves the basin through 
sub-surfacoflow not measured by the stream gauge after entering 
into the ground water storage from the surface. This value is 5 
therefore , taken as zero. 
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7. OFSS: 

This is obtained from the average of the measurements 
at randomly selected points on the basin. 4 measurements are 
taken from which tho average is calculated as described below? 



Average slope = 5.5 ft in 1 mile 0.001 ( = 0.1 percent) 

8. OF SI: 

OFSL can be estimated m different ways. In this study, 
it is represented as the reciprocal of twice the drainage density 
(obtained by dividing total lengths of all the streams within 
tho watershed m miles with the total area of the watershed in 
sq.miles). The calculations is shown below. 


Total arcs of the watershed = 650 sq. miles 

Total len th of all the streams 

inside the watershed = 456.0 miles 

456 0 

Hence, tho drainage density = 550^0 = 0.70 

01 3L - 2r ( d re mage density}" ~ 2 x 0.70 


0.715 miles 


3780 feet 
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9. OFMNs 

This value is taken from the table given below; 

Table 2s Manning's Roughness Coefficient for Overland 
Plow for Various Surface Types [9] 


Watershed Surface 

Manning's n 

Smooth Asphalt 

0.013 

Concrete (Trowel finish) 

0.013 

Rough Asphalt 

0.016 

r oncrcte (unfinished) 

0.017 

Snoobh barth 

0.018 

Firm Gravi. 1 

0.020 

Cemented Rubble Masonry 

0.025 

Pasture (Short grass) 

0.030 

Posture (High Grass) 

0.035 

Cultivated Area (Row Crops) 

0.035 

Cultivatea A rea (Field Crops) 

0.040 

Scattered Brush;, Heavy Weeds 

0.050 

Light Brush and Trees (Winter) 

0.050 

Light Brush and Trees (Summer) 

0.060 

Dense Brush (Winter) 

0.070 

Dense Brush (Summer) 

0.100 

Heavy Timber 

0.100 


OFM is taken for the basin as 0.035. 


30. ORMNISs 

This value is also taken from the above table. For tb 
present study ib is taken as 0.018. 
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11. CHOAP; 

The size of the channel at the gauging point may not he 
the size of the channel system as a whole. This can he measured 
from the hydraulic analysis of the profile and the cross section 
of the stream channel. Another method is to measure the gauge 
height of the hankful flow from the topographic map and then 
to read the GHCAP value directly from the stream gauge rating 
table. In this study CHCAP is indexed bo the outlet of the 
basin and is takon as more than the measured maximum discharge. 

12. DIY s 

There is no flow diversion into or out of the stream and 
hence this value is taken as zero for the watershed m the 
study. 

4.2.2 Meteorological Data : 

1. Evaporation Pats ; 

Pan evaporation data are available for one station 
near the basin. Evaporation data are recorded with mesh covered 
fixed point gauge, class 'A' pan evaporimeter. These values 
were m millimeters and are converted to inches hy multiplying 
with 0.03937. 

2 . Pan Coefficients s 

In this study these values are taken from the table given 
iii Pef. 10. Starting from the month of January upto December, 
there values are 1.00, 1.00, 0.90, 0.75, 0.70, 0.70, 0.75, 0.90, 
0.95, 0.90, 1.00 and 1.10. 



3. Daily Hourly Precipitation 1 Data ; 


These are recorded for one ram gauge station near the 
basin by a self recording rain gauge . These values are obtained 
again m mm and were converted to inches. 

4 . Daily Rainfall Data ; 

Daily rainfall data are also collected from an ordinary 
ram gauge station inside the basin. This ramgauge is called 
m the program as the secondary ram gauge. The values obtaned 
m mm were converted to inches. 

Evaporation, hourly rainfall, daily rainfall and discharge 
data are collected for a few years. All the data are not avai- 
lable tor the same place inside the basin. Two years data have 
been run but it is found that only m one particular year the 
rainfall data is m consistence with the discharge data. 

5 • Stream Flow Pat a . 

This data is collected from the drainage nallah of the 
basin. Discharge date are not available Tor the basin for the 
same place for different years. These are available m diffe- 
rent location m different years. 

The basin is epproxima;t ely rectangular m size and so 
the area covered by each stream gauge station is taken as 
linearly proportional to the distance of the stream gauge from 
the farthest end of the basin. 



4 . 3 Simulation A n alysis s 

4. 5.1 Initialization of Parameter Values ; 

OPSJ'J optimises 13 parameters altogether. These ere, 
LZC, Bi>!IR, 8UZC , E1LF, BUZC, SI1C, BIVF, BFRC, IFRG, CSEX, 
FbHX, JCi’RI, CiiCAP. 

Tho following ere the initial values of- first sin 
parameters inside the program. 

Debit; 3 • Initial Values of Paiameters (Ref. 3, p. 33 ) . 


Starting Values 


Parameter 

low 

Middle 

High 

LZC 

2.00 

12.00 

30.00 

'3 MIR 

0.20 

1.20 

4.00 

OUZO 

0.30 

1.30 

4.00 

iILF 

0.05 

0.25 

0.60 

3UZC 

0.20 

1.50 

5.00 

dliiC 

0.30 

0.90 

4.00 


For the consistent estimate of reasonable parameter 
"elaoj the medium starting values worked best because they 
' vo the best chance of being close to the final estimate rnd 
ess chance to produce out-of-range adjusted values whicn <iri 
, use the program to stop before the optimum point is re? ciioc 
Rej . 8, p. 34). 


The initial values of the other parameters are as follow.' 
BEVF = 0.0 , FSRX = CSRX = SEX = 0.98 

IFl’C = 0.1 , 8 FFlC = 0.9 


1 * . i pure 

£*!i I . library 
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4*5*2 Control Options ; 

1. COM OPT ; 

j-ho OPbiil program has the control options, expressed 
by COM 01 T (1), CONOPT (2) and CONOPT (3). 

uOifOiT (l) = 0, if daily evaporation data, are used. 

= 1, if evaporation is read by 10 day .period. 

= 2, if annual evaporation is read and the 

subroutine BVPDA.Y is used. 

CONG Pf (2) = 0, if stream routing is to he done with 

15 minutes time interval. 

- 1, if stream routing is to be done with 

hourly time interval. 

COMO PC (3) = 0, if no change in WSG- or SORT occur 

during the water year. 

= 1, if storage gauge was moved during the 

water jyear. 

in the test data. COMO IT (1), CONOPT (2) and CONOPT ("} 

'ere taken f*;; 2 V 0 and 1 respectively. Por the Indian data 

r. thin a Ludy tnoue are taken as 1,1, and 0. 

. /i.;kb , . , r".'H , and IflTRs 

these arc another set of control options which are given 

■. Low; 

HO a Minimum number of rough cycles to be made, 

NRi’K « Numb ox- of first trip to be run for a given year. 
li'L'K » Number of lost trip to be run for a given year. 
In ail th cases, the values of MNRC, NFTR and 11LTR 

r „ tr.iv n r.n 12, 1 and 5 respectively. 
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♦ 4 

Output 

-Parameters Using the Test Da-fen: 



x0S t 

deta [4] have been run through the OPSET program.. 

ompa 

rati vc 

study 

of these output parameters and that given in 

he report 

[8] is 

given in the following table. 



Table 4: 

Comparison of Obtained Parameter Values with 
those given in the Ref.. 8. 

31. 

ho. 

pc, 
X Ci 

r „ m tp Value given Value 

_„!l in_tho report from 

obtained 
the ^rogrom 

1 


LZC 

11.14 

11.76 

€ 


iJMIR 

4.23 

4.99 

3 


SUZC 

0.65 

0..65 

A 

*-* 


73 V LF 

0.15 

0.15 



<)Z-i 

1.04 

1.09 

8 


,j LTJ 

0.45 

0.45 

? 



0.00 

0.00 

re 

o 


OZKX 

0.935 

0.953 

9 


?-iRX 

0.935 

0.95S 

G 


UOTEI 

3- 

3 

1 


I/KC 

0.10 

0.10 

*3 

Pi-j 


UfhC 

0.895 

0.894 

V. 


9.. hn' 

100.0 

100.0 


In addition, the 1ZS value found in the report to he 

■)B axiu ti-rt o'bt-ineci from the program is 8.83. 

fiie synthesized md the recorded monthly flows and peal: 
• r y !•;/.( hn ;• re eon pa red and shown in Table 5 and Table 6 

; , t ,(* {; i vt; .],•/ on t in .• o feh*;.r pages . 


* 
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Table 5° Comparison of Recorded and Synthesized 
Monthly flows for the 'Test Data. 


Month 


Recorded flow 

Synthesized flow 

October 

1955 

15.5 

6.6 

November 

1955 

11.8 

8.2 

December 

1955 

21.5 

CO 

CO 

January 

1956 

65.9 

107.4 

February 

1956 

779.8 

719.9 

March 

1956 

469.4 

451.4 

April 

1956 

465.1 

456.9 

May 

1956 

160.5 

78.6 

June 

1956 

28.5 

50.0 

July 

1956 

202.5 

265.2 

August 

1956 

51.9 

49.5 

September 1956 

26.5 

117.5 

Yearly Total 

2294.5 sfd 

2511.5 sfd 


Table 6 ; Comparison of Recorded and Synthesized 
Hydrograph Peals for the Test Data. 


Date 


Recorded 

Synth 

esized 


Peakf low 
(cfs) 

Time 

Peakf low 
( cfs ) 

Time 

Feb. 

17 , 1955 

506.0 

11 PM 

458.4 

6.15 AM 

Mar. 

14, 1956 

264.0 

6 iM 

213.1 

6.00 AM 

Apr. 

6 , 1956 

153.0 

10 ill 

4 09.3 

9.15 All 

ipr. 

15, 1956 

227.0 

9 AM 

355.7 

8.30 ill 


From the results shown m the preceding tables , it is 
evident that the modol has boon implemented with reliable 
accuracy. Here the lowest value of SSQM = 
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Obtaining fairly accurate results with the test date 
the model was then applied to an Indian basin. 'The following 
are the results obtained for the basin. 


fable 1% Recorded and Synthesized Stream Flow 
For the Year 1971-72. 


Months 

Recorded 

Synthesized 

June 

20.0 

182.4 

July 

1306.0 

2849.6 

August 

12460.0 

9817.2 

September 

983-0 

585.9 

October 

0.0 

15.6 

November 

0.0 

0.5 

December 

0.0 

0.0 

J anuary 

0.0 

0.1 

February 

0.0 

2.2 

March 

0.0 

0.2 

April 

O 

* 

o 

0.3 

May 

0.0 

o 

• 

o 

Yearly Total 

14769.0 sfd 

13454.0 sfd 


Table 8; Comparison of Recorded and Synthesized 
Hydrogrcph Peaks for the Year 1971-72. 



Recorded 

Synthesized 

Date 

Peakflow 

(cfs) 

Time 

Peakflow 
(cf s ) 

Time 

July 30 
Augus t 5 
August 28 

110.0 

1950.0 

480.0 

12 AM 

12 AM 

12 AM 

121.0 

1898.8 

122.8 

0.15 AM 

1.00 iM 

0.15 AM 


Lowest SSQM = 2.674 
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Table 9; iiecorded and Synthesized Stresmflow for the 
year 1969-70. 


Months Recorded Plow Synthesized 


Plow 

June 1969 

0.0 

6.0 

July 1969 

5416.0 

554.5 

iugust 1959 

5850.0 

11171.9 

September 1969 

2191.0 

5464.9 

October 1969 

0.0 

557. 4 

■November 1969 

0.0 

41.0 

December 1969 

0.0 

11.8 

January 1970 

0.0 

8.8 

Pebruary 1970 

0.0 

19.7 

March 1970 

0.0 

12.0 

April 1970 

0.0 

1.5 

May 1970 

O 

* 

o 

1.4 

Yearly Total 

11457.0 sfd 

15710.9 sfd 


Table 10; Recorded and Synthesized Hydrograph Peaks 
for the Year 1969-70 


Recorded Synthesized 


Late 

Peak flow 
(cfs) 

Time 

Peak flow 
( cfs ) 

Time 

29th July 

'69 

425.0 

12 AM 

— 

* 

16th Aug. 

1 69 

540.0 

12 AM 

617.7 

9.15 AM 

15th Se t. 

? 69 

140.0 

12 AM 

125.2 

0.15 AM 


Lowest SSQM = 86.4999 
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Table 11 5 

Optimum Parameter Values Obtained for 
Indian Basin. 

Parameter 

Year 1969 

Year 1971 

1ZC 

12.0 

2.3582 

BMIR 

. 

11.6108 

SUZC 

1.3 

3.0 

3TLF 

0.092 

0.1793 

33UZC 

1.5 

0.2761 

SI AG 

4.0 

3.6 

BIVF 

o.o 

0.0 

BFRC 

0.9267 

0.8884 

I PPG 

0.10 

0.10 

CSEX 

0.9950 

0.90 

FSBX 

0.9950 

0.90 

IG l'RI 
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Table 13 i Synthesized Flow for Test Data at Wood Greek, London, Ky,USA, 1955-56. 
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Table 15s Synthesized Daily Flows (Indian Basin, 1971-72) 

(cfs ) 


DAT 

JUNE 

JULY 

AUGUST 

SEPT. 

OCT. 

NOV. 

1 

0.0 

70.8 

101.0 

70.2 

2.3 

0.1 

2 

0.0 

90.6 

841.9 

62.3 

2.0 

0.1 

3 

0.0 

185.8 

1277.3 

59.5 

1.8 

0.1 

4 

0.0 

194.6 

1115 . 2 

55.9 

1.6 

0.1 

5 

0.0 

208.3 

724.9 

49.6 

1.4 

0.0 

6 

0.0 

207.0 

591.2 

44.1 

1.3 

0.0 

7 

0.0 

192.5 

517.5 

39.2 

1.1 

0.0 

8 

0.0 

175.5 

458.5 

34 .8 

1.0 

C.O 

9 

0.0 

158.1 

407.0 

30.9 

0.9 

0.0 

10 

0.0 

141.4 

361.5 

27.5 

0.3 

0.0 

11 

0.0 

125.9 

321.0 

24 .4 

0.7 

0 « u 

12 

0.0 

111.9 

285.1 

21.7 

0.6 

0 . 0 

13 

0.0 

99.4 

253.4 

19.2 

0 . 6 

o . o 

14 

0.0 

88.3 

224.9 

17.1 

0. 6 

0.0 

15 

0.0 

78.4 

199.7 

15.2 

0.5 

C.O 

16 

0.0 

69.6 

177.4 

13.5 

0.4 

0.0 

17 

0.0 

51.8 

190.9 

12.0 

0.4 

0.0 

18 

0.0 

71.6 

174.1 

10.6 

0 . 1 - 

0.0 

19 

0.0 

79.5 

154.6 

9.4 

0.3 

0.0 

20 

0 . 0 

70.6 

137.3 

8.4 

0.3 

0.0 

21 

0.0 

62.6 

121.9 

7.4 

0.2 

0.0 

22 

0.0 

55.6 

108.3 

6.6 

0.2 

0.0 

23 

0.0 

49.4 

96.2 

5.9 

0.2 

0.0 

24 

0.1 

70.7 

85.4 

5.2 

0.2 

0.0 

25 

2.3 

69.0 

75.8 

4.6 

0.2 

0.0 

26 

20.8 

68.8 

79.9 

4.1 

0.1 

C.O 

27 

30.7 

65.9 

113.9 

3.7 

0.1 

0.0 

28 

29.2 

58.5 

112.8 

3.2 

0.1 

0.0 

29 

28.4 

94.8 

100.2 

2.9 

0.1 

0 . 0 

30 

27.8 

128.1 

89.0 

2.6 

0.1 

0.0 

31 

- 

113.8 

79.0 


0.1 

~* 


Flows m the other months are all zero. 
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Table 16s Recorded Daily Plows (Indian Basin, 1969-70). 
l£fsl 


DAY 

JURE 

JULY 

AUGUST 

SEPT. 

OCT. 

• 

0 

r— i 

1 

0.0 

0.0 

305.0 

102.0 

0.0 

0.0 

2 

0.0 

0.0 

305.0 

102.0 

0.0 

0.0 

3 

0.0 

0.0 

283.0 

100.0 

0.0 

c.o 

4 

0.0 

0.0 

262 . C 

100.0 

0.0 

0.0 

5 

0 . 0 

0.0 

242 . C 

100.0 

0.0 

0.0 

6 

0.0 

0.0 

242.0 

90.0 

0.0 

0.0 

7 

0.0 

0.0 

182 . 3 

85.0 

0.0 

0.0 

8 

0.0 

0.0 

150 . 

65.0 

0.0 

0.0 

9 

0.0 

45.0 

12 c. 0 

65.0 

0.0 

0.0 

10 

0.0 

40.0 

118.0 

110.0 

0.0 

0.0 

11 

0.0 

45.0 

115 . 0 

115.0 

0.0 

0.0 

12 

0.0 

50.0 

115. c 

120.0 

0.0 

0.0 

13 

0.0 

55.0 

113 .i 

134.0 

0.0 

0.0 

11 

0.0 

58.0 

134 . 0 

125.0 

0.0 

0.0 

15 

O.0 

60.0 

140.0 

110.0 

0.0 

0.0 

16 

0.0 

62.0 

327.0 

100.0 

0.0 

0.0 

17 

0.0 

66.0 

305.0 

10O.0 

0.0 

0.0 

18 

0.0 

66.0 

305.0 

90.0 

0.0 

0.0 

19 

0.0 

70.0 

262.0 

92.0 

0.0 

0.0 

20 

0.0 

66.0 

184.0 

86.0 

0.0 

0.0 

21 

0.0 

60.0 

262.0 

60.0 

0.0 

0.0 

22 

0.0 

50.0 

134.0 

45.0 

0.0 

0.0 

23 

0.0 

66.0 

150.0 

40.0 

0.0 

0.0 

24 

0.0 

66.0 

167.0 

30.0 

0.0 

0.0 

25 

0.0 

150.0 

160.0 

25.0 

0.0 

0.0 

26 

3.0 

327.0 

131.0 

0.0 

0.0 

0.0 

27 

0.0 

350.0 

134.0 

0.0 

0.0 

0.0 

28 

0.0 

4 24 . 0 

120.0 

0.0 

0.0 

0.0 

29 

0.0 

450.0 

115.0 

0.0 

0.0 

0.0 

30 

0.0 

100.0 

110.0 

0.0 

0.0 

0.0 

31 

— 

380.0 

102.0 


0.0 

_ 


Plows m the other months are all zero. 



Table 

17s Synthesized 

Daily Flows 

Lq£&1— 

(Indian Be 

'Silly 

1969-70) 

DAY 

JUNE 

JULY 

AUGUST 

SEPT. 

OCT. 

NOV. 

1 

0.0 

1.2 

39.9 

270.0 

30.3 

0.0 

2 

0.0 

1.2 

77.6 

250.2 

28.1 

0.0 

3 

0.0 

1.1 

96.0 

231.8 

26.0 

0.0 

4 

0.0 

1.0 

215.9 

214.8 

24.1 

0.2 

5 

0.0 

0.9 

258.4 

199.0 

22.4 

0.0 

6 

0.0 

0.9 

241.0 

184.4 

20.7 

0.0 

7 

0.0 

0.8 

223.3 

171.0 

19.2 

0.0 

8 

0.0 

2.7 

218.2 

158.7 

17.8 

0.0 

9 

0.0 

3.2 

295.4 

147.0 

16.5 

0.0 

10 

0.0 

3.0 

330.5 

136.2 

15.3 

0.0 

11 

0.0 

3.1 

307.9 

132.3 

14.1 

0.0 

12 

0.0 

3.3 

286.0 

127.8 

13-1 

0.0 

13 

0.0 

5.4 

265.2 

118.4 

12.1 

0.0 

14 

0.0 

6.7 

246.0 

109.7 

11.3 

0.0 

15 

0.0 

6.6 

262.0 

101.7 

10.4 

0.0 

16 

0.0 

7.4 

577.6 

94.0 

9.7 

0.0 

17 

0.0 

10.2 

594.7 

87.3 

8.9 

0.0 

18 

0.0 

9.2 

613.8 

81.5 

8.3 

0.0 

19 

0.0 

8.8 

627.2 

75.7 

7.7 

0.0 

20 

0.0 

8.2 

582.7 

70.2 

7.1 

0.0 

21 

0.0 

20.0 

539.9 

65 .0 

6 . 6 

0.0 

22 

0.0 

26.9 

548.5 

60.3 

6.1 

0.0 

23 

0.0 

40.4 

534.4 

55.8 

5.6 

o.c 

24 

0.0 

72.9 

495*5 

51.7 

5.3 

0.0 

25 

0.0 

67.7 

460.6 

47.9 

1.9 

0.0 

26 

0.1 

62.7 

426.7 

44.4 

4.5 

0.0 

27 

0.2 

58.1 

395.4 

41.0 

4.2 

0.0 

28 

1.6 

49.9 

366.4 

38.1 

3.9 

0.0 

29 

1.5 

53.8 

339.5 

35.3 

3.6 

0.0 

30 

1.3 

46.4 

314.5 

32.7 

0.0 

0.0 

31 

- 

43.0 

291.4 

- 

0.0 

— 


Flows m the other months are very negligible 



CHAPTER V 


DISCUS SI OHS AND CONCLUSIONS 

5 . 1 Discussion of the Results ; 

5.1.1 SSQH s 

SSQM is tiie sum of squares of deviations for 11 months 
and it is used to indicate which set of volume parameter values 
gives the best synthesis of flow volumes. The first month of 
the water year is excluded from this least square term because 
the flows m tho first month are too dependent on unlrnown 
initial conditions. Tho final value of SSQM gives a reliable 
indication as to how best the parameters a.re synthesized witn 
respect to the input data. 

5.1.2 Test Data ; 

After going through the results from the tost dat'yit is 
observed that the optimized 13 output parameters are closely 
co'iipt rablc with the reported values [8], The recorded daily 
stream-flow and the synthetic flow, as shown m the Tablesi2,13 
a.re also comparable. In addition, from Table 6, it is clear 
that peak hydrogrrph values also matched well. Hence the 
OPSEf program is well l m lamented , m tho computer system 
f variable with us. Monthly flows are compared m Pig. 9. 

5.1.3 Indian Data ; 

After analysing the output for 2 years for Indian date 
th^ following observations can be made. 
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Year 1969-70 ; 

From the table 16,3 7P-t is observed that the daily recorded 
stream flow did not match with the daily synthesized stream flow. 
Similarly the monthly flows and the yearly totals of recorded 
and synthesized stream flows have wide gap. Peak hydrograph 
values between the observed and the synthetic cannot bo compared. 
Hence the parameter values obtained on the basis of this year 
data arc not correct. SSQM value m this year is high. The 
daily recorded and synthesized flows are shown m Fig. 11. 

Y ear 1971-72 : 

From the figure 10 and the Tables 14 and 15 of tho 
daily stream flow hydrograph and the monthly totals respec- 
tively for the observed and the synthesized values , it is seen 


thee these two matched well. Again the recorded hydrograph 
peck values can be compared with the synthesized hydrograpti 


peak values. 

The SSQM 

value is 

very low m 

this year 


uho parameter 

values 

for the ye 

ar 1971-72 

arc taken 

as t ic 

optimum parameters for 

the basin 

These values are; 


JjZC = 

2.338 

, Bi-TIE 

= 11.610 , 

SUZC = 

3.00, 

FTLF = 

0.1793 

, BUZC 

= 0.2761 , 

SI AC = 

3.60, 

■ilVF - 

o 

9 

o 

, BFRO 

= 0.8884 , 

IFItC = 

CN 

o 

H 

« 

O 

CSRX - 

0.90 

, FSiU 

= 0,90 

NCTRI = 

25 and 

CHCAI : 

= 2000. 
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f i'-’l-l RFCOKUFD and synthesized daily flows of 

INDIAN BA Dili WATEP YEAR 1971-72 
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It is clear that the results obtained for the Indian 
basin a.re not really accurate. Also, the results from the 
different years are not equal to each other. This is basically 
due to the following factors. 

(i) All the rain gauge stations whose recorded hourly and 
daily rainfalls, used in the present work, are not situated 
within the basin. So the rainfall pattern fed to the program 
is different from the actual pattern. 

(ii) Streamflow records are not reliable. Some values are 
missing and some are not observed. 

(iii) Topographic sheet for the area is not available. 

(iv) Watershed area is large. 

(v) Unlike the test data for rainfall, the rainfall in the 
basin under study occur mainly for 3 months in a year. 

(vi) The peak flow discharge data are not ava.ila.ble. 

(vii) The streamflow gauging points are different in 
different years. 

However, the parameter values obtained from the data 
for the year 1971 are the bettor of the two. They can be taken 
rs very good approximation afor the parameters of the water- 
shed. This can be inferred from the following reasons s 

(i) In the year 1971 the rainfall and streamflow patterns 
are comparable. 

(ii) SSQM value is very low in this year. 
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(ill) With days passing, the data availability also is 
more. 'the data which are not available for the preceding 
years are available for the year 1971. Hence, the results 
also can be considered better. 

5 . 2 Conclusions end Scope of further Studies . 

through the OPSHT program very good prediction of annual 
streamflow and fairly good prediction of monthly streamflov ar' j 
obtained. Flood hydrographs are simulated less accurately. 

Natural phenomena being simulated are frequently too 
complex to permit an intricate analysis and to simplify these 
further research is necessary. 

Measurable watershed parameters should be measured 
very accurately and the proper authorities should be well 
informed m advance to keep the necessary data, for further 
studies. Peak flows are very much needed for estimating the 
channel routing parameters. 

OPSEi 1 program should be applied m different water -.hod 0 
varying m sizes, locations and climetological conditions. 

Effect m the changes m the watershed such as urbani- 
zation and agricultural land use should also be studied through 
this program. 

For larger watersheds , it is better to subdivide the 
total area into parts, each with its own set- of parameters, 
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and combine these results by stream routing. 

further oc bailed studies should be made to find out 
tho ground water re-charge through this model. 
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APPENDIX 


COMPARISON OF ANNUAL RUNOFF CALCULATED BY OPSET METHOD WITH 
THAT BY AN APPROXIMATE METHOD 

TEST DATAs Year 1956, Wood Creek, Ky, USA. 

At> proximate Method ; 

Runoff = Rainfall - Evapotranspiration 
In this year of test data 

Rainfall = 53.86 inches. 

Evapotranspiration (Potential) = 35*785 inches. 
Annual Runoff = (53.86 - 35.785) = 18.075 inches. 

Totalarea of the watershed = 3.89 squa^fe, miles " 


Annual Run off * 2 x 18 *° 7 5 = 1920 sfd 

Net evapotranspiration calculated hy OPSET program i 
29-094 inches. 


Hence Annual Runoff = (53.86 - 29.094) = 24.766 mcr ' . 

= x 24 766 = 2582 ° "0 

12 x 86400 X dq.fbb .l 

Annual Runoff calculated hy OPSET = 2311.5 sf d . 

It seems that the net evapotranspiration value gives toe 

tetter result m Approximate method. 

For Indian data the total annual evapotranspiration value; 

always exceeded the net precipitation m the year and hence the 

above method could not be used. Khosla’ s method could not be 


applied due to non-availability of the value of the constant for 
the catchment under consideration. 
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LISTING OF THE COMPUTER PROGRAMME 
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SHMOCC lu 
SHMOO02O 

SHM0003t 
SHM0004C 
SHM0005U 
SHMQ006Q 
SHMGOC 1C 
SHM00080 
SHM00090 
SHMOOIOO 
SHMOOl 10 
SHMO0120 
SHMOOl 30 
SHM00140 
SHMOOl 5G 
SHMOOl 60 

COMPUTER PROGRAM ON HYDRO LOO IQ SIMULATION SHM00170 

0PS6T : SHM00180-- 

A SELF-CALIBRATING VERSION OF THIRST AN FORD WATERSHEDMQDEL SHM0G190 

BASIC LOGIC OF INNER LOOP BASED ON STANFORD WATERSHED MODELS III AND XSHM00200' 
VERSIONOF NOVEMBER 12,1970 SHM00210 i 

J$ca$c## #### * #2$: 3^ $$$ $ HMQ022-0 : ' 

SHH0023O 
SHMQ0240 

MENTED FOR THE SYSTEMS IBM 7044 AND IBM 370/155 SHMOG250 
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SHMO027© 

****************************************************************** ***$HMQ0280 
IN THIS STUDY WATER YEAR STARTS FROM JUNE 1 
DAY * 274 IN THE PROGRAMME IS TAKEN FOR JUNE 1 
MONTHS IN THE TABULATED RECORDED AND SYNTHESIZED FLOWS ARE CHANGED 
IN THE NEW FORMAT ALSO DATES SHOULD BE READ AS FOLLOWS 


q**jM**************** JUNE 1 AS JUNE 1 
C* ******** *********** JULY 1 AS JULY 2 
q* ******** *********** AUG 1 AS AUG 2 
C* ************* ****** SEPT 1 AS SEPTl 
C* ************ ******* OCT 1 AS OCT 2 
C* ******** ************ NOV 1 AS OCT 30 
£******************** DEC 1 AS NOV 30 
C* ******************* jan X AS DEC 30 
C* ******************* feb 1 
C* ******************* MAR 1 
C# ******** ************ APR 
q* ****** ************** MAY 


AS JAN 30 
AS MAR 1 
1 AS APR l 
1 AS MAY 2 


2§S ^ 3§C J$i 3§g 3$C 3i|c 2 $ 5jt !^C^£ 3^ ^ 3|<I 3$S 3§t 9§t 3^ S^t 1 

2 §C j|C3§£.i^E3|r%5^^5 5|S3§C 3$: 3§C ^ 

nr nr ntf* ■***• *nr ‘■r* hf* ■hp ■*©* •**r *w' nir'nir; i , nrf' l 3 P w ' ;v-, : 
2 $C ^C3$C3^£^;3^^C^C3{c^C3^:^C3jt3|{'3^(^CaJC35 c ^r3^:3{lC3^3^^C^(C3§C3^^ 3^C 

♦ 3|£ ^ a^c 3|C 4e $ £ 4c £ $ £ $ £ $ $ # £ sfc ^ 

# ajc 5fc 3j!f 3jc ^ 3$C 5^ $ $ $ $ $ $ # $ * ^ * # # ♦ # 

3$C 3^ 3^5 3^^? jj? 3^£ jjt 3^ 3^f J^C 5^ 3^C 3^C3j< 3^C ^ 
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*£C al 3BE' a f K%£ 3^ •gaf ,!PT 3£X 30C- 33C jP*. *r ’»f T3& ■jfflE Im, ••.•■■aBRv, i lBtfi :;g. •■•• 

nr np mr nr ,> np < npmr nP “ vr np “ 'nr ,- T 

3jt aj( J^n 3^ 3^; jjj; ^ 3 ^; ^ ^ 5 ^ jj|{ ^ j|j Jjjj{ ^|g J$( ^ 3^ ^ J^C jgg 3 ^ 

^g^s^f,s|j;3|f3jC3^{3^C3{(: ^ ^ jg ^ ^ * ^f ^ gjj ff ^|f ^ ^ ^ jjt 

***************************** 


• «©#©•©«© ©»©•#■ 


*©©#©•««© ©•'•••©•••••©©•©©#©©©©©####©©#©## i ®(#© : ' 


SHM00290 

SHM./03OO 

SHM00310 

SHM/Q320 

SHM00330 

SHM00340 

mmmm 

SHM00360 

SHM00370 

SHM00380 

SHMC039r 

SHM00400 

SHM00410 

SHM0042C 

SHM' * : 

SHM0044U 

SHM0v450 

SHM0U46C 

SHM0C470 

SHM004B0 

SHM0049C.I 

SHM00500 

SHM0051U 

SH MOO 5 20 

SHMOO 530 


f 


DIMENSION BTRK99), C0N0RK5 ),CTRI ( 99) , DR6PH (366 ), DRHP ( 366,24 ) 

1 DR SGPC 3 66 ) » CPET < 366 ) , DRSFB66), DSSF(366), EHGWSU2), 

2 EMIFSC12), EMLZ$(12),EMSIAM{12) , EMUZC < 12 ) , EMUZ $ ( 12 ) » 

3 FPCMU2), HBF{ 5) , IDYB ( 5 ) » IDYE(5) , I HRB ( 5 ) » IHRE ( 5 ) » KPSH{5) 

4 LSHA ( £ ) , MEDC Y{ 1 2) ,M£OWY( 12 ) »RHPD( 5 ) ,RHPF£5 ) ,RHPH(5) 

DIMENSION RS6BF(20) f RSBD(20),RSBIF(20) t SBFRSt3*20) 

6 » THSF (24 ) » TITLEC20), TMBF( 123 ,TMIF( 12) , TMN£T(12), TMOF(12) , 

7 TMPETCX2), TMPRECC 12) , TMRTFU2), TMS6C12) , TMSTF< 12), 

8TMSTFI {12),UHFA(99),XMPFT(12),DUM(368),SIFRS{3,20),SSR<5,170) 

LOGICAL LBMIR* LBUZC, LETLF , LLZC, LNPR, LRC, LSHA, LSHP 
INTEGER CN, CCNOPT, DATE, DAY, DPY , EHSGD, HOUR,HRF,HRL, PDAY, 

1 PRD,RHPD»RHPH» RSBD, SGMD, SGRT, SGRT2, TRI P, YEAR , YR1 , 

2 YR2 

REAL I FPRC, IFRC, 1FRL,IF5,LZC ,LZRX,LZS,LZSR, MNRD, NHPT 


SHMD0540 
SHM( - ' 55 ’ 
SHMOO? 60 
SHMOO £7f 

SHMt C £8C 

3 HMD* §94 
SHM00600 
SHM0061C, 
SHM00620 
SHMOO 630 
SHM0064C 
SHM0G65G 
SHM0066C 
SHK0067C 
SHM006B0 


DATA MfcDCY/ (2, 31,59,90,120,151,181,212,243, 

DATA MEDWY/3C4, 334, 365,31,59,90, 120,151,181,212,243 
SPECIFY NUMBER OF STATION- YEARS INCLUDED IN COMPUTER RUN 

NSYC * 0 
READ 500* NS YT 
500 FORMAT *51 5:1 
100 NSYC NSYC + 1 
READ TITLE TO COMPUTER RUN 
READ 15,1) TITLE 


273, 

,273 


304,334/ 

J 


1 FORMAT! 20 A4) 

C READ CONTROL OPTIONS 


READ 101 , ( CCNOPT (I), 1*1, 3) 

101 FORMAT (31 2) 

READ 101 , MNRC »NFTR,NLTR 
C READ BASIC TIME-AREA HISTOGRAM 
DO ID 2 KI A * 1,99 
BTRKKIA) =0.0 

102 UHFA(KIA) = 0*0 
READ 500, NBTR I 

READ l r 3, ( BTR I(J ) , J=1»MBTRI ) 

1C 3 FORMAT ( SF 10* 5 ) 

C S C T INITIAL CONDITIONS 
I FT * 1 
LRC=.TRUE. 

LLZC = . FALSE. 

LBUZC =. FALSE. 

LBMIR FALSE. 

LETLF = .FALSE. 

LNPR = .FALSE* 

1 F ( CONOPT ( 2 ) ,EQ . 0 .AND. NBTRI.LE. 6) LNPR= .TRUE. 

KRC= 1 
KBRC * 0 
KFFC « 0 

SSSQM= 950.0 ■ ■ . : 

SGRT * 0 

C READ FIXED PARAMETERS 
1C 4 READ 103 , RMPF ,CHCAP 

READ 1»»3 , RGPM E,AREA , FI MR, FWT R 

READ 1* J, VINT PR, SUBHF,GWETF, OFSS,OFMN,OFMNIS,OFSL,DI 


SHM00690 
SHM00700 
SH MO 07 10 

SHM00T2G 
SHMOO? 30 
SHMOO 740 
SHMOO? 50 
SHMOO? 60 
SHM0O770 
SHM0078© 
SHMOO? 9© 
SHM.'CiS^ 
SHM* SH 
SHM00820 
SHM0083© 
SHMu*840 
SHMOftm 
$MM0086< 
SHMOO 870 
SHMOOtBO 


SHMOO 900 
SHM00910 
SHM00920 
SHM00930 
SHM00940 
SHM00950 
SHM00960 
5HM0O97© 
SHMO'0980 
SHMOO 990 
SHM010Q0 
SHM01O1G 
SHM® 1020 
SHM01O3G 
SHM01O4G 
SHM01050 
SH Mil 10 60 
SH MO 10 70 


C CALCULATE CONSTANTS SET BY FIXED PARAMETERS 
FPER = 1.0 -FlMP-FWTR 
IFCFPER .GT. C.01) GOTO 105 
TPLR * 100.0 
FPER « 0.01 
GO TO mb 

105 TPLR - ( 1.0~FWTR)/FPER 

106 VWIN * 26.8888*AREA 

wcfs - 2t»o*Via n 

RHFMC * 0*025 /WCFS 
SSRT « SGRTC'OFSS) 

OFRF * it>20.O*SSRT/{0FMN *OFSL) 
OFRFJS*1O20.O*$SRT/(OFMNIS*OFSL) 

EQDF* G.0G982*<CQFMN*OFSL/SSRT)**O.6) 
EQDFIS * Q.00982*( (0FMNIS*0FSL/SSRT)**0.6) 
RGPM “RGPMB 
C READ WATER YEAR 

READ 10||, YR1, YR2 ' 

DRY ■' 

IF(MdD(YR2,4I .EQ. 0) DPY = 366 

C READ 1 EVAPORATION OATA 

CTiOPTf 1) *NE* 1) GO TO 111 

1%I»CDPET(KRD)»KRD*2T4,360»10) 

9C KRD=274,360,10 

(KRD } « 0.03937*DPET(KRD) 

103 * C DPET(KRD) » KRD=1» 273, 10 ) 

; PRINT 103, (DPET1KRD), KRD= 1,273,10) 

DO 991 KRD * 1,273,10 
991 DPET(KRD) » 0.03937*DF6T(KRD ) 

DO 110 I DAY2- 1,9 
DO 109 1 DAY1= 274,360,10 
DAY® 1DAY1 +ICAY2 
109 DPET { DAY ) = DPET ( IDAYl ) 

DO 111* I DAY 1= 1,273, 10 
DAY=I DAY1 + IDAY2 
I F ( DAYcGT *> 273) GO TO 110 
DPET (DAY) =DPET (IDAYl ) 
lift CONTINUE 

DPET ( 366 J = DPET 1 59 ) 

DPET (365 ) = DPET (363) 

DPET (364) = DPET (363) 

GO TO 113 

111 I F ( CONOR T ( 1) « FQ« 2) GO TO 116 

DAY =274 

112 READ 103, DPET (DAY) 

IF ( DAY .EQ. 213) GO TO 113 
CALL DAYNXT( DAY, DPY) 

GO TO 112 

113 READ 103, (EPCM MONTH), MONTH- 1,12) 

PRINT 13, <F PCM (MONTH) , M0NTH=1,12) 

EPAET*0.0 

DO 115 DAY =1 ,0P Y 
115 EPAFT =EPAO +DPET ( DAY ) 

1 F ( K PC M ( 6 ) .NIE® 1.0) EPAET =0.7*EPAET 


SHM0108C 
SHMl,a09f 

SHM011G0 
SHMUlll 
SH MO 11 20 
SH MO 1130 
SHM01140 
SHM01150 
SHMC 116‘ 
SHM01170 
SHMOllSO 
SHM011 90 , 

SH MO 1200 j 
SHM0121O 
SHMC1220 
SH MO 1230 | 
SHMC 1246 
SHM01250 
SHM01260 
SHM01270 
SHM01280 
SHM0T290 
SHM01300 
SHM01310 
SHMO1320 
SHM0133C 
SHHttl34C 
SHMCU35C 
SHM01360 
SHMM37C 
SHMi)l'38f 

smimm 

SHM014O0 

SHH0141G 

SHM0142O 

SHM01431? 

SHMG144C 

SHM01450 

SHM01461 

SHM0147< 

SMM0148< 

SHM0149* 

SH MO 1501 

SHM01510 

SHM01521 

SHM0153C 

SHM0154P 

SHM01551- 

SHM0156L 

SHH0157G 

SHW0158I. 

SHM01590 

SHXK: 16* •. 

SHM0161 


GO TO 117 

PRINT 103, PPAET »MNRD 

116 R PAD 103, r.PAET,HNRD 

HMAFT =EPAFT* (365*0 + MNRD)/4O4.0 
CALL E VPDAVt DPET »EMAET ) 

C RFAD DAILY FLOW DATA 

117 OR S F ( 3 66 I * 0*0 
11R DAY *274 ' 

READ 1000, ( DLM( J) , J*l, DPY) 

PRINT lOOOf(Dt'H( J) ( J=l, DPY J 
I SEQ*1 

119 DRSF (DAY ) *DUM (I$EQ ) 

ISEG-ISEQ+1 

1000 FORMAT (13 F6. 1 ) 

CALL DAYNXT(DAY,DPY) 

I F ( DAY .NE* (396-DPY) ) GO TO 119 
IF( DIV .EQ. 0,0) GO TO 122 
'DO 121 DAY * 1,DPY f* 

IFCDRSF(DAY) *GT • DIV) GO TO 120 
DRSF(OAY) * O.C 
GO TO 121 or 

120 DRSF (DAY ) * DRSF(DAY) -DIV 

121 CONTINUE 

«|’‘f 2J WRITE! 6, 2 ) (TITLE(KTA), KTA =1,20) 

2 FORMAT (1H1»25X»2QA4) 
t ‘WRlfg., DAILY FLOWS 

CALL DAYS UM( DRSF ,MEDCY» DPY,RATFV,TMRTF) 

' WRITE (6, 3) 

3 FORMAT (1H0»42X, * RECORDED FLOWS * ) 

CALL DAYOUT (DRSF, MEDWY, DPY) 

WRITE(6*4) (TPRTF(KWD), KWD *1,12) ,RATFV 
A FORMAT (6X , *T0 TAL * , 2X, 12F8. 1 , 2X, F10. 1, 2X,3HSFD) 

C READ STORM HYDROGRAPH DATA 
READ 530, NRHP ,NHPT 
PRINT 510, NRHP, NHPT 
510 FORMAT ( 1 5 , F10 *5) 

I F ( NRHP o GO. 0 ) GO TO 124 
DO 123 KRD*1, NRHP 

READ 90, RHPD (KRD) ,RHPH(KRD) ,RHPF( KRDi 
90 FORMAT (2I5»F1C*5) 

123 WRIT€( 6, 5 ) KR E, NHPT, RHPD( KRD ) ,RHPH ( KRD) , RHPFf KRD ) 

5 FORMAT (//5X, ’ RECORDED HYDROGRAPH • , 13/ 10X , • HYDROGRAPH INTERVAL * 
1F5. 2, IX, 5HH0URS/10X, ‘CALENDAR DAY OF PEAK = * , 15, 5X, * HOUR OF DAY 
2, 14, 5X, ‘PEAK FLOW = • , F8« 1 , IX , 3HCFS ) 

C INITIALIZE PRECI PIT ION DATA ARRAYS • ' 

124 00 125 DAY * 1,366 

DRGPM( DAY ) =RG FMB 
DSSF(DAY) *0,0 
DRSGP(DAY) =0 «0 
DO 125 HOUR = 1,24 
125 DRHP (DAY, HOUR ) =0.0 
C READ AUXILIARY RAIN GAGE DAILY TOTALS 
READ 500, NSGR C 
I FUlSGRD • EQ* 0) GO TO 127 


90* SGRT 2, SGMD, W$G2 


98 CARO 


Rt;AD 3 ID, SGRT ,WSG 
1F(CONOPTO).CQ.1)READ 
DO 126 KRD * 1,NSGRD 
READ 133 f DRSGP (ISGRD) 

Rfc'AD 5 GO, ISGRC 

126 DRSGP( ISGRD) = f!U0393?*DRSGP {I SGRD ) 

C READ RECORDING RAIN GAGE HOURLY TOTALS 

127 READ 500, IWBG ,YEAR»MQNTH,DATEtCN 
C PUNCH NO NUMBERS AFTER CN ON YEAR .ECU 

I F{ Y EAR .GE. 98) GO TO 130 

HRF *12#{ CN» 1 ) +1 

HRL »12*(CN~1 ) + 12 

DAY « MEDCY{ MCNTH) + DATE + 122 

IFCDAY .GT. DFY) DAY* DAY-DP Y 

DO 2000 I JK * 1,12 

IF (DAY «>GT* PEDCYUJK)) GO TO 2000 

MONTH * UK - 1 

GO TO 2059 / 

CONTINUE 
CONTINUE * 

READ 1003* (ORHP (DAY, HOUR) ,HOUR=HRF, HRL) 

DO 1291 HOUR® HRF, HRL 

DRHPC DAY, HOUR ) = 0* 03937*0 RHP ( DAY, HOUR ) 

FORMAT (12F5* 2 ) 

IFCDPY .NE. 366 .OR. MONTH . NE. 2 .OR. DATE 
001:129 HOUR=HRF» HRL 
ORHP (366, HOUR ) = DRHP (60 » HOUR ) 

DRHP (60, HOUR ) *0.0 
GO TO 127 

C CALCULATE PRECIPJTAIQN 


2000 

2050 


1291 

1003 


129 


,NE. 29) GO TO 12? 


13'' 


131 


132 


0) GO 


WEGHTING 
TO 137 


FACTORS 


I F ( NSGRD • EQ 
PDAV = 274 
ROPT =0.0 
DAY *274 
EHSGD =SGRT 

I F ( SGRT. EQ. 0) EHSGD =24 
EHSGDF =EH$GD 

CONTINUE 

DO 136 HOUR = 1,24 

RDPT =RDPT +DRHP (DAY, HOUR) 

I F ( HOUR .NF. EHSGD) GO TO 


I F ( RDPT • LFo 
I F{ SGRT .EG. 
DRGPM (PDAY) 
l F ( CONOPT ( 1 ) 


136 

C.0) GO TO 133 
C) PDAY =DAY 

= (DRSGP (DAY ) *WSG +RDPT* ( 1.0-WSG ) ) /RDPT 
•NE. 0) DPET(PDAY) =0.5*DP£T(PDAY) 


,NF. .0 ) PDAY *DAY 


133 

134 

135 

136 


I F ( SGRT 
RDPT * 0.0 
GO TO 136 
IF(DRSGP( DAY) 

DO 134 KHOUR 
DRHP ( DAY, KHOUR) 

I F ( SGRT .NE. C) 
CONTINUE 

CALL DAYNXT(DAY,DPY) 


.LE. 0.0) GO TO 135 
=1, EHSGD 

»(WSG*DRSGP( DAY ) ) /EHSGDF 
PDAY=DAY 


S HMD 2 I 60 

SHHC 21?C 

SHM0218P 

SH MO 21 90 

SHMU2200 

SHM02210 

SHM02220 

SHMG2230 

SHM0224G 

SHM02250 

S HMD 2 2 60 

SH MO 2 2 70 

SHM02280 

SHM02290 

SH MO 2 3 00 

SH MO 23 10 

SH MO 2320 

SHM02330 

SHM0234O 

SHM02350 

SHM02360 

SHM02370 

SH MO 23 80 

SHM02390 

SHM02400 

SHMO2410 

SH MCI 24 20 

S HMD 2430 

SHM02440 

SHM02450 

SHMU2460 

SHMC24T0 

SHM *2480 

SHM02490 

SHMC250O 

SHMU2510 

SHM02320 

SHMO2530 

SH M0 2540 

SHMO2550 

SHM0256C 

SHM02570 

SHWJ2580 

SH MO 2 590 

SHMG260O 

SHMO2610 

SHM02620 

SH MO 26 30 

SHM02640 

SHM0265O 

SHttO 266ft 

SHM02670 

SHMG268D 

SHM0269CJ 


c 


c 


IF {DAY .HQ. 274) GO TO 137 
I F { CQNOPT C3 ) .EQ.O) GO TO 132 
I F { DAY e MEo SGMD) GO TO 132 
WSG =WSG2 
SORT * SGRT2 
GO TO 131 

ADJUST RAINFALL ANOMALIES 

137 MXTRH »2*NBTRI 

IF { CON OPT { 2) .EQ. 0) MXTRH = (2*NBTRI -l)/4 +1 

NATRH wMXTRH/ 2 

IF(NFTR .GP. 2) GO TO 138 

I F (NATRH .LT. 12) CALL PRECHKCDRGPM, DRHP ,DRSF , VWIN, SGRT , NATRH > 
SET INITIAL VALUES OF VARIABLE PARAMETERS TO BE OPTIMIZED 
LZC =12*0 
BMIR » 1*2 
SUZC * 1.3 
ETLF « 0.25 
SUZC *1.50 
SIAC =0.90 ■ 

BIVF =0.90 

138 I F ( N FT R .EQ.3 ) GO TO 139 
SRX -0.98 

NCTRI =NBTRI 


SHM027 Til.'; 
SHM02710 
SHM0272i 
SHMO 2730 
SHM‘,274! 
SHMO 2750 
S HMD 2? 60 
SHMO 2770 
SHMO 27 80 
SHM027 90 
SUM 12301 
SHMO 28 10 
SHMO 2 820 
SHMO 2830 
SHM02840 ' 

SHMO 28 50 
SHMO 2 8 60 
SHMO 2870 
SHM02880 
SHM02890 
SHM02900 
SHM0291C 
SHM02920 


CSRX -SRX 
PSRX **5RX 

CALL RECESS(DPSF,DPY,BFRC» IFRC, AREA, RSBO,RSBIF,NRS»RSBBF) 

I PC IFRC.GE.0. 3) GO TO 139 
WRITEC 6,6 ) I F PC 

» FORMATC/lOX, • REJECTED IFRC =«,F8.4) 

IFRC =0.1 
BIVF =0.0 

L39 IF{ NFTR.GE.2) READ 103, L ZC , BM IR, SUZC, ETLF , BUZC, $1 AC,B1 VF, LZS 
I F( NFTR.EQ.3 3 READ 540, CSRX ,F SRX, NCTRI ,CHCAP , I FRC ,BFRC 
>40 FORMAT (2F10.5 ,I5,3F10.5) 

LAO BFHRC =BFRC** (1.0/24.0) 

BFRL =• ALOGCBFHRC) 

CALL FIXTRHCTRI , BTR I, NBTRI , NCTR I) 

TRIP =NFTR 

srx =cspx • < \ : . \ 

KHYD =1 '■ ' ' ■■ * " ' V ' 

LSHP = .FALSE • 

DO 141 K I A = 1,5 ■ , \ , 

KPSH(KIA) =0 

141 HBF(KIA) =D.O 

POINT OF RETURN FOR NEW TRIP 

142 IFIKRC .LE. 5) FTX =1.0 

IF COPY .EQ . 366) ME DWY ( 5 ) =366 

pph 1 - ■ SMf’fl 'V 1 ■ 


IF COPY .EQ . 366) ME DWY ( 1 

PPH =1.0 

IFC. NOT. LRC ) PPH =3.0 
IFCTR1P .NE. 1) PPH= 4.0 
I PPH -PPH 
FHPP -1.0/PPH 
IFPRC -IFRC** { FHPP/ 24*0 ) 
I FRL *” ALOGC I FPRC ) 


i ff Ws 




SHM02930 
SHM0294© 
SHMO 2950 
SHMU2960 
SHMO 2970 
SHM0298C 
SHM02990 
SHMO3O00 
SHM03* 10 
SHM03020 
SHMCJ 36 30 
SHMO 3040 
SHMfiJGS© 

SNM63060 

SHMO30T0 

SHMO308O 

SHMO 30 90 

SHMOBXCO 

SHM0311S 

SHMG3120 

SHM0313C 

SHM03140 

SHM03150 

SHMO 31 60 

SHMO3170 

SHM03160 

SHM03190 

5HM03200 

SHM93210 

SHMO 3221 

SHM03230 


•ECU 0 ) NCTRH =(NCTRI - 1 ) /4 +1 
HOURS FOR ADJUSTING HYDROGRAPH SHAPE VARIABLES 
C .AND* TRIP • EQ. 2) CALL STRHRS < RHPD, RHPH, I DY8, 
IHRE»NHPT*MXTRH,DPY»NRHP» I BTPR) 


VINTCR =FHPP*VINTMR 
NCTRH “NCTRI 
IF ( CONOPT ( 2) 

DETfc RHINE STORM 
I FINRHP • NE. 

1 I DYE, I HRBi 
USE *0.0 
NRTRI «,) 

PH A I *0.0 
SP IF *0.0 
OFUS-0.D 
OFUSIS *0.0 
RHFO-CUO 
URHF *0.0 
AM IF *0.0 
AMNET *0*0 
AMPET *0.6 
AMPRftC *0.0 
AMBF *6.0 
AMSE *0.0 
KRS *1 
KDRS *405 
UZS *0.0 
IPS; *0.0 

IFINFTR • G£* 2) GO TO 
IFtKRG .N£. 1) GO TO 

-■ BYL2S *6.00 
LZS *B YLZ S 
GO TO 145 

143 I F ( EML ZS ( 11) .LT. 

LZS =LZS*LZC/ PLZC 
IF(LLZC) LZS * 

I F ( ABS (FTX -1.0) . 

LZS* FTX*BBYLZS*LZC/BLZC 
IFtLRC .AND. ( LZC-LZS )*LT< 

144 IFCTRIP .EQ. 3 .OR. KFFC 
KFFC*0 

145 0CT1BF = f>.D5*TMRTF(D 
IF(DRSF(274) .LT* Q.05*TMRTF t 1 > I0CT1BF = 

I F C DRS F( 276 ) .LT. 0CT1BF*BFRC**2 ) 0CT1BF 
BYGWS = 0CT1B F/{ WCFS*BFRL*SQRT{ BFRC ) ) 

GWS =BYGWS "■ ■ -' " 

8YLZS *LZS ’ 

BFNX =GWS*BFRL 
TFCFS *BFNX*WCFS 

WRITE {6, 7) TRIP »LZC ,BM1R,SUZC,ETLF,BUZC»SIAC,BI VF,BFRC, JFR'C, 

1 CSRX » FSRXf NCTRI ,CHCAP 

7 F0RMAT(1HI,3X, 'TRIAL RUN NUMBER* , I3/5X» ’PARAMETER VALUES ’/10X, 

1 5HLZC = ,3X, F8.4»2X,6HBMIR =, 2X, F8.4»2X,6HSUZC *,2X»F8.4,2X, 

2 6HETLF = , 2X»F8.4, 2X, 6H8UZC =, 2X»F8.4,2X,6HSIAC =,2X,F8*4/10X , 

3 6HBIVF =»2X»F8«4»2X, 6HBFRC =»2X» F8.4*2X, 6HIFRC *, 2X>F8.4, ZX, 

4 6HCSRX *,2X,F8.4, 2X,6HFSRX = » 2X, F8.4,2X,7HNCTRI = , 1X,I8/10X, 

5 7HCHCAP =, lXf F8.0 ) 

WRIT*-: (5,3) LZS, GWS 


145 

143 


LZS) LZS =6MLZS( 11 ) 

LZC~ ( LZC- LZS ) *{SATFV/RATFV 
LT . 0.02) GO TO 144 


2.0) LZC = LZS + 2.0 
, EQ* 1) LZS = BBYLZS 


) 


DRSFC274) 

= DRSF(276>/BFRC**2 


S HMD 2 24*.’ 

SHMU325L 

SHM(, 3261- 

SH MO 3 270 

SHM03280 

SHM0329C 

SHMC'3300 

S HM03r IT 

SHM0332f- 

SHM0333i 

SHM0334C 

S H MO 3356 

SHM03360 

SHM03370 

SHM03380 

SHM633 9C 

SBM0340C' 

SHM03410 

S HMD 34 20 

SHM03430 

SH MO 34 40 

SHM03450 

SHM03460 

SHM03470 

SHM03480 

SHM03490 

SHM03500 

SHM0 1 

SHM03520 

SWH03530 

SH MOB 540 

SHM0355© 

SHM03560 

5HM03570 

SHM03360 

SWM03390 

SHMO36O0 

SHK03610 

SHM03620 

SHN03630 

SHH03640 

SHM^3650 

SHM, 03660 

SHM03670 

SHM03680 

SHM03690 

SHM0370C 

SHM03710 

SHM63720 

SHM03730 

SH MO 37 40 

SHMU3750 

SHMC’3760 

SH MO 3 7 70 


6 


> FORMAT t/5X, ' INITIAL MOISTURE STORAGES, 
1 F9.4) 

AETX =24.0*EPAET/365.0 
AF; 96 = 1. 2*AETX 
A EX 90 =0.3*AETX 
SIAM *1 m 2. ♦♦ S f A(T ' 

UZC»SUZC*AFX9C+8UZC*EXP{"2*7 *L ZS/LZC ) 
IF(UZC.LT.O*25)UZC=0.25 
MONTH *i ' 

MDAY * 273 

I Ff TRIP .EQ. It GO TO 147 
.46 WRITE t6» 9) ( TITLE (KTAl » KTA=1,20) 

9 FORMAT! 2 5X» 2CA4 I 

WRITE! 6, 10) YR1,YR2 

19 FORMAT (03X»6iH3PTI MI ZATIGN OF MODEL If 
1R YEAR 19,12, 1H» , 12) 

WRITE (6, ID 
11 FORMAT (5H JUNE) 

BEGIN DAY LOOP 
,47 DAY *274 
,48 CONTINUE 

I F( TRIP *NE. 1) GO TO 149 
KDRS * KDRS +1 

1 Ft RSBD( KRS ) .ME* DAY) GO TO 149 


LZS = »,F9. 


* GWS 


MODEL INPUT PARAMETERS BASED 


DAY) 


KDRS *1 
KRS *KRS +1 
9 CONTINUE 

ADIF =0.0 
A DBF *0.0 
TDSF *0.0 
PET =DPET (DAY) 

I F ( CONOP T < 1 ) oNEo 2) f 
PETU =P£T 
TFMAX =0.0 
DO 190 HOUR, =1,24 
IFtTRIP .NF. 2) GO TO 
LOGICAL VARIABLE *LSHP* 


P ET*EPCM( MONTH ) 


TROE DURING DURATION OF RECORDED 


= », SHMi?37 80 
S HMDS 7 90 
SHM» 330( 
S HMD 3 3 10 
SHM03820 
SH MOB 330 
SHM0384C 
SHMOBSSli 
SHM0386C 
SH MO 3 870 
SHM0388G 
SHM03890 
SHM0390C 
SHM0391G 
N WATESHMO3920 
SHM0393G 
SHM03940 
SHM03950 
SHMO'3960 
SHM03970 
SHM03930 
SHMG3990 
SHM04000 
SHMO4O10 
SHM04O2O 
SHMO4030 
SHM04040 
SHMO403O 
SHM04060 
^ Stw 4 7 . 
SHMO408O 
SHM04'*9t'- 
SH MO 4 100 
SHMG4110 
SHM- 412 
SHM04130 
HYDR0-SHM0414C 


GRAPH SD SYNTHESIZED DATA MAY BE SAVED DURING 
IFtKHYD .GT.ARHP) GO TO 152 
I F { 1 DYB(KHYD) .EQ. DAY .AND. IHRB(KHYD) .EQ 
IFtKHYD «GE. NR HP) GO TO 150 
I F ( IDY&tKHYD+l) .EQ. DAY .AND. IHRBtKHYD+l) 

1 KHYD+1 

159 I Ft I DYE(KHYD) .NE. DAY .OR. I HRE ( KHYD) .NE. 

KHYD = KHYD + 1 
LSHP = .FALSE. 

151 IFt .NOT. LSFP) GO TO 152 
KPSHt KHYD )= KPSH(KHYD) +1 
I Ft KPSHt KHYD ) .LT. 171 ) GO TO 152 

WRI TF ( 6 » 12 ) 

12 FORMAT (5? X#*FL COD HYDROGRAPH ARRAY EXCEEDED, SHORTEN NHPT 
1T0 HOURLY ROUTING* ) 

GO TO 22b 

152 CONTINUE 


SAVED DURING CORRESPONDING PERIOD 


HOUR) LSHP 


► TRUE, 


■ HOUR) KHYD= 


HOUR) 


SHIFT 


SHM04150 

SHM04160 

SHM0417C 

SHMO4180 

SHMU4190 

SHMC4200 

SHM04210 

SBM*‘J4220 

SHMO4230 

SHMG424C 

SHM04250 

SHM04260 

SHM04270 

SHM0428G 

’SHMO4290 

SHM043OL 1 

SHM0431O 


IFUNSGRD .EC.C ) .AND. ( DRHBt DAY, HOUR) .NE. 0.0) .AND. (PET .EG. 

I PE'U) .AND. (CONOPT ( 1 ) .NE. 0)) PET = 0.5*P£T 
15:' I F { HOUR .60. SGRT+1) RGPM =DRGPM{ DAY) 

1FCH0UR „6Q. 9) HSE = ( FWTR*PET 1/12.0 

I F( HOUR .FQ. 21) HSE =0.0 
PRH =RGPM*DRhP{ DAY, HOUR) 

AMPREC -AMPREC + PRH 
ARHF =0.0 

C 15 MIN ACCOUTI NG AND ROUTING LOOP (60 MINUTES USED FOR ROUGH 
C ADJUSTMENT , AND 20 MINUTES FOR FINE ADGUSTMENT IN TRIP 1) 

DO 182 PRD= 1,IPPH 

I F ( L S HP .AND.CONOPT (2) »EQ® O.AND. PRD .NE. 1) KPSH(KHYD) = 

1 KPSHt KHYD ) +1 
PEB1 =0.0 
PPI *0*0 
OFR -0.0 
OFRIS =0.0 
WI = 0.0 

woifs =&.© 

PEP =FHPP*PRH 

IFtTRIP *GE. 2 • AND.LNPR) CALL PREPRD ( RGPM, DRHP, DAY, HOUR »DPY, PRD, 


SHM- A’-. 20 
SHM0433C’ 
SHM *43 4C 
SHM04350 
SHM04360 
SH Mi, 4 3 70 
SHM04S 80 
SHMu43 90 
SHM0440© 
SHM044 10 
SH MO 44 20 
SHM04430 
SHM04440 
SHM04450 
SHM0446O 
SHM04470 
SHM0448C* 
SHM04490 
SHMO45O0 
SHM045 10 
SHM04520 


1 PEP, PRH) 

I Ft PEP . GT. 0.0) GO TO 155 
I F t OF US .GT. 0.0) GO TO 15T 
I F t IF S .GT. C.P) GO TO 167 
; -IFtTRIP .FQ. 1) GO TO 181 
I Ft NRTRI .GT. D) GO TO 169 
TRHF = 0.0 

I F t • NOT. LSHP) GO TO 154 
KHPT =KP SH{ KHYD ) 

SSRtKHYD,KHPT) =0.0 
154 CONTINUE 

IFtPHFC .GT. 0.0) GO TO 178 
GO TO 181 

C RAINFALL UPPER ZONE INTERACTION 

155 IFtPEP . GE.V INTCR) 60 TO 156 
UZS =UZS +PE F*TPLR 

VINTCR = V INTCR- PEP. 

PPI =0oD 
PEB 1=0.0 

IFtOFUS .GT. C.D ) GO TO 157 

GO TO 167 

156 PPI = PE P-- VINTCR 

UZ-S =UZS+VI NTCR*TPLR 
VINTCR = 0.0 

LZSR =LZS/LZ C 

UZC =SUZC*AEX90 +BUZC*£XP {•» 2.7*LZSR ) 

IFtUZC .LT. C. 25) UZC =0.25 

UZRX= 2 .0*AB St UZS/ UZC- 1.0 ) +1.0 

FMP =(1.0/(1.0 +UZRX) ) **UZRX 

I Ft UZS .GT. LZC ) FMR =1.0-FMR 

PEBI=PPI#FMR 

UZS =UZS+PP I— PEB I 

1C LOWER ZONE AND GROUNDWATER INFILTRATION 


SHM04530 
SHM0454O 
SHM04550 
SHM04560 
&HM0457O 
SHM04580 
SHM04590 
. SHM04600 
SHtKMftlO 
SHMt/4623 
SHMG4630 
SHM04640 
SHM04650 
SHM04660 
SHM04670 
SHMU4680 
SHM04690 
SHM04700 
SHM04716 
SHM04720 
SHW0472M3 
SHMS4740 
SHM04750 
SHM<H76Ci 
SBM0A77C 
5HM0478O 
SMM0479O 
SHM04800 
SHM04810 
SHM<>4820 
SHM0483O 
SMM0484C 
SH MO 48 50 


157 


LZSR =L Z S/LZC 
FID =4© 0*LZ5R 
I F( LZSR .LG. 1.0) GO TO 158 
E ID =4.o+ 2. C*( LZSR®1. 0 ) 

I F ( L ZSR ©LE.2.D) GO TO 158 
HID* 6© 0 

CMIR =FHPP*SIAM*BMIR/{ 2.Q**EID) 

158 PEBI **PEBI +CFUS 

Cl VM *BI VF*2 *0**LZ$R 
IFtCIVM .LT. 1.0 )C1 VM *1.0 
PEAI =PEBI*PEBI/(2.C*CMIR*CIVM) 

WI » P6BI*PEBI/(2.0*CMIR) 

I F { Pi!B I ®GE. CMIR) Wf= PFBI-C.5*CMIR 
I F C P2=BI ©GE® CMIR*C IVM) PEAI =PEBI -Q. 5*CMI R*C I VM 
WE1FS *WI~PFAI 

IF ( C PEAI -OFUS ) .GT.O.O) GO TO 159 

EQD » (OFUS + PEAU/2.0 
GO TO HO , ■ z'.' 1 

159 EQD ■ 6QDF*( C FEAI-OFUS >**0.6 ) 

160 I F t C OFUS + PEa I > .GT. (2.0*EQD> ) EQD =0.5* (OFUS + PEAI ) 

IF( (OFUS+PEAI ) ®LE. 0.001) GO TO 161 

QFR=FHPP*OFRF*{{ (OFUS +PEAI ) *0.5) **1.67 )*( ( 1.0 + .6*( (OFUS + 

1 PEAI )/( 2.0*EQD) )**3.Q)**1.67) 

IF(0FR .GT. ( G.75*PEAI ) ) OFR =0®75*PEAI 

161 IFCFIMP .EQ® C.O) GO TO 165 

162 PEIS =PPI +OFUSIS 

f, Ft I PE IS* OFUS IS) .GT.O.O) GO TO 1*3 
EQDIS *( D FUS I S + PEIS)/2.0 
60 TO 164 

163 EQDIS* EQDFIS*((PEIS -OFUSIS) **0.6) 

164 I F ( ( OF US I S+P6IS) .GT. (2.0*EQDIS>) EQDIS =0.5*(0FUSI S +PEIS) 

I F ( ( OF US I S +PEIS ) .LE.0.01) GO TO 165 

OFR IS =FHPP*OFRFIS*( ( ( OFUSIS +PEI $)*0.5 )**1.67 ) * ( ( 1 .0 +©.6*( { 
i OFIJSIS +PEI S) / (2.0*EQDFIS) }** 3.0 ) **1.67 ) 

IF ( OFR IS . GT* PEIS) OFRIS =PEI$ 

163 OFUSIS =P E IS- CFR IS 
OFUS = PF A I -OFR 
IF(0FUS.GE.0*C01 )G0 TO 166 
LZS =LZS +OFUS 
OFUS=Qoii 

OFRIS =OFR IS 40FUSIS 
0FUSIS=C*0 

166 LZRX =1.5*ABS ILZS/LZC *1.0) +1.0 
FMR =(1®Q/(1.C +LZRX ) ) **LZRX 

IFCLZS • LT. LZC) FMR = 1.0- FMR* (LZS/LZC > 

PLZS=FMR* ( PGB I*=’WI ) 

PGW=(1© -FMR) *( PEBI ~WI)*( 1.0 -SUBWF J*Ff»ER 

GWS»GWS+PGW 

LZS =LZ$ +PLZS 

IFS=IFS+WEIFS*FP£R 

167 SPIF =IFRL*IFS 
AMI F =AM I F+SP IF 
ADI F*AOI F+SPI F 

I FS «IFS“SPIF 
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SHM04910 
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SH MO 49 30 
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SH MU 4 950 

SHMQ4960 

SHM04970 
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SHM04990 

SHM05000 

S HMD 50 10 

SHM05020 

SHM05030 

SHM0504D 

SHM05050 

SHM05060 

SHM05070 

SH MO 50 80 

SHM05090 

SHM05100 

S HMDS 11© 

SHMO5120 

S HMD 51 30 

SHM05140 

SHH05150 

SHF0S160 

SHM05170 

SHM05180 

SHMQS19C- 

SHM05200 

SHMOS210 

SHMU5220 

SHM05230 

SHM05240 

SHM05250 

SHM05260 

SHM0S270 

SHM05280 

SHM05290 

SHM0530C 

SHMO5310 

SHM0532D 

SHMCS3I& 

SHM05340 

SHMO5350 

SHMC5360 

SH MO 5375 

SHM0538D 

SHM05390 



168 


+FIMP*QFRI S +SPIF 


170 


INTER FLOW ENTERING STREAM 


=4»Q*URHF*WCFS 

=4.Q*UHFA(1)*WCFS 


IF( I FS.GE. 0.0001) GO TO 168 
LZS=LZS+I FS 

IFS=>' o ' 

LHFA ( 1 ) = FPER*OFR +PPI*FWTR 
IFCTRIP .ME. 1) 60 TO 169 • 

ARHF = ARHF +UFFA { 1 ) 

GO TO 181 
1 ROUTING 

169 I F ( CONOPT (21 .ME. 1) GO TO 
URHF =URHF +0 .25*UHF A( 1 ) 

I F( PRD.NE. 4) GO TO 178 
UHFA(1)=URHF 

C SAVES SYNTHESIZED DIRECT RUNOFF AND 
C DURATION OF RECORDED HYDROGRAPH 

170 I F ( .NOT. LSHP) GO TO 171 
KHPT *KPSH{ KH YD) 

I F ! CONOPT (2) .EQ.U SSRCKHYD ,KHPT) 

IFCCONOPT (2) *EQ. 0) SSR(KHYD»KHPT) 

CONTINUE ' . 

TRHF « 0.0 - W'.r- '• ■ 

KTRI •NCTRI 
URHF *UHFA( KTRI ) 

IFIURHF «LE. €.0) GO TO 1T4 
TRHF -TRHF +URHF*CTRI {KTRI ) 

UHFAiKTRI + II * URHF 
GO TO 175 

UH FA IKTRI +1) =D*0 
KTRI = KTR I -1 
IF (KTRI .GE. 1) GO TO 172 
IF (URHF oLE. C.O) GO TO 177 
NRTRI =NCTRI 
NRTRI *NR TRI -1 
UHFA(l) =0.0 
URHF *0. D 

I Ft TRI P • LE* 2 ) GO TO 179 
I F ( TFCFS .LE. 0. 5*CHCAP ) SRX = CSRX 

I F < ( TFCFS .GT. 0.5*CHCAP) .AND, (TFCFS . LT. 2.©*CHCAP)) 

1 +( FSRX -CSRX)* UTFCFS»0«5*CHCAP)/(1.5*CHCAP) )**3 

IFtTFCFS .GT. 2.0*CHCAP ) SRX =F$RX 
RHF1 =TRHF-SRX*(TRHF~RHFG) 

p Li C A -OUC 7 

I F( RHFO.LT « RHFMC) RHFO =0,0 
TFCFS =( 4.04R HF1 +CBF-HSE )*WCFS 
IFCTFCFS .LE. IFMAX) GO TO 180 
PRDF =PRD 
TDFP24 =HOUR 

I F ( PRO • LE.3 ) TDFP24=( TDFP24-1.C) 

TFMAX =‘i FCFS 
ARHF “ARHF + BHF1 

I F ( V INTCR .LT. FHPP* VINTMR ) VINTCR =VINTCR +DPET (DAY )/ t 24.0/ 
1 FHPP ) , ■ 

182 CONTIMUt 

C END OF 15 MINUTE LOOP 
C ADDING GROUNDERWATER FLOW 


171 


172 

173 


174 

175 

176 

177 


178 


SRX 


179 


+ 0,15*PRDF 


180 

181 


SH MG 5*6 00 
SH MO: 54 It 

SHM© 5420 
SHM054 31 
SHM0544:" 
SHM05450 
SHM0546G 
SH MO 5474 
SHM05480 
SHM0549C 
SHM05500 
SHM05510 
DURINGSHM05520 
SHM0553D 
SHM05540 
SHM0555C 
SHM05560 
SHM0557C 
SHM05580 
SHM05590 
SH MO 5 600 
SHM05610 
SH MO 5 6 20 
SHM05630 
SHM05640 
SHM05650 
SH«f66& 
SHM 56 7i 
SHM05680 
SHM0569C 
SHMO57O0 
SHM05710 
SHMO5720 
SHMO5730 
SHM05740 
SHM0#fS© 
= C SRX SHM© 5766 
5HM05770 
SHM15780 
SHMO5T90 
SHM© 5 800 
SHM© 58 1C 
5 HWJS 82 CI 
SHM0583O 
SHM05840 
SHM© 5 850 
SHM05860 
SHM05870 
SHMQ5880 
SHM© 5 ' 85)0 

SHHt :>9C { 
SHM© 5910 
SHM0592C- 
SHM05930 



183 CBF =GWS*BFRL 

IFtKHYD .GT® NRHP ) GO TO 164 

IFCLSHP .AND. (HBF(KHYD) .EQ. O.O) ) HBF(KHYD) 

184 GWS =GWS"CBF 
AMBF “AM8F +CBF 
THGR “ARHF +CEF 

: EVAPORATION FROM STREAM SURFACE 

185 IFIHSr, .GT. ThGR) HSE =THGR 

I F C CBF . GT. HSE) ADBF *ADBF +CBF -HSE 

AMSE »AMSE+H5fS 

THSF (HOUR ) « {THGR ~HSE)*WCF5 
IFCTFMAX • L6. 0*0 ) TFMAX =THSF (HOUR ) 

TDSF * TDSF + THSF (HOUR) 

C DRAINING OF UPPER ZONE STORAGE 
UZINFX =tUZS/UZC) - ( LZS/LZC) 

IF (UZINFX • LE • 0.0) GO TO 186 
LZSR «LZ$/LZC 

UZINLZ «f)®003*6MlR*UZC*UZINFX**3.0 

IFIUZ1NLZ .GT.UZS) UZINLZ = UZS 

UZS *UZ$ -UZINLZ 

LZRX < a l»5#ABS (LZSR ~1.0>+1.0 

FMR «{ X.O/C hO + LZRX ) ) **LZRX 

IFILZS .LT. LlC) FHR=1.0-FMR*LZSR 

PGW “U.3-PMR )*UZINLZ* { 1.0 -SUBWF ) *FPER 

PLZ5*FMR*UZINU 

LZfCeLZS +PLZS 

GWS =GWS +PGW 

C 4 PM ADJUSTMENTS OF VARIUS VALUES 

186 I F{ HOUR .NF. 16) GO TO 190 
A EX 90 *0* 9*{ AEX90 +PET ) 

AFX96 =0. 96* ( AEX96 +PET) 

C INFILTRATION CORRECTION 

SIAM = (AEX96/AETX) ♦♦SI AC 
I F ( SI AM .LT. €.33) SIAM =0.33 
IF { PET .EQ. 0.0) GO TO 190 
: EVAP-TRANS LOSS FROM GROUND WATER 

GWET =GWS*GWE7F*PET*FPER 
GWS =6WS“GWRT 
AMPgT =AMPFT +PET 
I F( PUT o GE. UZS) GO TO 187 
UZS =UZS -PET 
AMNET =AMNET -+PET 
GO TO 19D 

187 PET = PIT- UZS 
AMNFT =AMNET 4 UZS 
UZS =0.1 

LZSR = LZS/LZC 

IF( PFT . GE. E 7LF*LZSR) GO TO 188 

SET = PET* ( 1. 0- P£T/(2.0*ETLF*LZSR) ) 

GO TO lb9 

188 SET *0 © 5*TTL F*L ZSR 

189 LZS =LZS-SET 

AMNET =AHNFT +SET 

190 CONTINUE 


CBF*WCFS 


SHMC5 940 

SH MO 5 950 

SHM0596I; 

SH MO 5970 

SHM05 98C- 

SH MO 5 990 

SHM06000 

SHM06G1C 

SHM06020 

SHM06f-3C 

SH MO 60 40 

SHM0605I 

SHMf"‘6<;60 

SHM-. 6- IQ 

SH MO 60 80 

SHM06090 

SHM06100 

SHM06110 

SHM06120 

SHM06130 

SHM06140 

SHM06150 

SHM06160 

SHM06170 

SH MO 6 180 

SHM06190 

SHM06200 

SHM06210 

SHM06220 

SHM06230 

SHM0424O 

SHM06250 

SHM06260 

SHM0627O 

SHHO6280 

SHM0629G 

SHM06300 

SHM06310 

&HMQ6320 

SMM06330 

SHM06344) 

SHMO6350 

SHM06360 

SH MO 637 6 

SHM06380 

SHM06390 

SHM06400 

SHM06410 

SHMO6420 

SHM06430 

SHM06440 

SHM0643G 

SHM0646O 

SHM06470 



c 


END OF HOUR LOOP 

DSSF(OAY) =TD SF/ 24*0 
I Ft TRIP# HQ* 1 ) GO TO 192 
IFtTFMAX © LG© RMPF) GO TO 192 
I Ft DAY *EQ. 366) MDAY =337 
DATE “MOOt DAY, MDAY) 

WRIT?" t 6 , 13 ) DATE, ( THSF {HOUR ) ,H0UR=1 » 12 ) 

13 FORMAT tlH/ » IX / ,1X , I4,2X ,2HAM, IX ,6F8«1»3X,6F8.1) 

WRIT;H6, 14) { THSF (HO UR) ,H0UR=13,24) » DSSFtOAY) 

14 FORMAT (lH/,6X,2HPM,lX,6F8ol,3X,7F8*l ) 

I F ( T0FP24 «LT • 12.0) GO TO 191 
TDFP12 “TDFP24 -12.0 

WR IT2 1 6, 15 ) TFMAX* TDFP12 

15 FORMAT 1 1H/, 10X, 1 MAXI MUH» *, F8. 1, 2X»*C.F. S. « ,5X,' TIME* ,3X»F5. 2, 2X , 

1 4HP.M. ) 

GO TO 192 

191 WRITE { 6 , 16) TFMAX, TOFP24 

16 FORMAT (l H / , 10X,8HMAX IMUM=, F8.1 » 2X , 6 HC. F* S* ,5X,4HTIMf:,3X,F5*2,2X, 

1 4HA. M. I ' 

192 CONTINUE 

IFCTRIP .EG. 1 .AND* .NOT. LRC .AND. KDRS .L 6 . 3 .AND. I FRC .GT, 
1 O.i) SIFRS(KDRSfKRS-l) = ADIF*VWIN 
IFtTRIP 06 Q 0 1 .AND. KDRS • LIE. 3) SBFRSt KDRS ,KRS-1 ) =ADBF*VW IN 

MONTHLY SUMMARY STORAGE 

IF (DAY .NE. MEDWYt MONTH)) GO TO 206 
TMPREC (MONTH) = AMPREC 
AMPREC *»D.O 
* TMBFt MONTH) * AMBF 
AM 8 F =0.0 

TMI F ( MONTH ) = AMI F 

AM I F =0.0 

TMSE ( MONTH) “AMSE 

AMSO =0.9 

TMPET t MONTH) = AMPET 

AM PET =D • 0 

TMNETt MONTH) =AMNET 

AMNFT =0.0 

EMGWSl MONTH) = GWS 

UZC=SUZC*AFX9C +BUZC*EXP(-2. 7*LZS/LZC) 

IFtUZC .LT. 0.25) UZC =0.25 

EMUZC(MONTH)=LZC 

EMUZS( MONTH ) =UZS 

EMSI AM(MONTH) “SIAM 

EMLZSt MONTH) =LZS 

EMI FS(MDNTH) =IFS 

I F( MONTH .EG. 5) MEDWYtS) =59 

MDAY “MEDWYt MCNTH) 

IFtTRIP . EQ. 1) GO TO 205 

193 GO TO (194,195,196,197, 198, 199,200,201,202,203,204, 

1 2C 5 ), MONTH 

194 WRITE t 6, 17 ) 

17 FORMAT ( 1H/ , 4 H JULY) 

GO TO 215 

195 WRITE(6,18) 


SH MG 6481. 

S HMD 64 96 
SH MO 6 6 GO 
SHM06510 
SHMU652C 
SHM0653C 
SHM06540 
SHM0655C 
SHM0656& 
SHM0657G 
SHMC6580 
SHMO6590 
SHMC6606 
SHM0661© 
SHM06620 
SHM06630 
SHM06640 
SHM06650 
SHM0666C 
SHM06670 
SHM06680 
SHM06690 
SHM0670© 
SHM06710 
SHMQ6720 
SHM06730 
SHM06t40 
SHM06750 
SHM06760 
SHM06770 
SMMG6780 
SHM06790 
SHM06800 
SHM06810 
SHM06820 
SHM048$0 
SHM$fc840 
SHMQ6850 
SHM068AC 
SHM0A87O 
SHM06S80 
5HMO6890 
SHnmmm 
SHM06913 
SHM0692D 
SHM0693D 
SHM06940 
5HM0695C 
SHM0&96D 
SHM06970 

snmmm 

SHMD6990 

SHM07000 

SHMG7010 



18 FORMAT !1H/,6HAUSUST) 

GO TO 205 

196 WRXTc(6,19) 

19 FORMAT ( 1 H /, 9 H SEPTEMBER) 

GO TO 2.5 

197 WRITE (6, 20 ) 

20 FORM AT ( 1H/,7HCCT0BER ) 

GO TO 2D 5 

198 WRITL16, 21 ) 

21 FORMAT (1H/,8HN0V EMBER) 

60 TO 205 

199 WRITER, 22) 

22 FORMAT (1H/»8H DECEMBER) 

GO TO 2D 5 

20 f WRITE 16» 23) 

23 FORMAT (1H/,7H JANUARY) 

GO TO 2D 5 

201 WRITE (6,24) .«■ 

24 FORMAT UM/,8RFEBRUARY) 

GO TO 2S5 • 

202 WRITEC6.25) ' 

25 FORMAT ClH/» 5HMARCH) 

GO TO 205 

203 WRITE! 6, 26) 

26 FORMAT Cl H/»5HAPRIL) 
i ' ' GO TO 205 

204 WRITE! 6, 27) 

27 FORMAT (1H/»3HMAY ) 

205 MONTH -MONTH +1 
C END OF DAY LOOP 

206 CALL DAYMYT! D AY» DPYJ 

IF! DAY oNE.274) GO TO 148 
I F ( TRI P .NE. 2) GO TO 208 

C ADJUST BASF. FLOW FOR AVERAGE VALUE DURING STORM 
i IFtNRHP » FQo 0) GO TO 208 

DO 207 KH YD «l f NRHP ■■■■■' v 

DAY =1 DYB ( KHY E) 

IF ( DSSF! DAY) .GT. HBF(KHYD) ) GO TO 207 
HBF(KHYD) = ! H EF( KHYD ) +DSSF! DAY) )/2*0 

207 CONTINUE 

208 I F ( TRI P . NE* 1) WRITE(6,28) CTITLE(KTA) , KTA=l,2Q»l) 

28 FOR. MAT ( 1H 1 , 25 >» 20A4 ) 

C ANNUAL SUMMARY 
APREC =0.0 
ft 8 FV =0.0 
ASH V **0*0 
A NET =0*0 
APFT =0*0 
AIFV = 0*0 
DO 209 MONTH =1,12 

apr:;c =aprec +tmprec (month) 

ABFV=ABFV +TMEFC MONTH) 

ASEV=ASE V+- TM SEC MONTH) 

ANET =ANET +TMNET! MONTH) 


SHMQ702© 

SHM0703© 

SHMU704R 

SHM070 50 

SHH07C60 

SHM07C70 

SHM07080 

SHM0709C 

SHM07100 

SHMU711U 

SHM07120 

SHM07130 

SH MOT 140 

SHM07150 

SHM07160 

SHM07170 

SHMU7180 

SH MOT 190 

SH MO 7200 

SHM07210 

SHM07220 

SHM07230 

SHM0724© 

SHM07250 

SHM07260 

SHM07270 

SHMO7280 

SHM* 729* 

SHM07300 

SHM07310 

SHMf-7320 

SBM07330 

5HM07340 

SHM07350 

SHM07360 

SHM0737© 

SHMK738© 

SHM07390 

SHM07400 

SHM07410 

SWM07420 

SH MO 7430 

SHM07440 

SHMG7450 

SHM07460 

SHM©7470 

SHM07480 

SHMQ7490 

SHMO750O 

SHM0T510 

SHM0T520 

SHM07530 

SHM07540 

SHM0755© 



APET =APHT +TM PET (MONTH) 

2C9 AIFV=AIFV +TMIFC MONTH) 

WR I T(; ( 6* 29 ) 

29 FORMAT ( 1H// /44X, 23HS YNTHESIZED FLOWS) 

21! IFITRIP • EQ. 1) WRITE! 6t 30 ) 

:'C FORMAT ( //5X* * SUMMARY WHILE OPTIMIZING VOLUME VARIABLES*) 
211 CALL DAY$UM(D$$F,MEOCY,0PY,SATFV,TMSTF ) 

I F( TRIP .EO. 1) GO TO 212 
CALL DAYOUT ( D SSF ,MEDWY» OPY ) 

212 WRITS (6» 31 ) (TMSTF(KWD), KWD=1 , 12 ) » SATFV 
31 FORMAT (IX, 9H SYNTHET IC,3X, 12F8. 1 ,2X, FIO. 1 ,2X »3HSFD ) 

DO MONTH » T * 1 2 

213 TMSTFI (MONTH) = TMSTF (MONTH) /VWIN 
SATFVI *SATFV/VWIN 

WRITS (6, 32) (TMSTFI(KWD) , KWD=1 , 1 2) , SATFVI 


SHM0756L 
SHKu 757C 
SHMu7»er 
SHM0T39C 
SHMO76O0 
SHMO7610 
SHM0T620 
SHMQ7630 
SHM07640 
SHMO7650 
SHM07660 
5HM0767O 
SHM07630 
SH MO 76 90 
SHM07700 


3 2 FORMAT C IX, 5HT CTAL, 8X» 12F8. 3#4X * F7.3, 2X»6HINCHES ) 

00 214 MONTH =1,12 

TMOF (MONTH) «TMSTFI ( MONTH)-TMIF( MONTH )-TMBF( MONTH) + 

l TMSE (MONTH) 

.4 I F(TMOF( MONTH ) *LT. 0.0) TMOF( MONTH) =0.0 
AOFV=SATFVI**A IFV-ABFV +ASEV 
I F ( AOFV .LT.<H0) A0FV=0.0 
WRITE(6,33> ( TMOF(KWD) , KWD=1,12) , AOFV 

53 FORMAT (IX, * OVERLAND* , 5X» 12F8.3,4X,F7.3, 2X, * INCHES* ) 
WRITE(6,34) (TMIF(KWO), KWD= 1, 12) , AI FV 

54 F0RMAT(1X»9HI NTERFL0W*4X»12F8.3*4X, F7®3, 2X»6H INCHES ) 

WRITEf 6,35 ) (TMBF(KWO), KWD=1,12) ,ABFV 

55 FORMAT ( IX ,4HBASE ,9X,12F8.3,4X,F7.3,2X»6HINCHES) 

WRIT6(fc,36) C TMSE(KWD) , KWD=1,12), ASEV 

56 FORMAT (1X,9HSTRM EVAP,4X, 12F 8„3,4X, F7.3, 2X,6H1NCHES) 

WRITE- { 6, 37 ) ( TMPRECC KWD) , KWD=1 ,12 ) , APREC 

57 FORMAT ( lXf 6HPRECI P, 7X» 12F8. 2»3X,F8*2,2X,6HINCHES) 

WR I T C ( 6 , 3 8 ) ( TMNET ( KWD ) , KWD=1 » 12) , ANET 

58 FORMAT (IX , *EV, P/TRAN-NET * ,2X, 12F8® 3 »3X»F7.3*2X»* INCHES* ) 

WRI T& ( 6, 3 9 ) (TMPET(KWO), KWD=1,12) , APET 

59 FORMaTBX, 10H«-PDTENTIAL,2X,12F8.3,3X,F7.3»2X,6HINCHES) 
WRITr-;(6,40) ( EMU ZSIKWD) , KWD=1,12) 

tO FORMAT (1 X , 12H STORAGE S-UZS » 2X ,12F8. 3 , 1 2X, 6HINCHES > 

WRITE ( 6 , 41 ) ( EMLZS(KWD), KWD=1»12) 
kl FORMAT ( 10X ,3H LZS , 2X, 12F8.3, 1 2X, 6H INCHES) 

WRITFt 6,42 ) ( EMI FS(KWD), KWD = 1,12) 

42 FORMAT (10X, 3H IFS ,2X, 12F8.3* 1 2X,6HINCHES ) 

WRITE { 6, 43 ) ( EMGWS(KWD) , KWD=1,12) 

\3 FORMAT 1 10X f 3HGWS , 2X, 12F8.3 , 1 2X, 6HI NCHES) 

WRITE (6,44 ) ( EMUZC(KWD) , KWD=1,12) 

%4 FORMAT (IX, 12H IND ICES-UZC ,2X,12F8.3) 

WRITS (6,45) ( EMSIAM(KWD), KWD=1, 12) 

\5 FORMAT ( 9X , 4HS IAM, 2X, 12F8.3) 

AMBER =(LZS-BYLZS)*FPER+(UZS + IFS ♦ GWS - BYGWS>*(UO » f 
1 ) + SATFVI +ANET*FPER+ASEV-APREC , 

WR IIP ( 6, 46 ) AMBER 

<f6 FORMAT (1 H/ ,7H BALANCE ,5X,F10o4,2X,6HlNCHES ) 

ESTABLISH WHETHER MONTH IS PREDOMINATELY BASE FLOW OR DIRECT 

NOFM-O , ' ■' 


FWTR 


RUNOFF 


SHM07710 

SHM07720 

SH MO 7730 

SHM07740 

SHMU7750 

SHMG7760 

SHM07770 

SHMQ7780 

SHMO7790 

SHMO78O0 

SHM07810 

SHMC7820 

SHM07830 

SHM07840 

SHMC7850 

SHMU7860 

SHM07870 

SHM07880 

SHM07890 

SHM07900 

SHM07910 

SHM07920 

SHM0793B 

SHM07940 

5HMO7950 

SHM07960 

SHM07970 

SHM07980 

SHM0T990 

SHMO8OO0 

SHM0801G 

SHN08020 

SHM0803C 

SHM0804G 

5HM080S& 

SHM08S60 

SHH08070 

SHM08080 

SHM08U9G 



t-3 t_3 O 


M0NTH1=4 


215 


M0NTH1=1 

IFCFTX .LT. €.95) 

DO 216 MONTH-1,12 
XMPFT C MONTH ) *0.0 
I F( MONTH *LT* R0NTH1 ) GO TO 216 
I F( TMSTFl (MONTH ) .GT. €.OCi) GO TO 215 
XMPFT ( MONTH)® l.D 
GO TO 216 

I F ( TMBF( MONTH J/TMSTF I ( MONTH ) 

I F(TMOF( MONTH )/TMSTF I ( MONTH) 

NOFM “NOFM+1 


.GT . 
.LT, 


0, 

0, 


5) 

5) 


XMPFT ( MONTH) =1, 
GO TO 216 


Y 


MPFTCMONTH) =2.0 


SHMt.SlGC 

SHht-anv 

SHM**8120 

SHMi/8130 

SHM^814f; 
SHM08150 
SHMG8160 
SHMT'8170 
SHM08180 
SHM0819L 
SH MO 8 200 
SH MO 8 2 10 
SHM66220 
SHM08230 
SH MO 8 2 40 
SH MO 8 250 
SHM08260 
SHM08270 
SHM08280 
SHM08290 
SHM08300 
VALLESSHM08310 
SHM08320 
SHM08330 
SH MO 8 3 40 
SHNO63S0 


216 CONTINUE 
NATURE OF TRIPS 

TRIP 1 OPTIMIZE VOLUME VARIABLES WHILE BYPASSING ROUTING 
TRIP 2 SET FLOOD HYDROGRAPH VARIABLES. CSRX, FSRX » NCTR I » CHCAP 
TRIP 3 FINAL PUN WITH OPTIMIZED VALUES 

217 IFITRIP .EQ. 1) GO TO 218 
KRC *MNRC +1 

IFCTRIP .EQ. 2) GO TO 226 

GO TO 228 

C SYSTEMATIC ADJUSTMENT OF VOLUME VARIABLES CONVERGING ON OPTIMUM 

218 KRC* KRC +1 
K8RC =KBRC +1 
PLZC =LZC 
PBMIR =BM|R 
PSUZC *SUZC 

PgTLF -ETLF SHM0837Q 

PBUZC =BUZC SH.M08380 

PS T AC »SI AC > 5HPIQ0 , 3 ll 9 ! tlfc . 

C ADJUST FIVE VOLUME VARIABLES. L ZC, SUZC* ETLF , BUZC, SI AC SHM08400 

CALL SFTFVPIUZC, SUZC, ETLF, BUZC, SI AC, TMSTF,TMRTF,TMPREC,TMPET, SHM08410 

1 EMLZ S,S SQM, LRC , XMPFT, FTX, NOFM, LBUZC ,LETLF,LLZC, APREC,APET) SHM08420 

fC ADJUST INTERFLOW VOLUME CONSTANT DURING FINE ADJUSTMENT PHASE SHM0843® 

FNCTRH =NCTRH SHM08440 

IF (.NOT. LRC .AND. IFRC .GT. 0.1) CALL S ETB I V { B I VF , NRS , I FRC , R SB I F , SHM0845D 
1 SIFRS, FNCTRH) $HK08460 

C ADJUST INFILTRATION RATE CONSTANT 8MIR 
I Ft. NOT. LBMIR) GO TO 219 
BMIR =0.9*BMI R 
GO TO 220 

IF(ABS(FTX“1.C) .GT. 0.02 .AND, KRC. GT® 5) IFT =2 
CALL SFTBMI (BMIR ,NRS,BFRC,RS BBF , SBFRS, FNCTRH , I FT ) 

IF( (KRC .GT. 6 ) .AND. ( LZC . GT. 29.0) ) LLZC* .TRUE. 

IF( (KRC .GT. 6) .AND. (ETLF .GT. 0.59) ) LETLF * .TRUE. 

•GT. 3,9)) LBUZC* .TRUE. 


219 


220 


I F ( ( KRC .GT. 6) .AND. (BUZC 
I F ( .NOT. LLZC) GOTO 221 
LZC =PLZC*SAT FV/RATFV 
WRIT!. (6, 47) LZC 


47 FORMAT ( / 2X» ’ LZC WAS CHANGED T0NF6.2,* BASED ON ANNUAL RUNOFF 

iMf;* ) ■ ' ■ r-'M -'iK't. 

221 I F (KRC .LT. 6 
LBMIR * .TRUE; 

BMIR* 21.0 


.OR. BMIR .LT. 20.0) GO TO 222 


SHM08470 
SH MO 84 80 
'$HM0849© 
SHM08980 
SHM08510 
SH MO 85 20 
SHM08530 
SHM08540 
SHM08550 
SHM0856O 
SHMU857CI 
SHM0B58O. 

VCLU$HMQ8500 

SMMO'8600 

SHM0'8610 

■ : SHM;.,8 62 ; . 

SHM08630 



.GE. MNRC -AND* KBRC *GF. 2) .OR 


€•15 ©AND® LRC ) GO TO 223 


' f "Jf' 


GO TO 228 


222 !F{ SSSQM .LE. SSQM .AND. ( (KRC 

1 ( • NOT c LRC) )) GO TO 224 

I F ( SSSQM ,LF. SSQM) GO TO 142 
I F ( KRC .GE. MNRC) KRC-KRC-1 
BLZC-PLZC 
BBMI R-PBMIR 
B5UZC-PSUZC 
BXTLF- PETLF 
BBUZC-PBUZC 
BSIAC -PSIAC 
SSSQM- SSQM 
BBYLZS-BYLZS 
KBRC-0 
IFCSSQM .LT 
GO TO 142 
LRC-. FALSE. 

WRITE! 6,48) 

FORMATI/5X, ‘SHIFT TO FINE ADJUSTMENT BEGINNING AT BEST ROUGH 

1TMENT POINT ’ l 
SSSQM- 1000.0 
GO TO 225 
CONTINUE 

IFILRC) GO TO 223 

IFITRIP .GE. NLTR) 

TRIP -JRIP+1 
LZC -8LZC 
BMIR = BBMIR 
SUZC = BSUZC 
ETLF - BETLF 
SUZC -BBUZC 
SI AC = BSIAC 
KFFC = 1 
GO TO 142 

2 26 I F ( NRHP .FQ. C) GO TO 227 

; CORRFCT SYNTHESIZED RUNOFF TO RECORDED VOLUMES 

CALL ADJHYDl I CYB , IDYF, J HRB, 1 HRE, KP5H, DPY ,HBF ,NRHP, DSSF,DRSF»S SR , 

1 LSHA) 

; ESTABLISH STORM AND OVERALL OPTIMUM VALUES FOR SRX AND NCTRI 

CALL SETHPPCCTRI ,BTRI t WCFS ,CONOPT( 2 ) , HBF ,LSHA, S5R,NHPT ,KPSH, 

1 I B7PR, SRX, C SRX, FSRX»CHCAP, NRHP, RHPF, NCTRI, NBTRI ) 

I Ft NCTRI »EQ. 0) GO TO 228 

227 IFtTRIP • GE* NLTR) GO TO 228 
TRIP - TRIP +1 

GO TO 142 

228 CONTINUE 

IFtNSYC .LT. NSYT ) GO TO ICO 
STOP 

end 

SUBROUTINE AD JHYDtIDYB, I DYE, IHRB, IHRE,KPSH» DP Y, HBF, NRHP, DSSF, 

1 DRSF,SSR, LSFA) 


223 


48 


224 


225 


ADJUSTS SYNTHESIZED FLOW VOLUME TO MATCH RECORDED VOLUME FOR SETTING 
HYDROGRAPH ROUTTING PARAMETERS 
DIMENSION IOVBf 5) ,1 DYE (5 ) ,1HR8 ( 5 ) , KPSH(5 ) , S$Rf 5, 170 ) , 

1 DSSF (366) , DRSF (366 ) »HBF ( 5 ) »LSHA (5 ) , IHRE( 5) 


SHM08640 
SHM0865L 
SHMQ866L- 
SHM0867C 
SH MO 86 80 
SHMC869C 
SHMU87C C 
SH m 87 1C 
SHM0872D 
SHMQ8730 
SHM08740 
SHM0875C 
SHMG876C 
SH m 8770 
SHMC8780 
S HMD 87 90 
SH MO 8 800 
AD JU S S H MU 8 8 10 
SMM08820 
SHM08830 
SHH08840 
SH MO 8 8 50 
SHM08860 
SHM08870 
SHM08880 
SH MO 88 90 
SHM08900 
SHM08910 
SHM08920 
SHMQ893C 
SHM08940 
S HMD 8950 
SHM08960 
SHM08970 
SHMf'8980 
SHM08990 
SHM09000 
SHMU901C 
SHNC9&20 
SHM0903C 
SH MO 9040 
SHM09030 
SHM09Q60 
SHM09Q76 
SHM09080 
SHMd9090 
SHMfV9l©f» 
SHM09110 
SHMC9120 
SHM0913M) 
'SHH© 9X40 
SHMOnSO 
SH MO 9 160 
SH : M«91tB 


LOGICAL LSHA,LSHP 
INTEGER DAY,CPY 
LSHP = .FALSE. 

KRHP * I 
DAY * 274 
100 CONTINUE 

IF { LSHP) GO TC 102 

IF( IDYBt KRHP ) *NE. DAY) GO TO 107 

1 01 HTH« IHRB(KRHP) 

HBFM» 1.0 

IF ( DSSF( DAY) .LT. HBF(KRHP) ) HBFM= 0.0 

TSHV- (24.0-HTH) *( DSSFt DAY ) -HBF(KRHP) l*HBFM/24*0 

HBFM * 1.0 

1 F C DRSF( DAY) .LT.HBF (KRHP) ) HBFM= 0.0 

TRHV « (24.0 — HTH)*{ DRSFfOAY ) -HBF(KRHP) )*H8FM/24.0 

IF( IDYF(KRHP) .EQ« DAY) GO TO 104 

LSHP » .TRUE. 

GO TO 107 ■ ' t • r-‘ .. 

102 I F( DSSF(DAY) .LT* HBF(KRHP)) GO TO 103 
TSHV -TSHV * CSSFCDAY) •HBFC KRHP ) 

103 I F{ DRSF( DAY) .LT. HBF(KRHP) ) GO TO 104 
TRHV * TRHV + DRSF(DAY) -HBF (KRHP ) 

104 CONTINUE 


SH ML 91 BO 

SH Mu 9 190 
SHMC92C & 

V' 

SHM09220 

SH MU 9 2 30 

SHML924S 

SHML92S0 

SH MO 9 2 60 

SHM09270 

SHM09280 

SH MO 9 290 

SHM- 9>-' ■ 

SHMC931& 

SHM09320 

SH ML 933 L 

SHMG9340 

SHM0935C 

SHM09360 

SHM09370 

SHM09380 

SHM09390 

SHM09460 


IFUDYE(KRHP) .ME. DAY) GO TO 107 
HTH - IHRg(KRFP) 

TSHV -TSHV ~{ 24. O-HTH) * ( OSSF (DAY ) 
TRHV * TRHV- ( 24.0 -HTH) *(DRSF(DAY) 
LSHP “.FALSE. 

SHM“ TRHV/TSHsV 
LSHA(KRHP) ■ .TRUE. 

I F ( S HM t> GT® 8.0 *0R. SHM .LT. 0.12 
I F ( .NOT. LSHA(KRHP) ) GO TO 106 


•HBF ( KRHP ) ) /24.C 
-HBF(KRHP) )/24.0 


125) LSHA(KRHP) * • FAL! 


I F ( .NOT. LSHA(KRHP) ) GO TO 106 
KPCH =KP SH ( KRHP) 

DO 105 KHPT = If KPCH 

5 SSR( KRHP » KHPT ) “ SHM*SSR( KRHP, KHPT) 

6 WRITE! 6f 1 ) KRFPfSHM 

FORMAT (/ / 10X f ‘VOLUME ADJUSTMENT FACTOR FOR HYDROGRAPH* , 1 2, 
1 'EQUALS', Fl€. 4) 

KRHP = KRHP + 1 
I F ( KRHP . GT. KRHP) RETURN 

IF ( IDYBt KRHP) .£Q. I DYE ( KRHP— 1 ) ) GO TO 101 

07 CALL DAYNXT (DAY, DPY ) 

I F ( DAY oNG« 274) GO TO 100 
RETURN 

END 

SUBROUTINE DA YSUM ( DRSF , MEDCY »DPY , ATFV , TMTFWY ) 

SUMS DAILY VALUES TO GET MONTHLY AND ANNUAL TOTALS 

DIMENSION DRS F(366 ) , EMATF( 13 ), MEDCY ( 12 ),TMTFCY( 12) , TMTFWY ( 
INTEGER DAY, DPY 

SUM ANNUAL AND CLMULATI VE MONTHLY FLOWS 

SMATF(l) * O.C 

ATF = 0.0 

DO 101 DAY = 1,365 
ATF =ATF +DRSF(DAY) 




SH MO 94 10 

SHM09420 

SHM09430 

SHMC9440 

SHM094 50 

SHM09460 

SHM',;- 94 70 

SHM094B0 

S HMD 94 90 

SHM09S00 

SHM09510 

SHM09520 

5HH0953C* 

SHM^9540 

SHM09550 

SHM09560 

SH MO 93 70 

SHM- m 8( 

SHM!. «| fO 

SHM6960D ■ 

SHM09610 

SHM09620 

SHM09630 

SMM09640 

SHM0965L 

SHM09660 

SHMQ9670 

$HM09660 

SHM09690 

SHM09700 

SHM097 10 



DO 100 KMO = 2,12 

X F ( DAY oEQ» MEDC Y ( KMO ) ) EMATF(KMG) = ATF 

iXl CONTINUE 

EMATF< 13) =AT F 

ATFV =ATF +DR SF( 366 ) 

CALCULATE MONTHLY FLOWS 

DO 102 KMO = 1,12 

102 TMTFCY (KMO ) = EMATF(KMO + 1) -EMATF(KMO) 

TMTFCY (2) * T RTFCY ( 2 ) + DRSF (366) 

CONVERT MONTHLY FLOWS TO A WATER YEAR ORDER 
DO 103 KMO =1,9 

103 TMTF WY (KMO+3 ) =TMTFC Y( KMO ) 

DO 104 KMO = 10, 12 

104 TMTFWY (KMO-9 ) = TMTFCY ( KMO ) 

RETURN 

END 

SUBROUTINE FIXTRICCTRI »BTRI» NBTRI »NCTR1 ) 

FIX VALUES OF THE TIME ROUTING INCREMENTS TO MATCH REQUIRED TOTAL 
NUMBER OF VALUES 

DIMENSION AWSBIT (99) »BTRI (99 ) ,CTRI ( 99 ) 

IF ( NCTRI • GT* 99) GO TO 101 


SHMG972* ■ 
S BMC 977 C 
SH MO 9740 
S BMC 97 SO 
SHM09760 
SHMU 9770 
S BMC 97 80 
SHMC9790 
SHM09800 
SH MO 98 10 
SB MO 98 20 
SHM09830 
SHM09840 
SHM098S0 
SH MO 9860 
SHM09B70 
SHM09880 
SHM09890 
SH MO 9 900 
SHM09910 
SHM09920 


104 


99 USED* 


IF1NBTRI oNE. NCTRI ) GO TO 102 
DO 100 KRD * 1,99 

0 CTRI(KRD) = BTRI(KRD) 

RETURN 

1 WRITE! 6, 1 ) NCTRI 

1 FGRMAT(§X» "NCTRI OF* ,15, IX, ‘EXCEEDS MAXIMUM VALUE OF 99, 99 USEC* 
NCTRI * 99 

2 DO 103 KI A * 1,99 

3 AWSBIT (KI A) « 0,0 
FNTRI = NCTRI 

K81 = C • - ■ : • =- ,C: .';?*■ • LF 

KB2 = 1 

KB3 = 0 . • < 

4 KB1 =KB1 + 1 

I F ( KB1 o GT» NETRI) GO TO 107 

KS4 = 0 

WSBIT = BTRI(KBl) /FNTRI 
g K84 = KB 4 + 1 

I F ( KB4 cGT. NCTRI) GOTO 104 
AWSBIT (KB2 ) =AWSBIT(KB2) + WSBIT 

KB3 = KB3 +1 

IF ( KB3 .LT. NETRI) GO TO 106 

KB3 = l 

K82 = KB2 + 1 

6 GO TO 105 

7 DO 108 KBS = 1,9 9 

8 CTRI (KBS) = AWSBIT(KBS) 

RETURN 

CND 

SUBROUTINE PR ECHM DRGPM ,DRHP ,0R$ F, VWIN, SGRT ,NATRH) 

CHECKS PR! CIPITATION-STREAMFLOW ANOMALIES AND ADJUSTS PRECIPITATION 

~WHERG NECESSARY 

DIMENSION DRGPM ( 366 ), DRSF (366 ),DRHP( 366, 24) 


SHMQ9930 
SH MO 9 940 
SHM09950 
SHM09960. 
SHM49970 
)SHM09980 
SHM09990 
SHM 10000 
SHMlOOli) 
SHM1< 2 
' SHM10030 
SHM10040 
SHM IOC 50 
SHM 10060 
SHH 10 O 7 O 
SHM 16080 
SHM 10090 
SHM10100 
SHM 10 i 10 
SHM 10 120 
SHM 10130 
SHM 10 140 ■ 
SHM 1015 O 
SHMIC 160 
SHM 10170 
SHM 10180 
SHM 10190 
SHM l 1 .. 20! 
SHMIQ 210 
SHMlC 22 Ci 
SHMiO '230 
5 HM 1624 D 
SHM 10 2 56 



INTEGER DAY,HCUR,$GRT 
AHP = 0.0 
NR HA = 24--NATRH 
RGPM = DRGPM( SO) 

DAY=90 
RMWR«1*25 
50 DAY=DA Y+l 

I F{ DAY «GT«200 . OR. VWIN.GT. 750.0 ) RMWR=2.G0 
RFRlSE *{ DRSF (DAY) -DRSF{ DAY-1 )) /VWI N 

DO 101 HOUR = 1,24 

IF (HOUR • f:Q® SGRT+1) RGPM =DRGPM{ DAY ) 

AHP =AHP + DR t-P ( DAY, HOUR) *RG RM 
I F( HOUR • NE. fvRHA) GO TO 101 
RWRAIN =AHP 
AHP * 0.0 


SHM1A-26* f 
SHM li‘ 270 
SHM10280 
SH MIG 290 
SHM10300 
SHM 103 10 
SHM 1C 3 2* 
SHM 103 30 
SHM 10340 
SHM1035D 
SHM 103 60 
SHM10370 
SHM 10380 
SHM10390 
SHM10400 


101 CONTINUE 

IFtRFRISE .GT. RWRAIN .AND 
IF ({RWRAIN .GT.'RMWR .AND. 

1 .GT. 3.00 .AND. RFRJSE 

GO TO 108 

102 I F{ RWRAIN .GT. 0.051 GO TO 103 
RAA « RFRISE*2.0~RWRAIN +1.0 
DRHPtDAY, 13) = RAA 
WRITE(6,l> DAY, RAA 


RFRISE .GT. 0.1) GO TO 102 
RFRISE .LT. 0.02*RWRAIN) , 
.LT. 0.05*RWRAI N ) ) GO TO 104 


103 

104 
105 


FORMAT ( / lf)X» * FOR DAY*, 14, IX, 'RAIN ADDED 0F*,F7.2) 

GO 'TO 108 

RAM *2.0*RFRI SE/RWRA IN 
GO TO 105 

*R FRI SE/R WRAIN 
, 0.0) GO TO 108 
DAY, RAM, RWRAIN 

, ' FCR DAY* ,14, IX, 'RAIN ADJUSTMENT MULTIPLIER IS*,F8.4, 
F7.2) 


RAM = 10.0 
I F( RAM .LT 
WRITE (6, 2 ) 
FORMAT (/5X 
. IX 


, * RECORDED RAIN IS',! 

DO It 6 HOUR « 1, NRHA 

106 DRHPtDAY, HOUR ) =DRHP (DAY, HOUR) *RAM 
I F ( NATRH .£Q. 0) GO TO 108 

NFRHA *N R HA +1 
DO 107 HOUR =NFRHA,24 

107 DR HP { DAY- 1, HOUR) =DRHP ( DAY-1 , HOUR) *RAM 

; 108 IF (DAY .NE. 273) GO TO 100 

RETURN 

END 

SUBROUTINE RECES S( DRSF , DPY, BFRC, I FRC, ARE A, RSBD,RSBIF ,NRS , 
C ESTABLISHES RECESSION SEQUENCES 

DIMENSION DRS F(366) , LBFO(20) ,NDRS( 20 ) »R$BBF (20 ) ,RSBD ( 20 ) » 
1 RSBFRC ( 20 ) , RSB IF( 20) , RSI FRC ( 20 ), RSTFC 50,20 ) 

LOGICAL LBFO 

INTEGER DAY, DPY, RSBD,RSL 
REAL I FRC • 

MRSL * 9 

BFRC * 0.9 '■ 

I FRC * 0.05 

FRERS * 9. 1*5 CRT (AREA) 

DO 101 KSD ■ 1,50 ■ , ^CC;\CC, 


ICO 


SHM 104 IP 
SHM10420 
OR. (RWRAIN SHM1043O 

SH Ml 0440 

SHM10450 
SHM1046O 
SHM10470 
SHM10480 
SHM 10490 
SHM 10506 
SHM105I0 
SHM1O520 
SHM10530 
SHM 10 540 
SHM1O550 
SBM10560 
SHM 10 570 
SHM 10580 
SHM 105 90 
SHM 10600 
SMM1G610 
SH*$G620 
SHM 106 30 
SHM1064D 
SHM10650 
SHM1Q660 
SHM 106 70 
SHMLQ680 
SHM 10690 
sHMieiec 
SHMIDT 10 
SHM10720 
SHM10730 
SHM 10740 
SHM1015C 
SHM let 60 
SHM 16770 
SHM10780 
■ SHM1079P 


RSBBF) 



DO li.l KRS * 1,20 

101 RSTF ( KSD » KRS ) * 0.0 
KRS = 0 

DAY « 274 

BFGIN NSW SEQUENCE 

102 I F( KRS oGH* 20) GO TO 109 
KRS “KRS + l 

KSO * 1 

RSF1 * DR$F( DAY) 

CALL DAYMXT ( DAY, DPY ) 

IF (DAY .EQ. 274) GO TO 107 
RSF2 *DRSF(DAY) 

RSBD(KRS) = DAY 

IF ( RSF2 .LT. ( RSF 1+FRERS ) • AND. ( RSF2 

1 10.0)) GO TC 103 

KRS * KRS « 1 


SHM10800 
SHM 10 81© 
SHM10S2C 
SHM 10830 
SHM 10840 
SHM1Q65© 
SHM 10860 
SHM 10870 
SHM1G88© 
SHM 10890 
SHM 10900 
* SHM10910 

SHM 10920 

.GT. 0.4*AREA .OR. RSF2 . GT. SHM 10930 

SHM 1© 941 
SHM10950 


GO TO 102 

103 RSTF (1 ,KRS ) -8SF2 
RSFM * RSF2 
104 KSD « KSD + 1 

CALL DAYNXTC DAY, DPY ) 

I F( DAY .EQ. 274) GO TO 107 
RSFN =DRSF (DAY) 

IF( RSFN f LTo (RSFM +FRERS ) .AND. RSFN .GT. 

I F( K SD .GE. MRSL ) GO TO 102 
NOR'S (KRS > =0 
DO 105 KSD = 1, MRSL 
105 RSTF (KSD, KRS) =0.0 
KRS =KRS ~ 1 
GO TO It 2 

106 I F (RSFN .LT. RSFM) RSFM =RSFN 
RSTF (KSD, KRS) = RSFN 
NDRS (KRS ) =KSC 
I F { KSD . GE. 50 GO TO 102 

GO TO ID 4 

107 IF (KSD 6 GF. MRSL) GO TO 109 
NTRS =KRS - 1 
DO 108 KSD = 1,MRSL 

1 108 RSTF (KSD, KRS) = 0.0 
GO TO 11 a 
1C 9 NTRS = KRS 
lie CONTINUE 

IF { NTRS o GE. 3) GO TO 111 
I F ( MRSL . LT. 7) RETURN 
MRSL = 6 
GO TO 100 

WRI T c OUT ESTABLISHED ARRAY OF FLOW SEQUENCES 

111 WRIT! (6,1) 

1 FORMAT (/5X, 'FLOW SEQUENCES USED TO ESTIMATE 
DO 113 KRS * 1, NTRS 
NDRSC =NDRS( KRS) 

DO 112 KSD =2, NDRSC 
112 RSTF (KSD- 1, KRS) a R$TF( KSD, KRS) 

NDRS (KRS) *NDRS( KRS) - 1 


SHM10960 
SHM1097Q 
SHM1098© 
SHM10990 
SHM11000 
SHM 11010 
SHM lit* 2D 

0.0) GO TO 106 SHM1103G 

SHM 11040 
SHM11050 
SH*UO60 
sHMiioro 
SHMH080 
SHM11090 
SHMill©© 
SHM111T0 
sh M 11120 
SHM 1113© 
SHM1U40 
> SHM 11150 

SHMU1M 
SHM1H70 
SHMlllQfc; 

' -SHMIH9C 
3HM112Cv 

SBMU220 
'SHM I 1236 
SHM I 1240 ' 
SHM 11 250 
; SHM11260 
SH Ml 1270 

RECESSION CONSTANTS’) SHM l 128© 

'SHM1129© 
;SHM'113.CI0 
SHM 111 10 
: SHM 21 3 26 
SHM1133C 



NDRSC = NDRSC - 1 

WRlTfc. t 6, 2 ) KRS»{RSTF{K$D,KRS) , K$D=1 * NDRSC ) 

2 FORMAT t/lQX, I 2, 5 ( 10F8. 1 /12X ) ) 

L3 CONTINUE 

CETERMINL RECESSION CONSTANS FROM EACH SEQUENCE 
L 14 DO 116 KRS = 1,NTRS 

I F ( ( RSTF ( 1 »KR 5 ) .LT. 0.4*AREA) « AND® tRSTF(2»KRS) *GT* 0.8* 

1 { RSTFt 1 f KRS ))) ) GO TO 115 

LBFO(KRS) = .FALSE. 

CALL S£T2RC(RSTF, KRS, NDRS (KRS), RSI FRC (KRS), RSBFRC {KRS) ,LBF0(KRS 
I F t LBFOt KRS ) .OR. RSBFRCt KRS ) .GT. 1.2 .OR. RSBFRC (KRS ) «LT. 0. 
1 .OR. RSIFRC (KRS) .GT. 0.8 .OR. RSIFftC(KRS) . LT. (-0 .4) )G0 TO 1 

GO TO 116 

15 LBFO(KRS) = . TRUE* 

CALL SETIRCtR STF , KRS , NDRS ( KRS ), RSBFRCt KRS )) 

16 CONTINUE 

CALCULATE WEIGHTED AVERAGE RECESSION CONSTANTS 
BFRC * 0.0 


n 

6 

15 


17 

3 

18 


I FRC « 
ABFSL 
AIFSL 
DO 118 


0.0 

* 0.0 

» 0.0 
KRS = 


1.2 .OR. RSBFRC (KRS) .LT. 0.61 GO TO 117 


119 


1»NTRS 

I F ( RSBFRC t KRS ) . GT. 

RSL « NDRS t KRS) 

BFRC * BFRC + RSBFRC (KR$)*RSL 
AiFSl * ABFSL + RSL 
IFCLBFO(KRS) ) GO TO 118 
IFIRSL .GE. 2G.0 ) RSL = 20.0 
I FRC « I FRC + RSIFRC (KRS) *RSL 
AIFSL * AIFSL + RSL r 

GO TO 118' 

WRITE! 6*1 1 KRS 

FORMAT (IQX* * SEQUENCE '*13* IX* *OMITTED IN AVERAGING*) 

CONTINUE 

WRITE I 6« 4 1 ABFSL. AIFSL 

FORMATC10X, 'BASE FLOW DAYS =«, F5.0,2X, • INTERFLOW DAYS =*,F5.Q) 

BFRC = BFRC/ABFSL 

I FRC = IFRC/ AI FSL 

I Ft BFRC .GT. €.99) BFRC = 0.99 

IF (BFRC .LT. O.TO) BFRC » 0.70 

KSQ = 0 

DO 119 KRS = 1*NTRS 
I F ( LBFOt KRS) ) 60 TO 119 

CALL SC TRBFt RSTF* NDRS, KRS, BFRC, IFRC, CRSBIF,CRSBBF ) 

IFtCRSBIF .GT • 95000.0 .OR. CRS8BF .LT. O.C) GO TO 119 

IFtCRSBIF .LT. D.O) CRSBIF = 0.0 :•= V : ‘; 

KSQ = KSQ + 1 . “ 

RSBO (KSQ) = RSBOtKRS) . /I 

RSBIF(KSQ) = CRSBIF 
RSBBFt KSQ ) =CRSBBF 

CONTINUE ■■ f r 

NRS = KSQ 

DO 12» KSQ = 1* NRS 
DAY ■ RSBDtKS Q) 


IHfH; 


! *.&'■ 




SHM1134L 

SHM 11359 

SHM1136' 

SHM11371 

SHM11380 

SHM11390 

SHM 11400 

SHM11410 

SHM11420 

SHM11430 

SH Ml 144f 

SHM1145C 

SHM11460 

SHM11470 

SHM11480 

SHM 11490 

SHM115Q0 

SHM11510 

SHM 11520 

SHM1153C 

SHMU540 

SHM11550 

SHM11560 

SHM11570 

SHM 11 580 

SBM11S90 

SHMU600 

SHM1161G 

SHM11620 

SHMU630 

SHM l 1640 

SHM116SI 

SHM11660 

SHMU6T0 

SHM11&80 

$HMU*9C 

SHttUTOe 

SHMU710 
SHMU720 
SHMU730 
5HM11740 
SHM 1 ITS© 
SHMU760 
SHM 11770 
SHM11780 
SHM11790 
SHM 11800 
SHMU81G 
SHMUB20 
SHM 1 1830 
SHM 11840 
SHM11B50 
SHM11860 
SHM1187C 


CALL DAYNXT (DAY* DPY) 

L 20 RSBD(KSQ) * DAY 

DO 121 KSQ * 1»NRS 

CRSBTF * RSBIF(KSQ) + RSBBF(KSQ) 

121 WRITS(6»5 ) K$Q»RSBD( KSQ) ,RSB IF ( KSQ ) , RSBBF ! KSQ } , CRSBTF 
5 FORMAT </l0X» • REVISED FLOW SEQUENCE* » 13, IX, 'BEGINS ON DAY* ,14, 

1 * IX t * AT INTERFLOW “*,F7.2,1X, *CFS, BASE FLOW * * , F7 .2 , IX , *CFS, 

2 TOTAL = *»F?.2,1X,«CFS*) 

RETURN 

END - ’ r; :j . 

SUBROUTINE SETBI V{ BIVF»NRS,I FRC,RSBIF,SI FRS, FNCTRH ) 

SL : TS BEST VALUE CF BASIC INTERFLOW VOLUME FACTOR 
DIMENSION R5B IF( 20 ) , SI FRS ( 3, 20) 

REAL IFRC 


SHM11880 

SHM11890 
SHM 11900 
SHM11910 
SHM11920 
SHM1193L 
SHM11940 
SHM1195C* 
SHM11960 
SHM119TC 
SHM 1 1980 
SHM1199C 
SHM1200C 
SHM12010 


„ ;('i 


UP) )RSTR=3.€ 


im 

101 


ARSTR « 0,0 

DO 101 KRS “ I, NRS 
RIF » R56IF! KRS) /IFRC 
DO HtyKDY- 1,3 
RIF * RIF*IFRC 

* SIFRSCKCY, KRS) /IFRC ♦*{ FNCTRH/48.0) 

IFfRIF *GT • 0.0) RSTR = SIF/RIF 

IP! RSTR .GT. 3.0 .OR. C SI F . GT. 0.0 .AND. RIF .EG, 

ARSTR p ARSTR +RSTR 
WRITE 1 6,1 ) KR S»KDY*$ IF, RIF 

FORMAT ( 10X, *KRS= ' ,I3,2X, *KDY= * , I 2, 2X, * SI F= * , F7. 1 ,5X» *RIF=* , 

FT .1 ) 

CONTINUE . ■ - . -■ : - = . 

Continue 

TIRD * NRS*3 
PBIVF * 81 VF 
BIVF = 0.40 

I FC ARSTR .GT. 0.8) BIVF • {( Q8IVF-0. 40 )*TIRD ) /ARSTR + 0.40 
WRITE! 6, 2 ) PBIVF, BIVF 

2 FORMAT ( 5 X , * B I VF CHANGED FROM* , F6. 2, 2X, * TO* ,F6.2// ) 

RETURN 

END ■ ■ . ; ;■ V 

SUBROUTINE SETBMI { BMIR, NRS, BFRC, RSB8F, SBFRS, FNCTRH, I FT) 

C SETS BEST VALUE CF BASIC MAXIXMUM INFILTRATION RATE WITHIN WATERSHED 
DIMENSION RS36F(20)»SBFRS{3»20) 

ARSTR » 0.0 
DO 101 KRS = IFTjNRS 
R8F *5 RSBBF! KRS) /BFRC 
i DO 100 KDY « 1,3 
RBF = RBF*BFRC 


SHM12<’J2v 
SHM 12030 
SH Ml 2040 
SHM12050 
SHM 12060 
SHM12070 
SHM12080 
SHM12Q9D 
SHM12100 
SHM12110 
SHM12120 
SHM12131! 
SHM1214C 
SHM12150 
SHM12160 
SHM12170 
SHM12180 
SHM1219G 
SHM 12 200 
SHM1221C, 1 
SHM12220 
SHM 12230 
SHM12240 
SHM12250 
SHM 12260 
SHM12270 
SHM12280 
SHM l 2290 
SHM12300 


SBF * SBFRS! KEY, KRS) /BFRC**! FNCTRH/* 

RSTR « SBF/RBF ’ 

I F{ RSTR .GT. 3.0) RSTR *3.0 
ARSTR » ARSTR + RSTR 
WRITFJ6,!) KRS, KDY, SBF, RBF 
1 FORMAT (10X* * KRS= * , 13, 2X* *KDY S * , 1 2»2X» * 
1 FT m I ) 

100 COOT I HUE • • . ■ . 

101 CONTINUE 


" " ■ " ; ■ : ^ 

SHM1231D 


SHM12320 

0 ) 

hAsMyh v / . i . 1 

SHM 12330 

SHM12340 

C: ' ' " " • '• ' 

SHM12350 


SHM12360 


SHM12370 

SBF=*,F7.1,5X, 'RBF** , 

SH Ml 23 80 
SHM12390 
SHM12400 
SHM12410 


•RMFX/SMFX 


IBRD = (MRS + 1-IFT)*3 
ARSTR * ARSTR/TBRD 
ARSTR * ARSTR**1.3 
PBMIR »BMIR 
BMIR « PBMIR/ARSTR 
WRIT£(6»2) P8PIR»BMI R 

2 FORMAT (5X,*BM JR CHANGED FROM *, F6. 2 ,2X» * TO’ , F6. 2/ / ) 

RETURN < ' « 

END . vU-v ; . ■ •' . ' ■ 

SUBROUTj ME SETF0I ( MFDP*TMSTF , TMRTF ,SSQM) 

SETS VALUES OF FLOW DEVIATION INDICES 

DIMENSION MFDPU2) ,TMRTF( 12) ,TMSTF(12) 

REAL MFDP 
DO 101 MONTH = 1,12 
I F( MONTH .LE. 2) SSQM «= 0.0 
SMFX * TMSTF ( MONTH) + 20.0 
RMFX « TMRTFt MONTH) + 20.0 
MFDP (MONTH) * SMFX/RMFX - 1.0 
I Ft MFDP (MONTH ) «GT. 8.0) MFDRt MONTH) = 8.0 

ilFtMFDPt M onth j . lt. o.o) mfdp(month) = 1.0 

1 ’'/I'fPt MFDPT MONTH ) * LT. ( **’8.0) ) MFDP ( MONTH) * -8.0 

100 SSQM = SSQM + MF DP t MONTH) *MF DP ( MONTH ) 

101 CONTINUE 

WRTT'fet 6*1) (MFDP (MONTH), M0NTH=1,12), SSQM 
1 FORMAT (//2X, * MONTHLY DEVIATIONS* , /16X»12(F7.3» IX ) » *S$QM =*,F7.3) 

' * 'RETURN 

f * ■*<'« END 

SUBROUTINE SETFVPCLZC, SUZC,ETLF, BUZC, SIAC , TMSTF, TMRTF,TMPREC, 

1 TMP6T»£MLZS» SSQM,LRC»XMPFT» FTX»N0FM,L8UZC,LETLF,LLZC,APREC» APET 
; SETS BEST VALUES OF FLOW VOLUME PARAMETERS 

DIMENSION EMLZSt 12) , MFDPt 12) ,MXA ( 12 ) ,TMPET ( 12 ) ,TMPREC ( 12 ) , 

1 TMRTFt 12) , TMSTF ( 12 ) » XMPFT { 12 ) 

• LOGICAL LBUZC ,LETLF» LLZC,LRC 
REAL LZC , MFDP 

CALL SETFDI (MFDP, TMSTF fTMRTF , SSQM ) 

IF ( ( MFDP ( 2)+MFDP ( 3 ) ) .GT. 2.0 .AND. FTX .LT. 1.05) FTX = 0.9 
I F { ( MFDP ( 2 )+MFDP ( 3 ) ) .LT. (-2.0). AND. FTX .GT. 0.95) FTX=1.1 
C ADJUSTMENT OF LZC BASED ON MONTHS WHERE OVER HALF OF TOTAL 
C SYNTHESIZED RUNOFF IS OVERLAND FLOW, MINIMUM OF TWO MONTHS 
C WITH GREATEST RUNOFF USEC i* 

PL ZC = LZC 
FNOFM * NOFM 

' I F ( NOFM .GT. 2) GO TO 103 
MIR = 2 

M2R = l . 

IFtTMRTFt 2)*GT* TMRTF(l) ) GO TO 100 

, MIR = 1 bty " , 'H,- : , " . ■' ■. ■ i r - 

M2R « 2 

100 DO 102 MONTH =3,12 
IFtTMRTFt MONTH) .LT. 

IFtTMRTFt MONTH) .GT. 

M2R * MONTH 
GO TO 102 < ■ * 

101 M2R = MIR' 


gff:: 

yill 

...... 




r §1 ' $ 




TMRTF (M2R) ) GO TO 102 
TMRTFtMIR) ) GO TO 101 


SHM 12420 
SHM1243Q. 
SH Ml 2440 
SHM 12450 
SHM12460 
SHM12470 
SH Ml 2480 
SHM 12490 
SHM 12500 
SHM12510 
SHM1252D 
SHM 125 30 
SHM 12 540 
SHM1255D 
SHM 12560 
SHM1257U 
SHM12580 
SHM12590 
SH Ml 2690 
SHM1261G 
SHM1262U 
SHM12630 
SHM 12640 
SHM12650 
SHM12660 
SHM12670 
SHM12630 
SHM12690 
SHM12700 
SKM1271C 
SHM 127 20 
SHM12730 
SHM 1 274-; 
SHM 12*1 50 
SHM12760 
SHM12770 
SHM12780 
SH Ml 27 9(1 
SHM1280U 
SHM12B10 
SHM12820 
SKM12830 
SHM 12 840 
. SHM12850 
SHM12860 
SHM12870 
SHM12880 
SHM 12 890 
SHM12900 
SHM1291& 
SHM 1 2920 1 
SHM12930 
SHM12940 
SHM12950 



L 12 


PRM2 * TMPRSC (MONTH) 

CONTINUE 

FSUZC « MFDP(MSP) + MFDPt M2SP ) 
IF(ABS(XMPFT( ^ l*OI .GT* 0.2) GO TO 113 
FSUZC * FSUZC +MF0PU2) 


M12 


12 


'i ' & i ' ? * 
, ? ff * ’ 


113 IF( ABS CXMPFT t |1) - 1*0) .GT. 0.2) GO TO 114 
FSUZC - FSUZC MFDPtll) 

Mil • U , ' , 

1 14 IF( FSUZC .GT. 1.0) FSUZC = 1.0 

I F( FSUZC .LT. f-1.0)) FSUZC » -1.0 

IF{ FSUZC .GT. 0.0) SUZC » (FSUZC + 1.0>*SUZC 

IF ( FSUZC .LE. 0.0) SUZC ■ SUZC/ (1.0- FSUZC) 

WRITF;(6,3» SUZC,MlSP,M2SP t Mll,MI2 
3 FORMAT (6X» *SU«2C WAS CHANGED T0SF6.2,’ BASED ON MONTHS* ,413) 
IFCSUZC .IT. 0*3 .AND. LRC) SUZC * 0.3 
IF (SUZC .GT. 3.0 . AND.LRC) SUZC = 3.0 
ADJUSTMENT OF ft LF BASED ON SUMMER MONTHS OF RAINFALL EXCEEDING TWO 
INCHES OR NEED tO PREVENT MOISTURE BUILDUP 
I Ft EMC IS ( 12) .LT. PL ZC .OR. EMLZS(ll) .LI. PLZC .OR. APRFC .GT. 
GO TO 115 

MXAU) 13 

KWSM = 1 
GO TO ' 120 
115 SWSMD «* 0.0 
-'i.ltWSM «* 0 

DO 116 MONTH =1,12 
IF ( (MONTH .GT. MBWS .OR • MONTH < 

1 .AND. MONTH .LT. 9)) GO TO 116 
I F(TMPREC( MONTH) .LT. 2.0) GC TO 
SWSMD = SWSMD + MFDP (MONTH) 

KWSM = KWSM + 1 
MXA(KWSM) * MONTH 
CONTINUE 

I F ( KWSM .GE. 1) GO TO 117 

MXA(l) = MIR 
KWSM * 1 

FETLF » 5.0*MFDP (MIR ) 

GO TO 120 
WSM « KWSM 

I Ft .NOT. LFTLF .OR. KWSM .EG. 

EMFDP * 0.0 
DO 110 MONTH * 1»KWSM 
KM1 = MXA( MONTH) 

I F(MFDP( KM1 ) .LT. EMFDP) GO TO 118 
EMFDP = MFDP ( KM1 ) 

KM2 ■ MONTH I ‘ ■ 

CONTINUE • • : ' OOJ:, 

MXA( KM2) = ■ 0 
SWSMD * SWSMD -EMFDP 
WSM » WSM -i.C 
FETLF * 1.24SWSMD/WSM 
IFtFETLF .GT. 1.0) FETLF =1.0 


, GT. 2) .AND. (MONTH .LT. MBDS 


116 


116 


v 


i&n 


l) GO TO 119 


ir 4, 




t'A': 




IK 

fl! : v 

mmstm 




- * 

gipg ' ' • 





SHM 13500 

SHM1351C 

SHM1352C 

SHM13530 

SHM 13540 

SHM1355C 

SHM 13560 

SHM1357C 

SHM13586 

SHM13590 

SHM13600 

SHM 136 10 

SHM13620 

SHM13630 

SHM 13640 

SHM13650 

SHM13660 

SHM13670 

SHM1368D 

SHM13690 

SHM13700 

SHM 137 10 

SHM 137 20 

SHM13730 

SHM13740 

SHM1375& 

SHM13760 

SHM13776 

SHM 13780 

SHM13790 

SWM1380© 

SHM13813 

SMM13820 

SHM13830 

SHM13840 

SHM 13 8 5,0 

SHM13860 

SHM 13870 

SHM13880 

SHM13890 

SHM 13 90C 

SHM13910 

SHM 13920 

SHM13930 

SHM 13940 

SHM13950 

SHH13960 

SHM13970 

SHM 13980 

SHM13990 

SHM 14060 

SHM1401 S 

SHM 14020 

SHM 140 3G 


I F ( FETLF .LT • t-1.0 ) } FETLF s -1.0 
I F { FETLF .GT. 0.0) ETLF * { FETLF + 1.0)*ETLF 
IF(FSTLF .LT. 0,0) ETLF =ETLF/(1.0 -FETLF) 

WRJTFt 6,4) ETLF»«MfcMKM0J, KWO - i.KWSM) 

4 FORMAT (4X* »gTLF WAS CHANGED T0‘,F6.2t* BASED ON MONTHS 121 3) 

I F ( ETLF . LT* • Mi44tRC ) ETLF = 0 .05 

IF (ETLF *GT« €.6 .AND. LRCT ETLF = 0.6 
ADJUSTMENT OF BUZC BASED ON SEPTEMBER* NOVEMBER, AND DECEMBER 

kmi * :.z 

km 2 * 2 ;:ri; u , 

KM3 * 3 , . " . ' 

FBUZC • 3,4*t MFDPt 12 ) +MFDP( 2) + MFDP (3 ) ) 

1 Ft .NOT. LBUZC) GO TO 121 

FBUZC « 0.4* ( MFOP( 9) + MFDP(IO) + MFDP(ll) ) 

KMI » 9 

KM2 m W,'« . 

KM3 » U 

121 IF( FBUZC .GT. i.O) FBUZC * 1.0 

!E( FBUZC .IT* t-1.0) ) FBUZC = -i.O 
I Ff FBUZC .GT. 0.0) BUZC = (FBUZC + 1*0 ) *BUZC 

-sl‘f’4 0.0) BUZC=BUZC/( 1.0-FBUZC) 

) BUiZC,KMl,KM2»KM3 

FORMAT (4X,*BU>ZC WAS CHANGED TG*,F6.2, * BASED ON MONTHS*, 313) 
IF( BUZC '• LT. 0.2 .AND. LRC) BUZC = 0.2 
lElBUZC .GT. 4.0 .AND. LRC) BUZC = 4.0 
ADJUSTMENT OF SI AC BASED ON THREE FIRST MOISTURE EXCESS AND THREE 
FIRST MOISTURE DEFICIENT MONTHS 
KMI;* MBDS . ; 

KM2 - * MBDS + 1 

• -T'V- ; 


if » 


KM3 

KM4 

KM5 

KM6 

WFDX 


+ MFDP(MBWS+2))/3.0 
■ SIAO/O.4 




122 




MBDS 

©iltiM -h . 

© ; 3; s. I ;. % \ 

m , - rim ■ ' . 

o^Q 

IFtSIAC .GT. 1.0) GO TO 122 
WFDX ** (MFDP(MBWS) + MFDP( MBWS +1) 

IFtSIAC .GT. C.6) WFDX * WFDX*<1.0 

KM4 « MBWS ' ’ ' — . 1 

KM5 * MBWS + 1 
KM6 ■ MBWS +2 

SFDX = (MFDP(KBDS) + MFDP( MBDS+1 ) + MFDP { MBDS*”!) >/3 

FSIAC » 1.5*( SFDX-WFOX) 

I F ( F S I A C .GT. 1.0) FSIAC * 1*0 ' ■ - 4HV4 

I F ( F S I AC «L£. (-1.0) ) FSIAC * -1*0 
IFtSIAC .LT. 0.02) SIAC * 0.02 
IFtFSIAC .GT. 0.0) SIAC = (FSIAC + U0)*SIAC 
I Ft FSIAC .LE. 0.0) SIAC * SI AC/t 1.0-FSIAC) 

WRITE(6*6) SIAC, KH4 , KM5 1 KM6 , KM3 » KMI, KM2 
FORMAT t4X,» SI AC 
IFtSIAC .LT. C 

IFtSIAC .GT. 4.0 .AND. LRC) SIAC =? <r.v 

RETURN ■ >■? 

! HD ' ' 

SUBROUTINE SETHRPtCTRI »BTRI, WCFS»C0N0P2»HBF, LSHA,SSR,NHPT,KPSH 


?• *)t:l 

m 


m 




: . ^.4 ■ ■ ■ i/ :?*->! r ; '■ i 

■ : 

■ ■' ' 1 ' ' f 1 


WAS CHANGED TQ*,F6.2, 
.02 .AND. LRC) SIAC = 



,613) 


mm 


SHM14u4( 
SHM 14^50 
SHM14< 60 
SHM1407C 
SHM1406C 
SHM14* 90 
SHM1410G 
SHM14110 
SHM1412D 
SHM14130 
SHM14140 
SHM14150 
SHM 14160 
SHM141T0 
SHM14180 
SHM 14 190 
SHM1420U 
SHM 142 10 
SHM14221* 
SHM 14230 
SHM1424D 
SHM 14250 
SHM14260 
SHM 14270 
SHM 14 2c 
SHM14290 
SHM14SCI© 
SHM143 10 
SHM1432G 
SHM14S3',. 
SHM14B4C 
SHM143WO 
SHM 1436C 
SHM1437C 
SHM14380 
SHM14390 
SHM14400 
SHM14410 
SHM 144 2§ 
SHM 14430 
SHM 14440 
SHM14450 
‘ SHM, 14460 
SHM 1447 ft 
SHM14480 
SHM 14490 
SHM14SCH5 
SHM1451Q 
SHM14520 
SHM14530 
SHM 14546 
SHM14550 
SHM1456I 
SHM1457C 


1 IBTPRi S RX, C SRX » F SRX, CttCAP»NRHP»RHPF» NCTRI ,NBTRI ) 

•TS BEST VALUES OF HYOROGRAPH ROUTING PARAMETERS 

DIMENSION BTRIt99>fCTRI£99),HBF(5),HSRX(5),KPSH{5),LSHA(5>, 

1 HNTR I ( 5 I » RHm ® Lg | Wtt 5**7© ) » SSR 1 5 , 170 ) , T SRX £7) 

LOGICAL LSNA 

INTEGER CONG P2g'HNTR 1 1 SNTRJ ' 

REAL NHPT 

MHTP ** I ■ / O i,( ;MC' 

I F £ CONORt • »4 

I'IXTR! « 2*NBTRI 

MNTRI * ' 

TSRXU) » *095 

TSRXf 2 1 ; * S» 99 $ , 

TSRXC3) • ©•©•*'•>'*•• 

TSRXI4S ' 

T5RXI57 ' i 

'f 

DO M&mW” ltNRHP • 
g |rf ;# | ! W|' , -C&M;A<KHYDn GO TO 112 

fffTftt 1 ' ** '*} a .-t ^ v" w , ./ *’* r . •>r ‘ j 

ptktf CTRl»BTRI»N8TRI»NCTRI ) 

*»P‘| J ' 1 '*< •" V . 1 ^ ^ • ‘ 

■ , 

;; , '$$5 ' : ' 

SDRSP «’ 8 C. 0 *CHCAP 

SNTRI = MXTRI ; ■ ' : ... •, , , 

SRX a TSRX(KHl) ;; ; 

I F ( KH2 .EG., 21 LNIBRS =NIBRS 

WRITE(6, i ) MC TRI #SRX ,, , 

FORMAT £//15X» ‘TRIAL VALUE OF NCTR1 = *» I3» *» SRX * » F6#i) 

CALL TIMER! ( SSR *SRR»CTR I , NCTRI » KHYD, KPCH } 

CSRX * SRX 

CALL STORRTiS.BR. CSRX tFSRX»CHCAP»CONOP2t I BTPS »SHPF» KHYD.HBP(KHYD) 

1 NBPT* KPCH, IBTPRI 


1 11 i 

■HUp 1 

■HHp 

liiits 



LNTRX a ' NCTRI 
NIRTS « IBTPS- I BTPR*MHTP 
NlBRS m I ABS( KIRTS ) 

DRSP - ABSCSHPF - RHPF £ KHYD) ) 


IF(N1RTS .EG. 
1 RHPP(KHYO) 

I R IN I RTS ,GE. 
NCTRI * NCTRI 
IFl'NCTRI *LT* 
•IFCNCTM "*6T# 
CALL 
KH2 


."' 2 '' I 

TO 


0 .OR. CKH2 .EG. 

.GT* 1.2*SHPF) GO 
1) GO TO 109 
— NIRTS 

MNTRI ) NCTRI * MNTRI 

MXTRI) GO TO 106 

FIXTRUCTRI, BTRI»NBTRI, NCTRI > 
2 


AND* 

103 


NIBRS .GE. LNIBRS) .OR, 


Jfci, ‘ 






iFlDRSP.Gt. SDRSP) GO TO 108 

SNTRI. 1 * ‘ L'NTRI 
SDRSP'^'/ORSP' ■- 1 . ■ 1 1 


SHM 14 530 
SHM 145 9b . 
SHM146M. ! - 
SHH146.il. 
SHM1462',; 
SHM14630 
SHM1464G 
SHM1465C? 
SHM14660 
SHM14670 
SHM146 80 
SHM 14 £9', 
SHM147O0 
SHM14710 
SHM1472C 
SHM14730 
SHM14740 
SHM 14750 
SHM147 60 
SHM 14770 
SHM14780 
SHM 147 90 
SHM 148C*0 
SHM14810 
SHM14820 
SHM1493© 
SHMl4®'4© 
SHM14830 
SHM14S60, 
SMHI487G 
SHM14880 
SHM1489C) 
SHM 14 900 
SHM1491C 
SHM 14920 
SHM1493C 
t SHM 14940 
SHM14950 
SHM 14960 
SHM14970 
SHM 14980 
SHM 1.4996 
SHHL'5000 
SHM1501O 
SHM 150 2CI 
SHM 15630 
. SHM15046 
SHM15050 
SHM15060 
SHM 1567 0 
SHM1508O 
SHM 150 90 
SHM 1510b 
SHM 151 10 


I 

1 


KH3 = 2 

L kA KH1 = KH1 + 1 

IF(KH1 .FQ. 8) GO TO 105 
KH2 » 1 

GO TO 100 

L 05 HNTRI ( KHYD ) = LNTRI 
HSRX ( KHYD ) ■ SRX 

GO TO 111 

Wb IF(KH1 .GE. 2 .AND. KH3 «EQ. 2 .AND. DRSP .GE. SDRSP) GO TO 108 

NCTRI » MXTRI 

CALL FIXTRI ( CTRI , BTRI ,NBTRI » NCTRI ) 

IF(KH2 .EQ. 2 .AND. KH1 .EQ. 1 .AND. SHPF .GT. RHPF{ KHYD) ) GO 

1 TO HIT 

IF( KH2 . EQ. 2 .OR. KH1 .GE. 2) GO TO 109 

GO TO 1J2 

107 HNTRI ( KHYD ) = MXTRI 
HSRX (KHYD ) » 0.995 

GO TO 111 

108 HSRXCKHYD ) * TSRX(KHl-l) 

HNTRI (KHYD) * SNTRI 

GO TO 111 

109 IFtNCTRI .GT. MNTRI .AND. NCTRI .LT. MXTRI) GO TO 101 
IF (DRSP .GT. SDRSP) GO TO 110 


SHM15 1 20 

SHM1513© 

SHM1514& 

SHM15150 

SHM15160 

SHM15170 

SHM1518C 

SHM15190 

SHM1520C) 

SHM 1521© 

SHM 15220 

SHM15230 

SHM15240 

SHM15250 

SHM1526G 

SHM15270 

SHM15280 

SHM 15290 

SHM 15 300 

SHM15310 

SHM1532G 

SHM1533( 

SHM15340 


SORSP « DRSP 
SNTRI = LNTRI 
GO TO ID 4 

HNTRI (KHYD) = NCTRI 
HSRXCKHYD) * 0.995 

IFtKHl -GE. 2) HSRX ( KHYD) = TSRXCKH1 - 
I F { H SRX ( KHYD } .LT. 0.91 .AND. SHPF «LT< 
L * .FALSE. 

I F ( . NOT • LSHA (KHYD) ) GO TO 112 
WRITE! A, 2) KH YD, HNTR I ( KHYD ) , HSRX ( KHYD ) 


NCTRI =•, 13, 


FORM AT ( 17X , ' FCR STORM »,I2,' NCTRI =*,13 

CONTINUE 
KPA = 1 
ARHPF = 0.0 
APPKP = 0 . 0 

DO 114 KHYD = 1, NRHP 
I F ( .NOT. LSHA (KHYD) ) GO TO 114 
CHPV = HNTRI (KHYD) 

I F ( K PA .EQ. 2) CHPV = HSRX (KHYD ) 

APPKP = APPKP + CHPV *RHPF( KHYD) 

ARHPF = ARHPF + RHPF(KHYD) 

CONTINUE 

WAPV * APPKP/ ARHPF 
I F ( KPA . EQ. 2) GO TO 115 
NCTRI = WAPV + 0.5 
WRIT, ( 6 , 3 ) NCTRI 

FORMAT (/ / 10X , 'OPTIMUM NCTRI =',I3) 

I F ( NCTRI .EQ. 0) RETURN 

CALL FIXTRI (CTRI ,BTRI ,N8TRI , NCTRI ) 

WRITM' ,4) (CTRI (KTRI), KTRI = 1, NCTRI) 

FORMAT (18X, • CTRI ARE'/9(16X, 11F8.4/) ) 


HPV = HSRX (KHYD ) 
CHPV *RHPF( KHYD ) 
RHPF(KHYD) 


SHM15350 
SHM15360 
SHM153T0 
SHM15380 
SHM15390 
SHM 15400 

.5*RHPF(KHYD) ) LSHA( KHYD) SHM15410 

SHM1542SI 
SHM 15430 
SHM 154 4 

* SRX = * »F6.3 ) SHM15450 

SHM l 3460 
SHM 1547© 
SHM1548© 
SHM1549C 
SHM 1550© 
SHM15510 
SHM15520 
SHM1553G 
,, SHMlfMO 






iX, 'OPTIMUM NCTRI =* ,13) 





SHMl 
SHM15&&0 . 
SHMl 5640 
SHM15650 


THE USER MAY HAVE 
TO COMPENSATE FOR 


TO ADJUST THESE 
ROUNDING. * ) 


VALUES TC 


WR ITL ( 6, 5 ) 

5 FORMAT U.3X,' WARNING- 
I MAKE THEM ADD TO ONE 
KPA = 2 
GO TO 113 
115 SRX = WAP V 
CSRX = SRX 
FSRX = SRX 

CALL S K TSRP(CCNQP2,NRHP,LSHA,RHPF,HSRX,CHCAP,SSR,SRR,CTRI ,CSRX, 
1 FSRX, KH YD, I ETPS, $HPF,NCTRI ,HBF , NHPT , KPSH , I 8TPR 5 

SRX * CSRX 
RETURN 

END 

SUBROUTINE Sl : TRBFt RSTF *NDR$, KRS , BFRC, I FRC, CRSBI F »CRS BBF ) 

SETS VALUES OF INTERFLOW AND BASE FLOW AT RECESSION BEGINNING 
DIMENSION RSTF(5t)» 20) »NDRS(20) 

REALMS RA1»RA2,RA3,RA4,RA5,RA6 

REAL I FRC 


RA1 * 
RA2 « 
RA3 * 
RA4 » 
RA5 * 
MNDRS 


0.0 
0« 0 

0«. ? 

0.0 
0. 0 

* 12 


«LT © 
C.3) 


12) MNDRS = 
GO TO 101 


NDRS(KRS) 


100 


101 


I F( NDR$( KRS ) 

I Ft I FRC .GE. 

CRSBIF = 0.0 

DO 100 KSD = 1, MNDRS 

RA1 = RA1 + BPRC**(2*KSD) 

RA4 = RA4 + RSTF (KSD,KRS)*(BFRC**KSD) 
CRSO 8F = RA4/RAI 
CRSBBF = CRSO EF* BFRC 
RETURN 

CRSBIF = 1000 CO. 0 

DO If, 2 KSD = 1, MNDRS 
RA1 = PA1 + BFRC**( 2*KSD) 

RA2 = RAP + I FRC**(2*KSD) 

RA3 = RA3 + { 6FRC*IFRC)**KSD 

R.A4 = RA4 + R$TF(KSD,KRS)*(BFRC**KSD) 

RA5 = PA 5 + R STF ( KSD »KRS ) * ( I FRC**KSO) 


102 


CONT INUE 
RA6 = RAW A 2 
I F ( R A6 ,F«. fi 


-RA3**2 
0) RETURN 
= - ( R A3/R A6 ) #RA4 
* CRSO IF* I FRC 
= (RA2/RA6) *RA4 
= CRSO EF*BFRC 


SHM15660 
SHM 15670 
5HM15680 
SHM15690 
SHM 15700 
SHM 15? K 
SHM 15? 20 
SHM 15730 
SHM15740 
SHM15750 


+ ( RA1/RA6 )*RA5 


(RA3 /RA6 ) *RA5 




■C 


CRSO IF 
CRSBIF 
CRSOBF 
CRSBBF 
RETURN 
END 

SUBROUTINE SETSRP { C0N0P2 »NRHP» LSHA»RHPF » HSRX »CHCAP »SSR»SRR,CTRI t 
1 CSRX t FSRX, KFYDflBTPS, SHPF, NCTRI ,HBF, NHPT, KPSH, IBTPR) 

SETS BEST VALUES OF STORAGE ROUTING PARAMETERS 

DIMENSION CTRK99) ,HBF ( 5 ) ,HSRX { 5) , KPSH( 5 ) , LSHA (5 ) ,RHPF{ 5 ) , 

1 SRR ( 5, 170) , SSR (5,170) 


SHM15760 
SHM15770 
■ SHM15780 

SHM15790 
SHM15800 
SHM15810 
SHM 158 20 
SHM15830 
SHM15S40 
SHM15850 
SHM15860 
SHM15870 
SHM15880 
SHM15890 
SHM 15 900 
SHM15910 
SHM15920 
SHM15930 
SHM15940 
SHM1595C 
SHM 159*0 
SHM15970 
SHM1593CI 
SHM15990 
SHK16&00 
SHM 16,010 

mmmw 

smimm 

SHM160S© 

sHMimm 

SWMlkiTt) 

mntMi'm 

mmmm 

' SHM161F,..,- 
SHMltllv 
SHM161 2' 
SHM161V 
SHM1614' 
SHM 16 152 
SHM161&.- 
SHKU17& 
SHK16180 
SHM1619D 


LOGICAL LSHA 
INTEGER CONOP2 
REAL NHPT 
KLCCA = i 
SRX = CSRX 
EPS « O.OOOOOl 
NORHP * NRHP 


SHM 162N0 
SHM1621* ■ 
SHM16220 
SHM16230 
SHM 16 240 
SHM16250 
SHH162.60 


NRHP a NRHP-1 


2) GO TO 
LINE FOR 


103 

DETERMINING 


HSRX (KHYD ) 

RHPF (KHYD)*H$RX{ KHYD) 
RHPF (KHYD ) **2 


DO 100 KH YD « 1» NORHP 
IF( .NOT. LSHA(KHYD) ) 

100 CONTINUE 
I FC NRHP .LE. 

FIND REGRESSION 
RAI = U.O 
RA2 =0.0 
RA3 = 0.0 
RA4 =0.0 
FNRHP = NRHP 
00 101 KH YD * 1, NORHP 
IF (.NOT. LSHA (KH YD) ) GO TO 101 
RA1 = RAi + RFPF(KHYD) 

RA2 « RA2 + 

RA3 * RA3 + 

RA4 » RA4 + 

101 CONTINUE • 

AVRHPF » RAI/ FNRHP 
ASRX » RA2/FNRHP 

RSLP = ( RA3 -RA1*ASRX)/(RA4-RA1**2/FNRHP) 
IFfRSLP .LE. EPS) GO TO 106 
RINT = ASRS -RSLP*AVRHPF 

102 CSRX = RINT + R SLP* (0. 5*CHC AP ) 

I F { C SRX « GE. 0.99) RETURN 
IFtCSRX .LE. C.S) CSRX = 0.8 

GO TO i; 7 

103 K1AH =2 

DO 104 KH YD = 1, NORHP 
IF (.NOT. LSHA (KHYD) ) GO TO 104 
I F( K 1AH .EO. C) K1AH = KHYD 
I F ( K1AH .GTo C) K2AH= KHYD 

10 4 CONTINUE 

IF (NRHP .EQ. 

C FIT THE STRAGHT 

RSLP = ( HSRX ( K1AH 
I F ( RSLP .LE. EPS) 

R I NT = HSRX(KIAH) 

GO TO 122 
105 CONTINUE 

CSRX = HSRX(KIAH) 

FSRX = CSRX 
GO TO 115 
1C6 CONTINUE 
: CSRX a sax 

FSRX = CSPX 
! WI> IT - (otl ) 

t i FORMAT (//10X» 


CSRX, WHEN NRHP EXCEEDS 3 


1) GO TO 105 
LINE WHEN NRHP = 2 

-HSRX(K2AH) ) / ( RHPF{ K1AH )■ 
GO TO 106 
-RSLP*RHPF (K1AH ) 


RHPF( K2AH) ) 


SHM16270 
SHM 16280 
SHM16290 

SHM 16300 
SHM1631S 
SHM16320 
SHM16330 

SHM16340 

SHM16350 

SHM16360 

SHM16370 

SHM16380 

SHM16390 
SHM 16400 
SHM16410 
SHM16420 
SHM 16430 
SHM 16440 
SHM16450 
SHM16460 
SHM16470 
SHM16480 
SHM16490 
SHM16S00 
SHM 165 1-0 
SHM16520 
SHM16530 
SHM 16540 





•REGRESSION LINE HAS NEGATIVE SLOPE*) 


SHML6570 

SHM1656L 
SHM16590 
SHM16680 
SHM16610 
SHM 166 20 
SHM 16630 
SHM1664C 
SHM 16650 
SHM 16660 
SHM 16670 
SHM16680 
SHM1669* 
SHM 16700 
SHM167U 
SHM16720 
SHM 167 30 


GO TO 115 
L 07 CONTINUE 

BISRX = f).2*{Ci.99 - CSRX ) 

SISRX * 0.04* (0.99-CSRX) 

TFSRX = CSRX 
KISRX = 0 

1*8 KISRX » KISRX + 1 
FSRX * TFSRX 

WRITE ( 6» 2 ) KISRXtCSRX, FSRX, CHCAP 

2 FORMAT I//15X, *TRIAL* , 13, ’» CSRX *• , F8.5, ' » FSRX =»,F8.5, 
1 S CHCAP = *,F1O.0) 

SQPKO = 0.0 

ADRSP * 0*0 

00 1C 9 KHYD a If NORHP 

IF (.NOT. LSHA (KHYD) ) GO TO 109 

KPCH * KPSH(KHYD) 


SHM1674C 
SHM16750 
SHM 16760 
SHM 16770 
SHM16730 
SHM 167 90 
SHM16800 
SHM16310 
SHM16820 
SHM16S30 
SHM16640 
SHM 16850 
SHM16860 
SHM16870 
SHM16880 
SHM16890 


CALL TIMERT( SSRt SRR,CTRI t NCTRI , KHYD, KPCH ) 

CALL STDRRTf SRR, CSRX, FSRX, CHCAP »C0N0P2, I BTPS,SHPF, 

1 NHPTfKPCHfl BTPR) 

DRSP a SHPF » RHPF(KHYD) 

SQPKO * SQPKO + DRSP **2 
ADRSP = ADRSP +DRSP 

CONTINUE 

WR Iff (6,3) SQPKD 

FORMAT I/25X, * SQPKD =*, F14.0) 

I F{K ISRX .NE. 1) GO TO 110 
TFSRX * CSRX + B ISRX 
SSQPKD = SQPKC 
BFSRX = FSRX 
GO TO 1C* 8 

) I F ( SQPKD .GT. SSQPKD) GO TO 113 

I F ( K ISRX o FLU 6 o ANDe ADRSP .GT« 0.0) GO TO 111 
SSQPKD = SQPKC 
BFSRX = FSRX 

I F C KISRX . GF. 11) GO TO 114 
I Ft KISRX oLE. 5) TFSRX * TFSRX +BISRX 
I F( KISRX .GF. 6) TFSRX =TFSRX -SISRX 

GO TO 108 

KLCCA « KLCCA + 1 
I F(KLCCa .GF. 5) GO TO 112 


SHM1690G 

KHYD, HB F ( KHYD ) , SHM16910 

SHM16920 
SHM 16930 
SHM 16940 
SHM16950 
SHM 16960 
SHM16970 
SHM16980 
SHM1699D 
SHM 17000 
SHM 17010 
SHM17020 
SHM 17030 
SHM 17 4* 
SHM17050 
SHM17«60 
SHNlfOtO 
SHM17D8C 


CHCAP = 0.8*CKAP 

CSRX = RTNT + I 
GO TO IT 7 
CSRX = 5.990 
FSRX = 0.990 
GO TO 115 
I F ( KI SRX oGT. 
KISRX = 6 
SSQPKD = SQPKC 
BFSRX = FSRX 
TFSRX = TFSRX 
GO TO 138 
FSRX * BFSRX 


= RINT - 
1 7 

* 3 o 990 
= 0.990 
115 

oGT 


RSLP*(0.5*CHCAP) 


SHMltlfH 
17140 
SHM17340.. 
SHHltiSO 
SHMltlOO 
‘SHMilTIf C3 
SHM lft SO 
SH'M'If-tDO 


SHM 172 10 


•SISRX 


CSRX, FSRX, SSQPKD 


SHM1T24© 
SHM17250 
SHM 17260 
SHM17270 


*' FORMAT (// 10X, *CSRX =',F7.4,1GX, *FSRX ■* , FT, 4 , 10X , * SQPKD *»,F15. 

RETURN 

L 15 WRITE (6, 5) CSRX,FSRX 

5 FORMAT (//10X, *CSRX » * , FT. 4, 10X , 1 FSRX =«,F7.4) 

RETURN 

KND 

SUBROUTINE SGTIRCt RSTF, KRS, NDRSC, BFRC ) 

SETS BEST VALUE FOR ONE RECESSION CONSTANT 

DIMENSION RSTF (50* 20 } 

RA1 *0.0 
RA2 = f).3 

NDRSC1 * NOR SC -I 

DO 100 KSD = 1,NDRSCI 
RAI = RAI + RSTF(KSD,KRS)**2 
100 RA2 * RA2 + R STF {KSD, KRS ) *RSTF ( KSD+1 , KRS ) 

BFRC * RA2/RA3 
WRITEU,!) KRS, BFRC 

1 FORMAT USX, * KRS ** , 1 3, 5X, * BFRC =*,F8.4) 

RETURN 

END ■ ‘ 

SUBROUTINE SET2RC( RSTF , KRS,NDRSC, I FRC , BFRC, LBFO) 

SETS BEST VALUES FOR TWO RECESSION CONSTANTS 


2) SHM1728D 
SHM17290 
SHM 17300 
SHM17310 
SHM1732U 
SHM17330 
SHM17340 
SHM 17 350 
SHM1736C 
SHM17370 
SHM17330 
SHM17390 
SHM 17400 
SHM17410 
SHM 17420 
SHM17430 
SHM 17440 
SHM17450 
SHM17460 
SHM 17470 
SHM 17480 
SHM1749© 


01 MENS ION RSTF(50,20) 

LOGICAL LBFO 
1t£4£ I FRC 

REAL*8 RAI, R A2, RA3, RA4,RA5, CRSTFC 50 ) ,RA6, DBFRC, 
DO 100 KSD * 1, NDRSC 
ICO CRSTF(KSD) » RSTF { KSD, KRS ) 

NDRSC2 = NDRSC -2 

RAI * 0 

RA2 = ( .0 
RA3 = 0.0 

DO ID! KSD = 1,NDRSC2 


DIFRC,RA,RB,RD 


RAI + C RSTF (KSD) **2 

RA2 +CR STF ( KSD )*CRSTF( KSD+1) 

RA3+CRSTFI KSD) *CRSTF (K SD+2 ) 

RAI +CRSTF ( NDRSC- 1)**2 -CRSTF ( 1 ) **2 
RA2 + CRST F (NDRSC- 1 ) *C RSTF (NDRSC ) - CRSTF( 
RA4*RA1 “RA2**2 
■ « EQ. 0.0) GO TO 102 


RAI = RAI + C 
RA2 = RA2 +CR 
RA3 = RA3+CRS 
RA4 * RAI +CR 

RAS * RA2 + C 
RA6 = RA>*RAi 
I F (RA6 .EQ. 0 
P.A5 = RA5/RA6 
RA3 = RA3/RA6 
RA = R Al *RAB 
RB * R A4*RA3 
kD = RA**2 + 
IF (P.D oLT © 0. 


SHM1750G 

SHM17510 

SHM17520 

SHM17530 
SHM 1754V 
SHM17530 
SHM17560 
SHM 17570 
SHM1758C 
SHM1759G 
SHM1760C 
SH«i7&lC 


I)*CRSTF<2) 


SHI 

. i , 


Sill 


SHM1T640 


DBFRC 
DIFRC 
BFRC * 
I FRC = 
WRIT! 
FORMA! 
GO TO 
LB Ft) « 


- RA2*R A3 

- RA2*RA5 
4. 0*R8 

C) GO TO 102 
RD **Q. 5 ) /2.0 


= (RA +RD **0.5)/2 

• RA "DEFRC 
DBFRC 
DIFRC 

>,l) KRS, BFRC, IFRC 


,'KRS *• ,I3,5X, 'BFRC *• , F8®4,5X, * IFRC 


TRUE 


swMWH 

SHM17690 

mninm 
shmoti© ^ 

1 1 SHM 17720 1 
■ ' mnntm, 

SHMltf 40 : 

, SHM 17750 


wv 

. . ■ 


SHW1778©' 

*• ,F8*4) SHKttWO 

'SH Ml 7800 
SHM17810 


2 

Li* 3 

PE 


Wft I7F ( o» 2 ) KRS 

FORMAT { / 1 5 X f * IMAGINARY VALUES ENCOUNTERED IN SET 2RC» SEQUENCE 

1 13) 

RETURN 

END 

SUBROUTINE ST CRRTC $RR,CSRX, F SRX, CHCAP, C0N0P2 , I BT PS,SHPF, KHYD, 

1 CHBFfNHPTt KPCHtlBTPR) 

RFORMS CHANNEL STORAGE ROUTING 
DIMENSION A$RR(5,21) ,$RR(5,170) 

INTEGER CONOP 2»PRD 
REAL NHPT 

WR I T E ( 6 » i ) CHEF 

t 

CHBF 

NHPT 


me 


i*; f i 


102 


103 

1C4 

2 


FORMAT (/25X, * BASE FLOW «* ,F7.1, * CFS* J 


TFCFS 
INHPT 
MHTP » I 

I F ( C0N0P2 »EQ« 0) MHTP 
-> MHTP* INHPT 
0.0 : : 

0.9*SRR(KHYD,1) 

'0 V— ■' ' ■■ 

KHPT * li KPCH 


INHPT 

SHPF * 

RHF0 * 

KAFH » 

DO 102 
PRD 

PRD * PRD + 1 

mm T* »'$'RRCKHVDf KHPT) 

I F{ TFCFS .LE. 0 . 5 *CHCAP) SRX = CSRX 

I F € (TFCFS .GT. D. 5 *CHCAP) .AND. (TFCFS .LT. 2 . 0 *CHCAP) ) SRX = 
1 CSRX + ( FSRX -CSRX )*( (TFCFS - 0 . 5 *CHCAP )/( U 5 *CHCAP ) )**3 
I F( TFCFS .GE. 2 . 0 *CHCAP) SRX - FSRX 
RHF 1 * TRHF- SRX*(TRHF «RHFO) 

RHFO * RHF 1 

TFCFS * RHF1 4 CHBF 

I F( TFCFS .LT. SHPF) GO TO 101 

SHPF = TFCFS 

1 BTPS = KHPT - ‘ ‘ 

I F ( PRO oLEo 3 .AND. CONOP2 . EQ. 1) GO TO 100 

KAHP * KHPT - IBTPR*MHT P 4 - 5 *INHPT 

I F { KAHP .LT. C) GO TO 102 

IF(MOD(KAHP t INHPT) .NE. 0 ) GO TO 102 

KAFH * KAFH 4- 1 

ASRR ( KHYD t KAFh) * TFCFS 

CONTINUE 

IF (KAFH • FQo 21) GO TO 104 
KAFH = KAFH + 1 
DO 1U3 KIA » KAFH, 21 
ASRR (KHYDfKIA) » O.D 

WRITE (6, 2 ) KHYD, NHPT, (ASRR(KHYD,KWD) 


KWD * 1,21) 


J ; 




p 

• 1 


FORMAT (/25X, ' SYNTHESIZED HYDROGRAPH' , 13, • INTERVAL «*,R5.2, 

1 'HOURS' /*»( 2 2 X,TF 10 . 1 /) ) 

WR I T 11 ( 6t *4 ) SHPF 

FORMAT ( 25 X, 'FLOOD PEAK *' ,FlC.i, • CFS* ) 

RETURN . / ' v », 

SUBROUTINE STRHR S( RHPD , RHPH, IDYB,I DYE, IHRB, I HRE , NHPT , MXTR H, DP Y , 


SHM17820 

SHM 1783 D 

SHM17340 

SHM 17 850 

SHM17860 

SHM1787Q 

SHM17880 

SHM17890 

SHM17900 

SHM17910 

SHM 17920 

SH Ml 7936 

SHM1794D 

SHM17950 

SHM1796G 

SHM17970 

SHM17980 

SHM1799D 

SHM18OO0 

SHM18010 

SHM 18020 

SHM1803O 

SHM l 8040 

SHM18O50 

SHM18080 

SHMljMWt) 

SHM1808D 

SHM 18090 

SHM181O0 

SHM 18 110 

SBM181 20 

SHM 18 147: 
SHM 181 SO 
smiBim 
SHM 18 170 
SHMieiec 
SHM18190 
SHM182W 
SHM13210 
SHM 18 220 
$HM18330 
SHM1824D 
SHMISZSD 
SHM 18200 
$HM18270 
SHM18280 
SHK18290 
SHM 18300 
'SHM 18 31© 
SHM18320 
SHM 18330 
SHM Ife 340 
SHM18350 


1 NR HP, I BT PR } 

S;:TS BEGINNING AND END TIMES OF RUNOFF ENTERING RECORFD HYDROGRAPHS 
DIMENSION RHPC{5),RHPH(5),IDYBC5) , I DYE ( 5 ) , IHR8 (5 ),IHRE(5 ) 

INTEGER DAY,DPY,RHPD,RHPH 

REAL NHPT 

ESTIMATE HOURS EACH WAY FROM PEAK 
I NHPT = NHPT 
IBTPR = 5 * I N H PT + MXTRH 
I PTE * 15*INHPT 

DETERMINE TIME OF BEGINNING AND ENDING FOR EACH STORM 
DO 106 KRHP * IfNRHP 
KHBCK » IBTPR » RHPH(KRHP) 

IFCKHBCK . LT. 0) GO TO 101 
KDBCK = KHBCK/24 + 1 
IHRBCKRHP ) * 24*KDBCK- KHBCK 
DAY » RHP DC KRHP) 


SHM 18360 
SHM1B37C 
SHM 18460 
SHM 18390 
SHM18400 
SH Ml 94 10 
SHM 1 8420 
SHM 18430 
SHM18440 
SHM18450 
SHM 18460 
SHM18470 
SHM18480 
SHM18490 
SHM18500 
SHM18510 


DAY * DAY -1 

IFCDAY.EQ. 59 .AND, DPY .£Q< 

I F( DAY «6Q. 365) DAY = 59 
IFCDAY .BQ* 0) DAY * 365 
KDBCK * KDBCK - 1 
IR KDBCK .GT. 0) GO TO 100 
IDYB(KRHP) =* CAY 
GO. 10 132 
IDYB (KRHP) * RHP DC KRHP ) 
1HRBCKRHP) * RHP H (KRHP )- IBTPR 
KHFOR = IPTE + RHPHC KRHP) 
IFCKHFOR .LE. 24) GO TO 105 


366 ) DAY = 


KDFOR = KHFOR /24 
IHR6CKRHP) = KHFOR -24*KDF0R 
I FC I HRE(KRHP) .NE. 0 )G0 TO 103 
KDFOR = KDFOR - 1 

l\ 3 DAY = RHP DC KRHP) , * '.4?N 

1 ( 4 CALL DAYNXTCDAY, DPY) 

KDFOR = KDFOR - 1 ' / '' hr,} 

IF (KDFOR .GT. 0) GO TO 104 if 

I DYE (KRHP) = CAY . • 

GO TO 11*6 

145 I DYE ( KRHP ) = RHPDC KRHP) 

IHRFCKRHP) = RHPHC KRHP ) + IPTE 

106 CONTIMUh ' . • ; E ' ' : ■ 

C ELIMINATE HYDRO GRAPH OVERLAPPING 
NRHP1 * NRHP - 1 

; IFCMRHPl .PO. 0) GO TO 109 

1 DO 1C 8 KRHP = 1, NRHP1 . . -,. : m T|. 

I F( (IDYF ( KRHP ) . GT. IDYBCKRHP+1) .AND. (.NOT. C CIDYECKRHP) .G|. 

1 274 .AND. I CYB ( KRHP+1 ) .LE. 273) .OR. I DYE (KRHP) .EQ. 366)5) .OR 

2 (1DYRKRHP) *FQ. I DYB { KRHP+1 )* AND. IHRE (KRHP ) .GT. IHRBC KRHP+1 ) 

3 ) ) GO TO 101 . v * .. -0 ' * 

GO TO IDE 


SHM18520 

SHMI8530 

SHM18540 

SHM1855C 

SHM18560 

SHM18570 

SHM18580 

SHM1859C 

SHM18600 

SHM18610 

SHM18620 

SHM18630 

SHM18640 

SHMIM50 


*i 9 * * <4* w 1 » 1 

SHM 18680 
SHMt 8,690 
SHM18700 , 
SHM1871S 

SHM18720 
SHM 18730 
SHM1874T 
SHM 187:5# 
SHM18760 
SHMiSITft 
SHM 187 BO 
SHM1879D 
%mu%m ' 

S MM 1 SB TO 


SHM18830 

SHM18.840. 


IDYL (KRHP) 
IHRE(KRHP) 
CONTINUE. 
IFCIOYB(l) 


IDYBCKRHP+1) 
IHRBC KRHP+1) 


.AND. RHPDC1 ) .GE. 274 


RHPDC 1) 




SHWT8860 ; 
, SWM1887© 
SHM 18886 
366)SHM18890 




AND. IDYE(NRHP) 


1 GO 10 1,10 

GO TO 111 

IDYD(1) * 274 
IHRB(l) = 1 

I F{ I DYE( NRHP > *GE« 274 * AND. RHPD(NRHP) 0 LE„ 273 .AND. I DYE ( N! 
1 .HE. 366) GC TO 112 

GO TO 113 

1 DYE (NRHP) * 273 
I HRE (NRHP ) * 24 

CONTINUE 

DO 114 KRHP * It NRHP 

WRI TE( 6» 1 ) KR HP, IDYB ( KRHP ) » I HRBIKRHP ) , I DYE (KRHP) , IHRE ( KRHP ) 
FORMAT (5X, 'RUNOFF CONTRIBUTING TO RECORDED HYDRO GRAPH* *12/ 10X 
1 • BEGINS ON CAY*. 14,* AT HOUR' ,I3/10X, * AND ENDS ON DAY* .14. 


HOUR' 


CONTINUE 

RETURN 


AT HOUR' ,I3/10X, * AND ENDS 


SUBROUTINE TI MERTC SSR, SRR,CTRI , NCTRI , KRHP, KPCH ) 

PERFORMS CHANNEL TIME ROUTING 

DIMENSION SSR (5, 170 ) ,SRR( 5 ,170) ,CTRI(99) 

DO 100 KHPT * 1, KPCH 

10 $RR( KRHP, KHPT ) « 0.0 

■ ■ - ; : • • 

11 KdMTWUr : >. ' 

OO 102 KHPT - KTRI , KPCH 
* NRTRI » KHPT -KTRI + 1 


n'*l 

DAY* 


NE DA YNXT(DAY»DPY) 

NUMBER OF NEXT DAY OF 
DAY,DFY 


DAY = 1 
, AND. DRY 
DAY=60 


YEAR 


366) DAY = 


SHM1890Q 
SHM1891L 
SHM18920 
SHM1893I * 
HP) SHM 18940 
SHM1895* 
SHM18960 
SHM18970 
SHM 18 980 
SHM18990 
SHM 19000 
SHM1901C 
SHM 19020 
SHM 19030 
SHM 19040 
SHM19050 
SHM1906C 
SHM 190 70 
SHM19080 
SHM 1909© 
SHM19106 
SHM 191 Hi 
SHM19120 
SHM19130 
SHM19140 
SHM 191 50 

wmnm 

SHM 191 TO 
SHM19180 
SHM 191 90 
SHM 19200 
SHM19210 
SHH19220 
SHM1923C 
SHM19240 
SHM 192 50 
SHM1926C 


102 SRR( KRHP, KHPT ) « CTRKKTRI )*SSR( KRHP, NRTRI ) + SRR( KRHP, KHPT ) SHM19170 

KTRI * KTRI + 1 "-v.-: . . SHM19180 

I F ( KTR I .LE. NCTRI) GO TO ' 10 1 7 

RETURN ' • ■ , t#Wt©2©0 

iimd ■ . ■ 

SUBROUTINE DA YNXT (DA Y, DPY) • « - , , . . ‘>§IWt©lS0’ ' 

C DETERMINES NUMBER OF NEXT DAY OF THE YEAR ' ’ ■’ ;' : SHMi : 92»0 

INTEGER DAY, DFY ' • ■■iV .'$HN192©0. 

DAY = DAY + 1 v SHM19250 

IF (DAY .EQ. 366) DAY * 1 ' ; ' ;.§©f»l©i0® 

I F ( DAY ,r*Q. 60 .AND. DPY .EQ. 366) DAY = 366 - '.O'.: ■ " ■; 

I F ( DAY oEQ.367) DAY=60 . /$M©t9i©0 / 

RETURN rS'HM19290’ 

HMD l| $H«M©00 

SUBROUTINE £ VPDA Y( DPET , EMAET ) ' SHM19310 

C DETERMINES DATED PAN EVAPORATION TOTALS SHM1932C 

DIMENSION DPET ( 366 ) SHM19330 

RETURN SHM 19340 

END SHM1935© 

; SUBROUTINE DA YQUT( VDCY , MFDWY ,DPY ) - SHM 19' 6, 

C PRINTS TABLE OF DAILY VALUES SHM19S7© 

i; DIMENSION MEOVY( 12 ) , VDCYl 366 ) , VDMD( 12 ) SHM19380 

) : IT . G K DATE , CAY , DPY SHM 1 9391 

| 1C' writ r ,1) ‘ . 

I 1 FORMAT ( 7X» 3H0 AY, 7X »3HJUN» 5X» 3HJUL, 5X» 3HAUG»5X»4HS£PT , 5X, 3HQCT , 5X, SHM19410 
1 3HN0V, 5X »3HCEC , 5X, 3H J AN, 5X »3HFEB»5X»3HMAR, 5X,3HAPR»SX ,3HHAY ) SHM 19*20 

I MEDWYI3)-': $HMI943 : M 

K •mmmm y * \ , r , , 


DO 1 4 DATE = 1,28,1 

1 F ( MOD (DATE, 5 ) .NE. 1) 60 TO 102 
DO ID! KMO = 1,12 
DAY = MEDWY{ KMO) + DATE 
1C 1 VDMD(KMO) » V CCY (DAY ) 

WRITE (6,2 ) DATE, VDMD ( 12} , ( VDMD(KWD) , KWD=1, 11) 

2 FORMAT <1H9,3X, 16, 3X,12F8.1> 

GO TO IDA- 

102 DO 103 KMU»i, 12 
0AY-MEDWY1KM0 ) + DATE 

103 VDMD(KMO) * VCCY(DAY) 

WRITE<6,3) DATE, VDMD (12 ) , ( VDMD(KWO) » KWD*1, 11 ) 

3 FORMAT { IX , 3X , 16, 3X, 12F8. 1 ) 

104 CONTINUE 

I F ( DPY .HE. 366) GO TO 106 

DATE « 29 

VDCY (60 ) « VDCY( 366 ) 

00 105 KMO * 1,12 

DAY » MEDWY(KMO) + DATE 

105 VDMD (KMO) * V CCY (DAY) 

WRITE! 6,3 ) DATE, VDMD (12), (VDMD (KWD), KWD=1,11) 

GO TO 107 
1Q6 CONTINUE ' • . j, i • 

-Nf|iti*6»4) VDCY( 302) ,VDCY{ 333) , VDCY (363) ,VDCYt 29) , VDCY( 88) , 

1 VOCY (119) , VDC Y( 149 ) , VDC Y( 1 80 ) ,VDCY ( 210 ) , VDCY ( 241 ), VDCY ( 272 ) 

4 F0RMAT(1X,7X,2H29,3X,4F8.1,8X,7F8.1) 

107 CONTINUE 

108 WRITE! 6, 5) VDCY( 303 ), VDCY ( 334) , VDCY (364) , VDCY ( 30 ) , VDCY ( 89) , 

1 VDCY (123) , VDCY (150), VDCY( 181), VDCY (211), VDCYt 242), VDCY( 273) 

5 PHRMAT (1X»7X» 2H30 ,3X ,4F8. 1 ,8X,7F8. 1 ) 

WRITE (6, 6) VDCY( 304) , VDCYt 365) , VDCY (31 ), VDCY (90 ), VDC Y ( 151 ) , 

1 VDCY (212) »VDCY(243) 

6 FORMAT (lH/,7X,2H31,3X,F8«l,8X,2F8wl,8X,F8el,8X,F8*l,8X,2F8«l) 

MEDWYC3) * 365 ' lipl 

RETURN ■' • ' '■ fill 

cmd ' 

SUBROUTINE PR EPRD( RGPM, DRHP, DAY, HOUR, DPY, PRO, PEP *PRH) 

; DIVIDES HOURLY PRECIPITATION TOTALS AMONG PERIODS FOR SMALL BASINS 
DIMENSION PP4F(4) , DRHP( 366 , 24) 

INTEGER DAY, DPY, HOUR, PRD 

PFP=0.0 ' ■ ' .v ' 

1F(PRH .HQ. 0.0) RETURN 

IF (PP.D .EU. 1) GO TO IOC ■ - ,J ? ‘ 

PPP = Pi 4P(PRC) 

RETURN 

ICC LHOUR = HOUR- 
LOAY = DAY 
I F( LHOUR « GE« 

LHOUR * 24 
LDAY»D/.Y-1 
IF( LOAY • 6Q- 
[ 1 F ( LDAY .EQ« 

IF( LDAY . HQ. 




1 


1 ) GO TO 101 


C) LDAY*36S 
365) LDAY=59 
59 .AND. DPY 


EQ. 366) LDAY=366 


101 PRLH = kGPM*DRHP (LDAY, LHQUR) 


SHM 19440 

SHM19450 

SHM19460 

SHM 19470 

SHM19480 

SHM1949C 

SHM1950D 

SHM19511 

SHM19320 

SHM19530 

SH Ml 9540 

SHM1955C 

SHM1956C 

SHM 19570 

SHM19580 

SHM19590 

SHM19600 

SHM 190 115 

SHM1962U 

SHM 19630 

SHM19640 

SHM19650 

SHM 19660 

SHM19670 

SHM19680 

SHM19690 

smimrn 

SHM 19710 
SHM19720 
SHM' 19730 
‘SHM 197 40 
SHM19750 
SHM 1 1,9760 
SHM1977C 
.SHM2.97.dO 
SHM19790 
SH Ml 9800 
SHM19810 
SHM19S2-: 
SHM19S3 
SH Ml 9840 
SHM198J0 
SHM19860 
SHM19870 
‘ SHM 19880 
SHM19890 
SHM199-;,'. 

SHM 19910 

SHM 199 
SHM19930 
SHM1994< 
SH Ml 9950 
SHM 19960 
SHM19970 


NHOUR = HOUR +■ 1 
NDAY » DAY 

I F ( NHOUR .LE. 24) GO TO 102 

NHOUR* 1 

CALL DAYNXT { NCAY »DPY ) 

1C 2 PRNH *R GP M*DRH F(NDAY, NHOUR ) 

I F { PRH * GT* PRLH .AND. PRH . GT. PRNH) GO TO 103 

GO TO m 
103 PC4P (1) » 0.10 
P r :4P { 2 ) * 0.28 

PF4P (3 ) * ft. 46 
PE4F (4) * ft. 16 
GO TO 108 

1C 4 I F ( PRH .LI. PRLH .AND. PRH .LT. PRNH) GO TO 105 

GO TO 106 

105 PE4P { 1 ) * ft* 28 
PF4P (2 ) * 0* 10 
PE4PO) - 
PE4P(4) « 0.46 
GO TO 108 

106 IF(PRNH *GE. PRLH) GO TO 107 

PF;4P(l) * 0.46 

pup hi * o.i8 

PE4Pt3)- 0.2e 

fltfWW * 0.10 

GO TO .108 

107 PE4PU) = 0.1C 
PE4P (2 ) * 0.28 
PF4P (3 ) * 0.16 
PF4PC4) * 0.46 

108 DO 1C 9 KPRD = li4 

109 PE4P (KPRD ) * PE4PC KPRD ) *PRH 
PEP=PE4P { 1) 

RETURN 
END 


SHM199SI 

SHM1999 

SHM 20001 

SHM2001& 

SHM2002L 

SHM20030 

SHM20040 

SHM20050 

SHM 20060 

SHM20070 

SHM2G080 

SHH20O9O 

SHM2D100 

SHM20110 


SHH20120 
SHM20130 
SHM20140 
SHM20150 
' SHM2O160 
SHM2G17G 
SHH20180 
SHM20190 
SHM202Q0 
SHM20210 
SHM2O220 
SHM2-23. 
SHM2024C- 
■ SHM2025* 
SHM20 26 i 



